Quantum-critical theory of the spin-fermion model and its application to cuprates. 

Normal state analysis. 



O 

o 



O 
(N 

O' 
o ; 

u • 

Or 

-*-* . 



-a 
c 
o 
o 



o 



c3 



i 

c 
o 
o 



X 



Ar. Abanov 1 , Andrey V. Chubukov 1 , and J. Schmalian 2 

1 Department of Physics, University of Wisconsin, Madison, WI 53706 

2 Department of Physics and Ames Laboratory, Iowa State University, Ames, I A 50011 

(February 6, 2008) 



We present the full analysis of the normal state of the spin- 
fermion model near the antiferromagnetic instability in two 
dimensions. This model describes low-energy fermions inter- 
acting with their own collective spin fluctuations, which soften 
at the antiferromagnetic transition. We argue that in 2D, the 
system has two typical energies - an effective spin-fermion 
interaction g and an energy w 3 / below which the system be- 
haves as a Fermi liquid. The ratio of the two determines the 
dimensionless coupling constant for spin-fermion interaction 
A 2 oc g/uj s f. We show that A scales with the spin correlation 
length and diverges at criticality. This divergence implies that 
the conventional perturbative expansion breaks down. We de- 
veloped a novel approach to the problem - the expansion in 
either the inverse number of hot spots in the Brillouin zone, 
or the inverse number of fermionic flavors - which allowed us 
to explicitly account for all terms which diverge as powers of 
A, and treat the remaining, 0(A) terms in the RG formalism. 
We applied this technique to study the properties of the spin- 
fermion model in various frequency and temperature regimes. 
We present the results for the fermionic spectral function, spin 
susceptibility, optical conductivity and other observables. We 
compare our results in detail with the normal state data for 
cuprates, and argue that the spin-fermion model is capable 
to explain the anomalous normal state properties of high T c 
materials. We also discuss the non -applicability of the con- 
ventional </> 4 theory of the quantum-critical behavior in 2D. 



I. INTRODUCTION 

The 15 years after the discovery of superconductivity 
in cuprate oxides O witnessed a large number of efforts 
to understand the mechanism of superconductivity and 
numerous unusual normal state properties of cuprates. 
Parent compounds of cuprates (e.g., YBiC^O^) are insu- 
lators and also best examples of 2D Heisenberg antiferro- 
magnets pi. Upon hole doping, the antiferromagnetism 
disappears, and the system instead becomes a supercon- 
ductor. Above optimal doping (which is defined such 
that superconducting T c is the largest), superconducting 
features, such as the gap, A, in the spectral function, dis- 
appear at T c ||. Below optimal doping, the gap survives 
above T c , and disappears only at a higher temperature 
T* |J. In addition, T c and A(T = 0) become progres- 
sively uncorrelated with underdoping - T c decreases as 
doping decreases, while A(T = 0) gets larger [||(|. This 
clearly contradicts BSC theory in which A(T = 0) scales 
with T c . The region between T c and T* is called the 



pseudogap region. 
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FIG. 1. A schematic phase diagram of cuprates. 

A schematic phase diagram of cuprates is presented in 
Fig.pl. There is general agreement among the researchers 
that the superconducting state predominantly has d x 2_ y 2 
symmetry J7j. In this pairing state, the superconducting 
gap has nodes for diagonal directions along the Fermi sur- 
face, i.e., for k x = ±k y . Much less, however, researches 
agree at the present time about what is the origin of the 
pseudogap and also the anomalous normal state prop- 
erties of cuprates. These two features (the pseudogap 
and the anomalous normal state behavior) are challeng- 
ing for researchers and keep high T c problem afloat over 
such large period of time. 

This paper is devoted to the normal state properties 
of cuprates, so to keep our discussion focused we discuss 
only the experimental evidence for anomalous normal 
state behavior. The key evidence collaborated by very 
detailed photoemission ||[|, optical jl(|, transport |l|] 
and other studies is that near optimal doping, the in- 
verse quasiparticle lifetime is nearly linear in both tem- 
perature and frequency over a range of T between T c 
and lOOOX, and at frequencies roughly between few tens 
and few hundred meV. The Fermi liquid theory on the 
other hand predicts that the inverse quasiparticle life- 
time is quadratic in u> and T at the lowest energies. In 
overdoped cuprates, the low- frequency behavior of the 



quasiparticle lifetime is more consistent with the Fermi 
liquid theory Q. 

While there is no experimental proof at the moment 
that the measured non-Fermi liquid behavior at optimal 
doping extends down to u), T = (i.e., that there is no 
crossover to the Fermi liquid behavior at the lowest ener- 
gies), the very fact that the system behavior is markedly 
different from that in a Fermi liquid over a wide range of 
energies implies that theoretical description of cuprates 
should necessary involve strong coupling effect. 

From theoretical perspective, there exist two very dis- 
tinct proposals about what may cause non-Fermi-liquid 
behavior. First proposal is that this behavior is generic to 
doped Mott insulators in two dimensions and is due to the 
fact that strong Coulomb interaction not only makes the 
picture of weakly interacting electrons invalid, but elec- 
tron even ceases to exist as a quasiparticle (its residue 
is identically zero), and should instead be viewed as a 
convolution of two other objects one of which (spinon) 
carry spin but no charge, while the other (holon) carry 
charge but no spin [f!3[ . The theories based on the idea 
that spin and charge excitations in cuprates are separated 
have been worked out in great detail |1J] and substan- 
tially deepened our understanding of the physics of very 
strongly coupled electrons. 

A second proposal is based on the assumption that the 
anomalous normal state behavior is caused by the close- 
ness to a quantum phase transition of one type or an- 
other |la,nfl,n8 17 20]. Near phase transitions, quantum 



fluctuations are generically enhanced, and this should 
give rise to deviations from Fermi- liquid behavior pl| . 
Moreover, in low D systems, the region of Fermi liquid 
behavior progressively shrinks to lower and lower ener- 
gies as the system approaches quantum criticality. In 
distinction to the first approach, the idea that normal 
state behavior is due to a near quantum criticality neces- 
sary implies that at least on one side of a transition, the 
system behaves as a Fermi liquid at the lowest frequen- 
cies. 

The idea of quantum criticality emerged in the early 
days of high T c era as a way to justify the successful 
phenomenology |15| which used the linear energy depen- 
dence of the quasiparticle scattering rate as an input (the 
Marginal Fermi Liquid theory (MFL)). Another early 
idea of quantum criticality was focused on magnetic prop- 
erties of cuprates, and was based on the assumption that 
the behavior of weakly doped antiferromagnet is close 
to that of undoped systems near the quantum transition 
between magnetically ordered and disordered states P2| . 

At present, there are three suggestions as to what kind 
of quantum criticality may cause the observed non Fermi- 
liquid behavior. The first is a phenomenological idea that 
there exists a quantum critical point somewhere around 
optimal doping, which separates Fermi-liquid and pseu- 
dogap phases p3] . In this phenomenology, superconduct- 
ing region emerges as a dome on top of this point. The 



second, also phenomenological, suggestion is that quan- 
tum criticality and pseudogap behavior may be caused 
by a competition between superconductivity and some 
other yet undiscovered ordering (specific suggestion in- 
volves d— density wave 24 1). The third suggestion is 
that non Fermi liquid physics is due to the closeness to 
a magnetic quantum critical point which separates para- 
magnetic and antiferromagnetic phases p7|p(| . 

In this paper we explore the third possibility and 
discuss what kind of critical behavior one can expect 
near antiferromagnetic transition. The magnetic sce- 
nario emerged from the early days of high T c era |L7 25 
and was primarily based on the experimental observa- 
tions that although antiferromagnetism observed in par- 
ent compounds quickly disappears upon doping, short- 
range antiferromagnetic correlations persist up to much 
larger doping concentrations as is evident from NMR |27]] 
and neutron scattering (2^j2^] data, and these concentra- 
tions possibly cover the whole doping range of a super- 
conducting behavior. Another favorable feature of a mag- 
netic scenario, known even before high T c , is that near 
half-filling, the spin fluctuation exchange gives rise to an 
attraction in a d x 2_ y 2— pairing channel p5[ and hence 
can cause d— wave superconductivity p9| ] . A simple way 
to understand this is to assume that spin-fluctuations 
in high-T c materials play the same role as phonons in 
ordinary superconductors and create an analog of a de- 
formation potential in the form 



Tisf =9 1 dkdq c: ka a a ^c k+q ,pS. 



-<r 



(1) 



where a l a a are Pauli matrices and g is the coupling con- 
stant. Using this Hamiltonian, one can obtain an ef- 
fective pairing interaction between fermions which am- 
plitude is proportional to the propagator of the bosonic 
field S q , i.e., to the dynamical spin susceptibility. In a 
weak coupling BSC theory, spin dynamics is irrelevant, 
and the gap equation in a spin singlet channel has the 
form ||30tl 



A fc = -|g 2 / x(k-p)A p 



tanh(e p /2T) 
2Z 



dp 



(2) 



In distinction to phonons, the r.h.s. of this equation has 
an extra minus sign due to a summation over spin com- 
ponents |29f| (a projection of a^po^s onto a singlet spin 
channel yields — {2>/2){8 a p8 1 s — SasSp-y) whereas a con- 
ventional spin-independent interaction yields no overall 
minus sign). Because of a sign change, an s— wave so- 
lution A/j = const is impossible. However, since x(*l) 
is peaked near an antiferromagnetic Q = (7r,7r), pairing 
interaction relates the gap at momenta k and k + Q. In 
this situation, one can eliminate the overall minus sign 
by using an ansatz A k = — A^+q. For tetragonal lat- 
tice, this ansatz implies d x i_ y i symmetry of the pairing 
gap pi], i.e., spin fluctuation exchange gives rise to a 
d— wave superconductivity. 
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FIG. 2. Two different theoretical paths to study the 
physics of underdoped and optimally doped cuprates. In the 
first path, labeled as "theory I" , the point of departure is the 
antiferromagnetic state at half-filling. In the second path, la- 
beled as "theory II" the point of departure is the conventional 
metallic state in the overdoped regime. 

The studies of the normal state properties of cuprates 
within the magnetic scenario could also be divided into 
two subclasses (see Fig. ||). In the first class, which we 
labeled as "theory I" , the point of departure is Heisen- 
bcrg antiferromagnetism at half-filling, and the issue ad- 
dressed is how the system evolves with increasing dop- 
ing (I|,^9|,||j32H],||,||,|6|. These studies mostly con- 
sider small concentrations of electrons at which either 
long-range or quasi-long range magnetic order keeps ma- 
jority of electrons localized, and the Fermi surface con- 
sists of separated patches (hole pockets) . In the second 
class, labeled in Fig g as "theory II" , the point of de- 
parture is the Fermi liquid metallic state at large dop- 
ings, and the issue addressed is how this Fermi liquid 
behavior modifies and eventually disappears as the sys- 
tem approaches the magnetic transition. In this second 
approach, the physics is governed by the interaction be- 
tween electrons and their collective spin bosonic excita- 
tions which become soft modes near the transition. This 
model is often called the spin-fermion model. Indeed, 
both approaches describe the same strongly interacting 
electronic system near magnetic instability and should 
yield the same results for the same regions of doping and 
temperatures, unless there are instabilities at large en- 
ergies (comparable to e.g., fermionic bandwidth), which 
cannot be detected in the low-energy theory and have to 
be taken as inputs. From this perspective, the choice of 
the model is somewhat subjective and is mostly deter- 
mined by which starting point is thought to be closer to 
real situation at a given doping. Very near half-filling, the 
system clearly displays Heisenberg antiferromagnetism, 
and the Fermi surface should consist of hole pockets. In 



this limit, one has to work rather hard within the spin- 
fermion model to reproduce the features which are readily 
obtained if one departs from antiferromagnetic state at 
half-filling. On the other hand, we believe that the spin- 
fermion model is a better starting point if one wants to 
describe the physics at and somewhat below optimal dop- 
ing mostly because at present the majority of photoemis- 
sion experiments even in heavily underdoped cuprates do 
not show the emergence of the shadow Fermi surface (the 
extra piece of the Fermi surface shifted from the original 
one by antiferromagnetic Q), which is the necessary in- 



gredient for the formation of hole pockets 1 38 1 . 

The present paper summarizes our analysis of the spin- 
fermion model in the normal state, or, more precisely, 
when the pairing correlations can be neglected. We ar- 
gue that in 2D, a dimensionless coupling constant for the 
spin-fermion interaction A scales as inverse magnetic cor- 
relation length £ and diverges as the system approaches 
a quantum-critical point (QCP) of a transition to a mag- 
netically ordered state. We show that due to a divergence 
of A, there exists a wide region in the T, A plane where the 
system is in the quantum-critical regime near the mag- 
netic transition. We will show that in this regime, the 
system behavior is qualitatively different from that in a 
Fermi liquid. 




T 


b) 


\ Quantum 
\ critical 




afm\ /'" 


FL 



doping 



doping 



FIG. 3. A cartoon of the system behavior near antiferro- 
magnetic quantum-critical point 

Our point of departure is the "zero-order" phase di- 
agram with antiferromagnetism included as an input ( 
Fig.||a). The deviations from the magnetic transition 
are measured by the doping dependent spin correlation 
length £(x). We assume that all other features on the 
phase diagram, such as superconductivity, pseudogap 
regime, and the crossover from a Fermi-liquid to not 
Fermi liquid behavior in the normal state are low- energy 
phenomena which should be obtained within the model. 
Alternatively speaking, we assume that antiferromag- 
netism is produced by fermions with energies comparable 
to the bandwidth, while other phenomena are produced 
by fermions near the Fermi surface. In cuprates, the sep- 
aration of scales is not very strong (T* for strongly under- 
doped materials is a fraction of an exchange integral J), 
i.e., to some extent high-energy fermions do contribute 
to pseudogap phenomenon and may also contribute to 
anomalous normal state behavior. This contribution is 



inevitably lost in our approach. However, we believe that 
the fundamental physics can be understood by separat- 
ing the scales perhaps more strongly than in reality, and 
the inclusion of fermions with energies comparable to the 
bandwidth will change the results quantitatively but not 
qualitatively. 

The phase diagram which emerges from our studies is 
presented in Fig. rab. We argue that the region near 
the QCP is divided into a Fermi liquid regime and a 
quantum-critical regime where the system behavior is the 
same as at the critical point. The upper boundary of 
the quantum-critical behavior is roughly located at fre- 
quencies comparable to spin-fermion coupling constant g. 
The crossover from a Fermi liquid to a quantum-critical 
behavior on the other hand occurs at energies of order 
uj s f ~ gX 2 where A, which we already introduced above, 
is the dimensionless coupling constant in the problem 
(A ex g/vp£ , where vp is the Fermi velocity). At weak 
coupling, ui s f exceeds g, and the quantum-critical behav- 
ior is not realized. However, at strong coupling, to s f <C g, 
and the system displays quantum-critical behavior in a 
wide range of frequencies. This strong coupling behav- 
ior necessary occurs very near the transition as A ex £, 
but it can also be reached at intermediate £ when the 
spin-fermion interaction g increases. 

In the Fermi liquid regime, fermionic self-energy be- 
haves as S" ex to 2 /io s f. In the opposite limit, when 
uj s f = 0, we found that it scales, up to logarithmical 
corrections, as E" ex y/gE, where E — max^T). At 
a finite Lo s f, we find that there also exists a wide inter- 
mediate region, roughly between 0.5w s / and 6 — 8w s /, 
where E" ex E, i.e., it is linear in both frequency and 
temperature. This behavior of the self-energy causes a 
cascade of crossovers in the fermionic spectral function 
and conductivity measured as functions of temperature 
and frequency. 

These are the key results of the paper. We derive them 
in Sec. |l| -jvj below in a formal l/N expansion where 
N can be regarded either as the number of hot spots 
in the Brillouin zone (crossing points between the Fermi 
surface and the magnetic zone boundary), or the number 
of electron flavors. 
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FIG. 4. A schematic phase diagram of cuprates obtained 
in fl3£(] in the Eliashb erg- type formalism 

As we said, in this paper we focus on the normal state 
properties of the spin-fermion model. For completeness, 
we also present in Fig. ^ the phase diagram which we ob- 
tained |3£] in the Eliashbcrg-type formalism (extended 
compared to that for phonons J42J to include the feed- 
back from pairing on magnetically mediated interaction) 
by using the normal state results as a base and adding 
the interaction in the pairing channel. Out key point 
here is that the two scales u> s f and g which we find in 
the normal state are also present in the pairing problem. 
Namely, the pairing instability temperature Ti ns tends to 
a finite value ~ g (more accurately, g/N) at the quantum- 
critical point. In a conventional situation, this would 
imply that superconducting properties gradually evolve 
below Ti ns . ln our situation, this, however, is not the case 
as at strong coupling, when Ti ns (xj> w s /, the pairing 
comes from frequencies which well exceed uj s f, and thus 
involves incoherent fermions [E0[. This pairing is quali- 
tatively different from that in the BCS theory. We found 
that in this situation, immediately below Ti„ s electrons 
form singlet pairs, but the pairs still behave incoherently 
and do not propagate. Only when the temperature is 
reduced well below Ti ns , the feedback from superconduc- 
tivity eventually "clears up" fermionic excitations, and 
the system behaves as a conventional superconductor. 
We found that the conventional superconducting behav- 
ior with e.g., a sharp quasiparticle peak in the spectral 
function) emerges only when the temperature becomes 
smaller than the fraction of {ujsfg) 1 ' 2 ~ g/X <C g. Phys- 
ically, this second scale (modified u) s f) is the gap in the 
spectrum collective spin excitations which become prop- 
agating, magnon-like in a d— wave superconductor |41j| . 

The cuprates are not the only objects which motivate 
our study. Over the last few years, there is a growing 
interest in the behavior of various heavy-fermion materi- 



als near the magnetic instability which can be achieved 
either by chemical variation (i.e., doping) or by applying 
external perturbation such as pressure pjjE4LE7|. The 
magnetic instability can be either ferromagnetic or an- 
tiferromagnetic, depending on material. There is a con- 
vincing experimental evidence that at least in some of 
heavy fermion materials, there exists a dome on top of a 
quantum critical point where the system is in a supercon- 
ducting phase Q . From this perspective the behavior of 
heavy-fermion materials may be not very different from 
that of cuprates. There is also a number of experimen- 
tal data which show that the normal state behavior near 
magnetic instability likely deviates from that in a Fermi 
liquid down to very low energies Jl5|,|46| . These deviations 
have been observed in both fermionic and spin proper- 
ties [jj7f . We caution, however, that heavy-fermion mate- 
rials are mostly three-dimensional |46[| while our analysis 
is valid in two dimensions. Still, we believe that the com- 
putational scheme we developed for 2D can be applied 
also to 3D materials at strong coupling. 

Our approach to the spin-fermion model is an exten- 
sion of the earlier detailed studies by D. Pines and co- 
workers [ pTp9| , |4q | . The only qualitative difference is that 
we do not assume a'priori the form of the dynamical sus- 
ceptibility but rather derive it. From this perspective, 
our work should be regarded as an attempt to put the 
studies of the spin-fermion model on a solid theoretical 
basis and to understand whether one can attack the prob- 
lem analytically, in a controllable way. Previous studies 
of the spin-fermion model were mostly performed numer- 
ically. 

Also, our approach and computational technique are 
in many respects similar to the studies by B. Altshuler, 
A. Millis and L. Ioffe fl4J3,|5(| of fermions interacting with 
low-energy collective modes. Our studies also agree, in 
most part, with numerical studies of the spin-fermion 
model by Benemman and co-authors |pl| . Our analysis 
of the effective bosonic theory of the model at critical- 
ity agrees in most parts with the studies of Lercher and 
Wheatley @. 

Finally, our analysis bears some parallels with large 
D studies of fermionic systems [[34]]. Just as in D — ► oo 
theories, our self-energy predominantly depends on fre- 
quency ( at the lowest energies the renormalization of e& 
in the fermionic propagator is by log A/A, smaller that the 
renormalization of lu). The difference, as we understand 
it, is that large D theories were chiefly applied to explain 
the effects, such as metal-insulator transition, for which 
fermions with energies comparable to the bandwidth play 
an important role. We, on the contrary, focus on univer- 
sal features of the system behavior at the lowest frequen- 
cies, when the fermionic density of states can be approx- 
imated by a constant, i.e., the bandwidth can be set to 
infinity. The paper is organized as follows. In the next 
section, we discuss the spin-fermion model and its rela- 
tion to lattice models of cuprates. In Sec |n| we present 



the calculations at T — 0. We first discuss ordinary per- 
turbation theory and then 1/N expansion in the strong 
coupling (A ^> I) limit. We then proceed in Sec. fV\ with 
the calculations at N — ► oo. In the next Sec |y| we dis- 
cuss 1/N expansion at T = 0. In this Section we also 
compare our results at the QCP with the conventional 
4> 4 theory of such transition |5^,|56| . We argue that spin- 
fermion model at criticality is not described by a standard 
<^> 4 theory because of complex momentum and frequency 
dependences of the effective four-boson vertex made out 
of fermions. In the next Sec VI we present the results 
of our calculations at finite temperature. Here we ob- 
tain the full form of the fermionic spectral function which 
can be d irectly compared to the photoemission data. In 
Sec. VII we discus in detail f the behavior of optical con- 
ductivity a = o\ + io~2, effective scattering rate 1/r and 
the effective quasiparticle mass. These expressions can 
be directly compared with optical and transport data. 
In Sec VIII we present the results for other observables: 
fermionic density of states, dynamical spin susceptibility, 
and Raman intensity. Finally, in Sec IX we discuss the 
behavior of the pairing vertex and show that although 
pairing fluctuations in the normal state generally affect 
fermionic propagator only to order 1/N, the ladder cor- 
rections to the pairing vertex are not small in 1/N. We 
argue that the ladder series give rise to a d— wave pairing 
instability. The results of our analysis are summarized in 
Sec fxl In this last section we also compare our results for 
the fermionic spectral function and conductivity with the 
photoemission and optical data for near optimally doped 
cuprates. 

Some of the results presented in the paper have already 
been discussed in short communications 1 20 3^{I(lp7] p2| . 



II. SPIN-FERMION MODEL 

In this section, we present the microscopic justifica- 
tion for the spin-fermion model. We use renormalization 
group (RG) arguments and general considerations to ob- 
tain the effective Hamiltonian which describes the inter- 
action between low-energy fermions and their collective 
spin degrees of freedom. We then develop a new method 
to compute the forms of the fermionic propagator and 
dynamical spin susceptibility for strong spin-fermion in- 
teraction. 

The first-principle approach to cuprates demands 
that one chooses a microscopic model which describes 
fermions on a lattice with a short-range, Hubbard-type 
interaction. If the choice is correct and if one man- 
ages to solve this model and understand the result, there 
is no need to introduce any low-energy phenomenology. 
Unfortunately, this microscopic model is unknown and 
presumably is rather involved. The Hubbard model 
with on-site interaction is a undoubtly a reasonable first 



choice [p3| , p3| , but one has have in mind that it in- 
deed gives only an approximate description of the actual 
physics in cuprates. Furthermore, even for the single 
band Hubbard model no exact solution has been yet ob- 
tained in two dimensions, and one has to cither study the 
model numerically |63J33], or use an approximate theo- 
retical scheme such as small U expansion J64[ , or an ex- 
pansion in the artificially extended number of electronic 
orbitals per site J65[ . 

In the low-energy approach which we advocate, we are 
not interested in the behavior at high energies (compa- 
rable to the fermionic bandwidth W), and want to get 
only a limited information about the system, namely, how 
it behaves at energies smaller than some cutoff energy 
A < W. In the RG sense, this effective theory can be re- 
garded as obtained by integrating out fermionic degrees 
of freedom with energies between W and A. 
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It is generally believed that this procedure does not 
take the system away from the basin of attraction of the 
Fermi- liquid fixed point. We are not aware of any solid 
calculations which would indicate the opposite. In other 
words, we assume that the fermionic interactions at high 
energies do not destroy the Fermi-liquid behavior. This, 
however, does not imply the throughout applicability of 
the conventional Landau Fermi liquid theory because of 
possible infrared singularities. In fact we show below that 
in some range of parameters the fermionic self-energy cor- 
rections are indeed singular, and give rise to a non-Fermi 
liquid behavior. 

Strictly speaking, the RG procedure is justified only 
if the spin-fermion coupling does not exceed fermionic 
bandwidth which we assume to hold. If this condition 
is not satisfied, the separation between low-energy and 
high-energy excitations becomes problematic. We be- 
lieve that our results are valid, with some minor modifi- 
cations, also in the limit of very large couplings provided, 
indeed, that there is no destruction of the Fermi liquid 
due by fluctuations at lattice scales. However, universal 
results which we report here are only valid if the coupling 
are smaller than the bandwidth, and typical fermionic 
momenta are located near the Fermi surface, when the 
fermionic dispersion can be linearized. 

In general, the low-energy behavior of fermionic sys- 
tems is governed by degrees of freedom which have low- 
energy excitations. One such degree of freedom is obvi- 
ously given by the fermion itself since it possesses an ar- 
bitrary low energy near the Fermi surface. Potential can- 



didates for other low energy excitations are the bosonic 
collective modes of fermions. In a general situation, these 
collective modes are gaped already at energies compara- 
ble to the bandwidth. However, if the fermionic system is 
close to an instability of some kind, then the correspond- 
ing bosonic collective mode has a much smaller gap which 
vanishes at the instability point. The closeness to anti- 
ferromagnetism naturally makes spin excitations relevant 
candidates for such low-energy degrees of freedom. 
Throughout this paper we assume that 

1. spin fluctuations can be treated as low-energy ex- 
citations in the wide range of dopings, 

2. one can neglect the effects of all other low-energy 
collective degrees of freedom, independent of spin 
excitations 

These are two basic assumptions of our approach. 
Their applicability is solely determined by comparisons 
with the experimental data. We believe that the NMR 
and neutron scattering data in cuprates indicate that for 
all doping ranges (including overdoped materials), the 
spin susceptibility possesses a substantial momentum and 
frequency dependence near the antiferromagnetic mo- 
mentum (7r, 7r), and the typical energies associated with 
these dependencies are at least by order of magnitude 
smaller than the fermionic bandwidth. We also believe 
that there is no clear experimental indication of the pres- 
ence of the low-energy excitations in other interaction 
channels. In principle, this does not preclude moderately 
strong effects from other channels which can account for 
moderate values of the Landau parameters J66|]. These 
effects may change numbers but certainly do not change 
the basics physics. Note also that our approach does not 
rule out the effects in other channels which are secondary 
in the sense that they emerge within the spin-fermion 
model. For example, charge excitations may acquire a 
low-energy dynamics due to a coupling with spin fluc- 
tuations. This coupling may eventually lead to phase 
separation []67|J68| and, potentially, (when lattice effects 
are included) to the formation of stripes |69j . 

The treatment of spin fluctuations as separate bosonic 
degrees of freedom requires some extra clarifications. The 
point is that in distinction to heavy-fermion materials, 
there is only one spice of fermions in cuprates, and spin 
fluctuations just result from a multiple interactions be- 
tween particles and holes, i.e., are collective modes of the 
fermions. From this perspective, the introduction of spin 
excitations as an extra low-energy degree of freedom is 
just a convenient way to separate the energy scales: we 
assume that there exists a single dominant channel for 
fcrmion-fcrmion interaction at energies smaller than A, 
and introduce a (spin) collective mode which mediates 
this interaction. At the same time, the static propagator 
of this collective mode (and, in particular, the value of the 
spin correlation length) is obviously determined by high 



energy fermions about which we have little information. 
We only expect, as we already said, that the integration 
over high fermionic energies does not give rise to any sin- 
gularity in the bare spin susceptibility. The latter then 
should have a regular Ornstcin-Zcrnikc form 



xo(q.w) = 



Xo 



(3) 



£- 2 + (q-Q) 2 -(^M) 2 ' 

The input parameters in Eq. (m are the spin correla- 
tion length, £ and the spin velocity v s which is obviously 
of order vp as spins are made out of fermions. The over- 
all factor xo can be absorbed into the renormalization 
of the coupling constant (see below) and should not be 
counted as an extra variable. 

Note that the bare Xo(Qi w ) is rea l- This is an essential 
point in our consideration. We argue that the imaginary 
part of the susceptibility is determined by fermions with 
energies smaller than A and has to be computed within 
a low-energy theory rather than taken as an input. In- 
deed, since the spins have no source of damping other 
than to decay into a particle hole pair, the inverse life- 
time of a spin fluctuation coincides with the imaginary 
part of the fully renormalized particle-hole bubble. Due 
to the necessity to conserve energy and momentum in a 
decay process, it involves only low-energy fermions with 
frequencies smaller than the external spin frequency. 

In earlier phenomenological studies of the spin-fermion 
model, the full spin susceptibility, including its imaginary 
part, was considered as an input |Ffl . The spin dynam- 
ics was assumed to be a simple relaxational one with 
X _1 (Q,w) oc 1 — iuj/uj s f. The value ofu) s f was taken from 
experiments and assumed to vary with doping indepen- 
dently of other input parameters. The phenomenological 
assumption about the form of the full susceptibility is 
roughly consistent with our findings. At the same time, 
we will see that ui s f is expressed in terms of other input 
parameters and cannot be varied independently of the 
spin-fermion coupling. 

Note that the damping term in the spin susceptibil- 
ity only appears if the Fermi surface contains hot spots 
(points separated by Q) |7Q] . For a Fermi surface without 
hot spots the spin decay is forbidden at low frequencies 
due to energy constraint. In this situation, the full spin 
susceptibility is real at small frequencies and differs from 
(||) only due to effects of quantum criticality. On the 
other hand, for a Fermi surface with hot spots, the spin 
damping is permittable down to the lowest energies and 
overshadows the quadratic in uj term in the bare suscep- 
tibility which therefore becomes irrelevant. 

The topology of the Fermi surface in cuprates is non- 
universal and in principle can vary from one material to 
the other. ARPES data indicate that at least in the best 
studied Bi-based cuprates, the Fermi surface is centered 
at (it, 7r), i.e., it is an open electron Fermi surface PJnfl. 
Applying Luttinger theorem which states that the area 
of electron states, measured in units of the Brillouin zone 



area, equals to the density of electrons J72[ , one can im- 
mediately make sure that such Fermi surface necessarily 
crosses the magnetic Brillouin zone boundary and there- 
fore contains hot spots. As we just said, in this situation, 
the uj 2 term in the bare spin susceptibility is irrelevant, 
i.e., Xo(q> w ) can be approximated by its static part. This 
in turn implies implies that the spin-fermion model can 
be described by the Hamiltonian. This Hamiltonian in- 
volves fermions which live near the Fermi surface, spin 
fluctuations, and the interaction between the two degrees 
of freedom, and can generally be written as 



k,a q 

+9 Yl 4+q.a °"«,/9 Ck,/3 • S_ q . 
q,k,a,/3 



(4) 



Here, c£ is the fermionic creation operator for an elec- 
tron with crystal momentum k and spin projection a, at 
are the Pauli matrices, and g is the coupling constant 
which measures the strength of the interaction between 
fermionic spins and the collective spin degrees of freedom 
described by bosonic variables S q . 

Alternatively to the Hamiltonian in Eq. [|, we can also 
write down the effective action for the model 

S = -f dr f dT'Y^L(r)Go\KT-T')c k<T (T') 

+\( dr f dr'^Xo 1 (q)S q (r)-S. q (r') 
L Jo Jo 

+g f dr^s q (r).S_ q (r). (5) 

Jo q 

where Gq (k, r) = d T — vj,(k — kp) is the bare Fermionic 
propagator. In general, this representation is more ad- 
vantageous as it allows a time dependence of the bare 
spin propagator Xo(Qi t )- In our case, as we said, this 
time dependence is irrelevant, and one can therefore use 
a conventional Hamiltonian description. 




FIG. 5. A schematic picture of the Fermi surface. The 
notation h.s. stands for "hot spot" - the point on the Fermi 
surface separated by another point on the Fermi surface by the 
antiferromagnetic vector Q. The arrows denote the velocities 
at hot spots. 0o is the angle between the velocities at hot 
spots separated by Q. 

The input parameters in (0) are £, the fermionic veloc- 
ity Vfc and the spin-fcrmion coupling g. We demonstrate 
below that the physics associated with the closeness to 
an antiferromagnetic instability mostly involves fermions 
near hot spots. For Q = (it, it), there are 8 hot spots in 
the Brillouin zone with the same vf = Wk\- The direc- 
tion of Vfc at two hot spots separated by Q differs by an 
angle <pQ. This angle, however, should not be counted as 
a relevant input parameter as we will see that the physics 
in insensitive to the actual value of </>o as long as </>o 7^ it, 
i.e., there is no nesting at hot spots. The results for the 
nested Fermi surface are different, and we will not dis- 
cuss this special case in the paper. Experimentally, in 
cuprates the hot spots are located not far away from the 
corners of the Brillouin zone (i.e., from (0, it) and sym- 
metry related points). In this situation, <fio is obviously 
close to it/2. For practical purposes, it is more convenient 
to introduce, instead of vp and 4>o, the two velocities v x 
and v y via e k = v x k x + v y k y and e k+Q = -v x k x + v y k y , 
where k measure a deviation from a hot spot. Obviously, 
v f = v x + v y an d 0o = 2t&n~ 1 v x /v y . The limit (fto = tt 
corresponds to v y — 0. (see Fig. ||) 

The spin-fermion coupling constant g is a fully renor- 
malized irreducible vertex for the particle-hole channel. 
In the Born approximation, g equals to the Hubbard U. 
Beyond Born approximation g acquires a strong doping 
dependence and rapidly decreases with increasing dop- 



ing. For example, an RPA renormalization of g in the 
particle-particle channel yields g — U/(l + (U/t)F(x)) 
where F(x), subject to F(0) = 0, is an increasing func- 
tion of the doping, x. This form shows that g ~ U only at 
small doping while at larger x it progressively decreases 
down to g ~ t. This form of g is indeed only an ap- 
proximate one as there is no justification to restrict with 
one particular renormalization channel at high frequen- 
cies. We can only quite generally assume that g is some 
doping dependent coupling constant which increases as 
the system approaches half-filling. 

To summarize, the relevant input parameters for our 
model are £, vp, and the effective coupling g = g 2 Xo, 
which is a combination in which g and xo appear in per- 
turbation series. The first two parameters can be merged 
into an energy scale Vf£ . Out of the two energies, g 
and vp£~ l , one can construct an overall scale and a single 
dimensionless ratio 

A = 3.g/(4 7 n; F r 1 ) 

which, as we will see in the next section, determines the 
strength of both fermionic and bosonic self-energy cor- 
rections (the numerical factor in A is introduced for fur- 
ther convenience). Physically, A measures the ratio of 
the effective coupling constant and a fermionic energy 
at a typical fermionic \k — kp\ ~ £ , which sets the 
momentum range for spin-fcrmion coupling. When A is 
small, fermions are nearly decoupled from spin fluctua- 
tions and behave as an almost ideal Fermi gas. On the 
contrary, when A is large, the bare fermionic dispersion is 
almost completely overshadowed by the interaction, and 
this may give rise to a non-Fermi liquid behavior. 

The fact that a large number of experimental data 
for cuprates already near optimal doping differ from the 
predictions of the Fermi liquid theory indicates that if 
spin-fcrmion coupling a relevant mechanism for devia- 
tions from the Fermi-liquid theory, A should be large al- 
ready for optimally doped samples, and increase with 
decreasing doping. Our estimates of A using the pho- 
toemission §@0, neutron f|, and NMR data 
yields A ~ 2 near optimal doping (see Sec |x|) . This ob- 
servation is a challenge to a theory as real systems turns 
out to be outside the basin of applicability of a conven- 
tional, weak coupling perturbation theory and one should 
instead try unconventional procedures in a search for a 
strong-coupling solution. 

We now proceed with the evaluations of the fermionic 
and bosonic self-energies. 



III. PERTURBATION THEORY FOR THE 
SPIN-FERMION MODEL AT T = 

In this Section, we analyze the spin-fermion model at 
T = 0. Our strategy is the following: we first present the 



results of a formal perturbation theory for both fermionic 
and bosonic propagators. We show that perturbation ex- 
pansion for the fermionic self-energy and the spin-fermion 
vertex holds in powers of A and obviously breaks down for 
A > 1. We next show that at A > 1 the self-energy cor- 
rections to the spin propagator also become relevant and 
change the dynamics of spin fluctuations at frequencies 
which are mostly relevant for fermionic self-energy and 
vertex renormalization. We then demonstrate that even 
for large A, one can construct a variant of perturbation 
theory in which dominant, O(A) fluctuation contributions 
are included into new zero-order fermionic and bosonic 
propagators for which we obtain explicit answers. The 
remaining fluctuation corrections are logarithmical in A, 
and we analyze them in the one-loop RG formalism. 



A. A formal perturbation theory 

We begin with the direct zero-temperature perturba- 
tion expansion in powers of A. Direct perturbative ex- 
pansion means that all diagrams are evaluated using the 
bare values of the spin susceptibility and the fermionic 
propagator. 



1. fermionic self-energy 

We start with the fermionic self-energy, S(fc, O) re- 
lated to the fermionic propagator by G~ 1 (k,cu) = CI + 
S(fc, Q) — efc. Let's hrst compute £(fc, ft) near a hot spot 
expanding to first order in frequency and in the quasipar- 
ticle energy. We will demonstrate that this supposedly 
straightforward computation has to be performed with 
more care than one might expect. 



where 



q,oo 



k.Q 




FIG. 6. The lowest-order diagram for the fermionic 

self-energy. 

The lowest-order self-energy correction involves a sin- 
gle spin-fluctuation exchange and is presented in Fig. |6|. 
In the analytical form we have in Matsubara frequencies 



S(k,O m ) = -3. 9 2 J 



cPqdiUn 



(2tt) 3 
x G (k + q, uj m + fi m )xo(q, U) m ) . 



(6) 



G n 1 (k,U> m ) = lUJm - Vp(k— kfc s ). 



(7) 



Subtracting from (ft) the self-energy at fi m = and 
k = kfc s which can be absorbed into the renormalization 
of the chemical potential, and expanding to first order in 
k — kfts and in frequency, we obtain after simple manip- 
ulations 



S(k,O m ) = (iO m - e k+Q ) 7(k,f2 m ) 



(8) 



where 
I(k,Q m ) 



(2tt) 3 



d 2 qduj r , 
1 



£-2 + q2 + ( W / Us )2 

1 



id. 



e k +Q + iuj m - v F q x iu m - v F q x 



(9) 



Here q = q Q, and x-axis is chosen along vp at k + Q. 
In the spirit on a conventional expansion in O m and 
tk+Q, it is tempting to evaluate J(k, O m ) right at k = 
k/ ls and Q m = 0. The corresponding computations are 
performed in Appendix A, and the result is 



/(kfc„0) 



-A 



Vf 



(10) 



Substituting this into (g) one would obtain that the self- 
energy only accounts for the renormalization of the quasi- 
particle residue, and the magnitude of the renormaliza- 
tion depends on the ratio v s /vf- Apparently this is the 
whole story. However, it turns out that Eq. (Q) is not 
the full result, and 



lim /(k,O m )^J(k hsi 0) 

-+0.K — >kfo a 



(11) 



To show this we evaluate /(k, O m ) keeping both Q m and 
v^k in the integrand in (g). Analyzing the integrand, 
we observe that there exists a tiny region of frequencies 
where the poles in the two fermionic propagators are close 
to each other, but still are located in different halfs of a 
complex q x plane. For Ct m > 0, this region is sandwiched 
between u m > — fl m anduj m < 0. The integration over q x 
in (H) in this range of frequencies, and a subsequent inte- 
gration over Lo m yields an extra contribution to I(k, Cl m ) 
in the form 



lan\K-i *'m) — ^ 



ifl r 



0, 



£ fc+Q 



(12) 



We see that J an (k, i7 TO ) comes from the integration 
over internal frequencies u> m which are smaller than the 
external fl m . Clearly, this anomalous piece could not be 
obtained in a conventional perturbation expansion over 
n m as in the latter one assumes that the internal energies 
are much larger than the external one. This assumption 
is generally justified by a simple phase space argument, 
i.e., the contribution from \u) m \ < Vi m is normally small 



due to the smallness of the integration range. Here, how- 
ever, the smallness of the phase space is compensated 
by the smallness of the denominator in (pf) as without 
f2 m and ek+Q, the product of two Green's functions con- 
tains a double pole. At finite fi m and tk+Q, the double 
pole splits into two single poles, but the energy differ- 
ence between them is still only ifl m — ek+Q, i.e., is of the 
same order as the integration range over frequency. The 
same reasoning is used to extract the effects due to chiral 
anomaly in quantum chromodynamics | J77[ . 

Substituting the total expression for /(k, $7 m ) into the 
self-energy, we obtain 



the fermionic self-energy (—1 instead of of 3). Expand- 



E(k,fi) = A 






(iO m — £fc+Q)- 



V F 



(13) 



We see that both terms in E are of order A, i.e., there is 
no strong distinction between "anomalous" and "normal" 
contributions to the fermionic self-energy. This could 
be anticipated as, e.g., for v s = oo, Xo(q.jW) is purely 
static, and E(k, Q) should only depend on e^+Q as the 
dependence on external f2 m is eliminated by shifting the 
internal frequency in the fermionic propagator in m) . We 
indeed reproduce this result: for v s = oo Eq. (TF3) yields 
E(k,n m ) = Xe k+Q . 



2. vertex renormalization 

We next compute the lowest-order correction to the 
spin-fermion vertex for fermions at hot spots. 



k hs +q,co 




k hs +q+Q, co 



FIG. 7. The lowest-order diagram for the correction to 
the spin-fermion vertex. 

The lowest-order vertex correction diagram is pre- 
sented in Fig. [?]. The vertex correction is obviously the 
largest for a fermion at a hot spot. For k = k/j S and 
u> = we obtain g = g + Ag where 



9 

<hj m d 2 q 

3-x(q,w m j 



Ag . 
T =ty J'(2n 

x G{uj m ,^-hs + q) G(w. 



+ Q + Q) (14) 



Observe that the overall factor (the result of the sum- 
mation over spin components) is different from that in 



ing, as above, near hot spots as £k hs +q 



v x q x 



)%■> 



e k hs + q +Q = ~v x q x + v y q y where v 2 + v 2 = v F , and per- 
forming the computations described in Appendix A, we 
obtain 



A Vf 



A5 = = 

g 3 v s /3yl 



1 



In- 



(vT 



f>) 



a 2 ^/((3 2 +a 2 )(l-a 2 )-(3 



(15) 



where a — v y /v s and (3 — v x /v s . For purely static bare 
susceptibility, i.e., vf/v s — > 0, Ag/g is a smooth function 
of v T and v,, . 



Ag _ Xvf_ . , _i Vx_ 
g 3 v x ' v y 

In the opposite limit v s ^C vf, we obtain from ( p5| ) 

Ag A v s Vp 
9 



3 V X Vy 



. ± Vx 

tan — 

v„ 



(16) 



(17) 



Observe that the correction to the spin-fermion vertex 
vanishes when the velocity of spin excitations tends to 
zero. 

We also verified that the correction to the spin-fermion 
vertex does not contain a singularity, similar to what we 
found for a self-energy correction. In other words, the 
limiting value of Ag at vanishing frequency coincides with 
Ag evaluated right at zero frequency. The absence of a 
singular correction to a vertex is associated with the fact 
that the two internal fermions in the vertex correction 
diagram have different directions of the velocities, and 
hence there is no double pole which might give an extra 
piece after regularization. 



3. higher-order diagrams and the structure of the 
perturbation theory 

We see that the perturbative expansion for the 
fermionic self-energy and the spin-fermion vertex holds 
in the dimensionless parameter A. Higher-order diagrams 
can be easily estimated. They all scale as higher powers 
of A. To this end, diagrammatic, perturbative approach 
clearly fails if A is large. 

There is, however, a hint already at this stage that 
there are some peculiarities in the perturbation theory 
at strong coupling. Namely, we expect (and we show be- 
low) that at strong coupling, spin fluctuations are com- 
pletely overdamped, i.e., the frequency dependence of the 
spin susceptibility is dominated by iuj term instead of uj 2 
term. Crudely, this effect can be modeled by reducing 
v s in ( |l3| ) and ([To]). We see, that when vp/v s becomes 
very large, the regular part of the fermionic self-energy 
vanishes leaving only the anomalous piece in (U3). Si- 
multaneously, Ag/g, i.e., vertex correction also vanishes. 
One can easily make sure that this tendency persists to 
all orders of the perturbation theory. 
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Alternatively speaking, at v s — ► 0, the vertex correc- 
tion and the regular piece of the fermionic self-energy 
both vanish leaving the anomalous piece in the fermionic 
self-energy as the only term one has to evaluate. 



4- spin polarization operator 

To verify that the spin-fermion interaction reduces v s 
and makes spin fluctuations softer, we now compute the 
bosonic self-energy It comes in our model from a vir- 
tual processes in which a spin fluctuation decays into a 
particle-hole bubble, and hence coincides, up to an over- 
all factor, with the fully renormalized spin polarization 
operator n(q, u). We define II(q,u;) via x _1 (q, u>) = 
Xo' 1 (q,w) - n(q,w)/(xo£ 2 ), or 



x(q»w) 



x e 



l + ^ 2 (q-Q) 2 -n(q,^- 



(18) 



The contributions to the particle-hole bubble from high- 
energy fermions are already absorbed into the bare sus- 
ceptibility, and in our consideration we have to consider 
only fermions with energies smaller than A. Also, as 
static xo(q) is assumed to be peaked at q = Q, it is 
sufficient to compute the polarization operator only at 

q = Q 



k,Q,p 




k+Q, H+co, a 



FIG. 8. The lowest-order diagram for the spin polarization 
operator. 

Let's again start with the lowest order perturbation 
theory in the spin-fermion coupling. The lowest-oder di- 
agram for II(Q,w) is presented in Fig. g. In analytical 
form, we have in real frequencies 



n(Q,w) 



2§e 



cPkdQ 

3- G (k + Q, ft + u) G (k, ft) 



(2tt) 3 



(19) 



(a factor of 2 comes from a summation over spin com- 
ponents). Expanding both fermionic energies at k and 
k + Q to linear order in the deviations from hot spots, 



integrating over momentum and frequency in ( |l9| ) and 
multiplying the result by the number of hot spots N = 8, 
we obtain, 



n(Q,w) = i uj/oj s { 



where 



" sf= N 



47T v x v y £, 2 _ j5_ v x v y v F £, 1 
a = N vl 



A 



(20) 



(21) 



Instead of explicitly inserting N = 8, we will keep N as 
a variable - later we perform a formal expansion in 1/N. 

The independence of ImII(Q, w) of the fermionic cutoff 
A is indeed the consequence of the energy conservation 
requirement which confines fermions to hot spots. Ob- 
serve, however, that Eq. ( p?j| ) is the full expression for 
II(Q,w), not only its imaginary part. In other words, 
as long as one restricts with the linearized fermionic dis- 
persion near the Fermi surface, ReU(Q,u>) = 0. This 
last result implies that there is no universal correction to 
£ from low energy fermions. The corrections to £ only 
appear when we expand quasiparticle energies further in 
deviations from the Fermi surface. It is easy to check that 
these corrections scale as A/W and are therefore just mi- 
nor leftovers from the contributions from high fermionic 
energies. 

We see that the lowest-order spin polarization operator 
has a simple relaxational, linear in u> form. This implies 
that at low frequencies, the renormalized spin suscepti- 
bility more strongly depends on frequency than the bare 
susceptibility. This confirms our early conjecture that the 
effect of the bosonic self-energy can be crudely modeled 
by increasing v s . 

We next estimate what bosonic momenta and frequen- 
cies mostly contribute to the fermionic self-energy in dq) 
and to the vertex correction in (O ) . Simple power count- 
ing shows that the integrals in (|6|) and (JT3) are dominated 
by |q — CJ| ~ £ _1 and u ~ u.f£ -1 . Comparing the full 
and bare susceptibilities at these momenta and frequen- 
cies, we immediately see that for A > 1 EI(Q,w) ~ A 
dominates over bare xoi c li LU )- This implies that for 
large A, spin fluctuations which mostly contribute to the 
fermionic self-energy, are completely overdamped due to 
a strong decay into a particle-hole pair. Clearly, in this 
situation, one cannot neglect the transmutation of the 
spin dynamics in evaluating the fermionic self-energy and 
the correction to the spin-fermion vertex. 



B. renormalized perturbation theory 

Our next strategy is the following. We assume for a 
moment that Eq. d2G) for II(Q, u) with cu s f given by 
( pl| ) is valid for all couplings, and re-evaluate the lowest- 
order fermionic self-energy and vertex correction using 
the renormalized form of the spin susceptibility. We 
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then check the self-consistency of this procedure by re- 
evaluating II(Q, uj) with renormalized fermionic propa- 
gators and vertices. 



1. fermionic self-energy 

Consider first the fermionic self-energy. Substituting 
the renormalized %(q, ui) instead of xo(q I < jJ ) into Eqjq, 
subtracting, as before, the self-energy at uj = ek+Q = 0, 
and expanding to linear order in k — k^ s , we obtain 



E(k > n) = (tfi-e fc+< ,) J(k,w) 



(22) 



but now 



I(k,n)=3ge 



1 



d qduj 

T^F i + m 2 + H/^f + #* 
i i 



ifl - e k+Q +iu — v F q x iuj - v F q x 



(23) 



Consider, as before, the limit fl, ek+Q — * and com- 
pute first the regular piece /(kfc s ,0). Without ek+Q and 
ft, the last two terms in (£3|) are equal and produce a dou- 
ble pole at v F q x = iuj. Closing the integration contour 
over q x by a semi-circle in a half-plane with no double 
pole, we find that we only have to consider the pole in 
the spin susceptibility. This pole is located at a much 
larger q x oc \fu~j than the double pole. In this situation, 
one can neglect uj in the Green's functions in comparison 
with q x ex y/uj. The remaining computation of I(kh s , 0) is 
straightforward. The details are presented in Appendix 
A, and for X^$> 1, the result is, with logarithmic accuracy 



IregyK-hsi UJ 



l2v x V y 

„ r 9 log A 

ttNv 2 f 



(24) 



We see that the regular piece in /(k, ft) only logarithmi- 
cally depends on A, in distinction to a direct perturbation 
theory where it was of order A. This result is indeed con- 
sistent with our earlier observation that l reg (kf ls ,0) in 
(|l(j) vanishes in the limit v s — > 0. The still presence of 
the log A factor in ( p4[ ) is a consequence of II(Q,w) oc ui, 
and could not, indeed, be anticipated in a perturbation 
theory with a bare spin susceptibility. 

For completeness, we also present the expression for 
^reg(k/i S ,0) at arbitrary A. The calculations are pre- 
sented in Appendix A. The result is particularly simple 
for v s — > oo, i.e., for purely static bare susceptibility. We 
obtained 



A-eg (k/ ls , 0) — —A 



1+T 



(25) 



where a — (vf£, 1 /^sf) — (NX/3) (v F /v x v y ). For a <C 
1, we recover the result of a direct perturbation theory, 
-A. For a > 1, Eq. © reduces to (|J). 



^reg [K-hsi tlj 



We next compute an anomalous term in /(k, CI) which, 
we remind, is a contribution from a regularized double 
pole in /(k, Q). Integrating, as before, in Eq. fl23[ ) over 
a tiny region — Q < uj < where the poles in the two 
fermionic propagators are located in different half-planes, 
we obtain 



i an (k,n) = \- 



ifi 



fl- 



e k+Q 



(26) 



The fact that I a n(k, Q) remains the same as in a direct 
perturbation theory can be easily explained. Indeed, the 
integration over uj in I an (k., O) is confined to vanishingly 
small frequencies for which the dynamical term in the 
spin susceptibility is small no mater whether the dynam- 
ics is ballistic or relaxational. 

Substituting the full form for I(k, fl) into (g2|) we ob- 
tain E(k, ft) = T, reg (k, ft) + E on (k, ft), where for A > 1 



T, an (k,ft) = iXfl 
E reg (k,fi) = (ifi-e k+Q ) % y log A 

7TiV Dp 



(27) 



We see therefore that the fermionic self-energy con- 
sists of two parts - a regular part which accounts for the 
renormalization of the quasiparticle residue, and a sin- 
gular part which depends only on uj. The regular part 
is strongly reduced compared to the direct perturbation 
theory, and scales as log A rather than A. This reduction 
is due to the fact that at frequencies which determine 
the fermionic self-energy, the renormalized spin excita- 
tions obey relaxational dynamics. On the other hand, 
the singular part of the self-energy comes from very low 
internal frequencies, comparable to the external one, and 
is totally insensitive to the change of the spin dynamics. 
At strong coupling, the singular, nonperturbative piece 
in E is much larger than the regular piece and obviously 
should play a central role in all our analysis. 

In more general terms, the appearance of I an (k., ft) is 
the reflection of the singularity in the particle-hole bubble 
at small momentum and frequency transfer. This follows 
from the observation that for small VFq x and uj' the spin 
susceptibility, integrated over q y , can be approximated 
by a constant, and the anomalous contribution to I(k, uj) 
can be reexpressed as 



Uk,fi) = An(k,n) ; 



(28) 



where 



n ( k . o , = - ; / ^^ G (q x ,uj) G (k + q x , uj + ft) 



{2-Kf 



(29) 



is the particle-hole polarization bubble with a small mo- 
mentum/frequency transfer (it differs by an overall con- 
stant from II(k, uj) introduced in (|l8|)). At vanishing k 
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and f2, this polarization bubble is formally ultraviolet 
divergent, and its value hence depends on how the reg- 
ularization is performed |30| , i.e., on the order of the 
integration over frequency and over momentum. Doing 
frequency integration first, one obtains that n(k, f2) van- 
ishes, while doing momentum integration first, one ob- 
tains that 



n(k,o) 



n 



Q — Vi?k 



(30) 



In our case, the regularization of the self-energy is im- 
posed by the the fact that /(k, uS) is actually a convolu- 
tion of the polarization bubble and the spin susceptibility. 
In the above calculations we neglected the dependence of 
xG w ) on Qx and u>. In reality, however, x(q, w) indeed 
vanishes when either q x or u> diverge. The general rule of 
the ultraviolet regularization procedure is that it has to 
be performed such that to avoid generating extra poles 
at energies smaller than the energies of the poles in the 
particle- hole bubble, Eq. @. The pole in @ is lo- 
cated at oj oc |q|. The pole in the spin susceptibility 
x(q, w) is located at w oc q 2 , i.e., at a much smaller fre- 
uency, and at a much larger q x than those for the pole in 
29|) . Clearly, then, the correct way to regularized (£9|) is 
to first integrate over momentum in (|29|), and only then 
integrate over frequency. This sequence of integrations 
yields the anomalous piece in TI(k, SI), Eq. (|30|), which 
in turn gives rise to the anomalous piece in the fermionic 
self-energy (|28|). 

Note by passing that in a direct perturbation theory 
(which operates with the bare susceptibility), a typical 
frequency scales with a typical momentum. In this sit- 
uation, the result of regularization strongly depends on 
the ratio of vp and v s (see Eq. (|l3|)). In particular, the 
absence of the frequency dependence of the self-energy at 
v s = oo is caused the fact that in this situation \o (Qi w ) 
is purely static, and the regularization obviously has to 
be performed by doing frequency integration first. 



2. vertex correction 

We now re-evaluate the vertex correction with the re- 
laxational form of the spin susceptibility. The analyt- 
ical expression for Ag/g is the same as in (14), but 



with x(q, uj) instead of xo(q)- Linearizing, as before, the 
fermionic dispersion near the Fermi surface, and perform- 
ing momentum and frequency integration as described in 
Appendix A, we obtain for A ^> 1, with logarithmical 
accuracy, 



where 



Ag_ 
9 

Q(v) 



Q(y) 

N 



log A 



tan 



(31) 



(32) 



is a smooth function of the ratio of velocities interpolat- 
ing between Q = 1 for v x = v y , and Q — 2 for v y — > 0. 
The last limit corresponds to almost nested Fermi surface 
at hot spots. A similar logarithmical form of the vertex 
correction has been obtained by Altshuler,Ioffe and Mil- 
lis @. 

We see that the change of the spin dynamics has a 
strong effect on the strength of the vertex correction: in 
a direct perturbation theory, it was of order of A, now it 
scales only as log A. This result is fully consistent with 
our earlier observation from a direct perturbation theory 
(Eq. (|T7|)) that Ag/g vanishes when v s — v 0. As with 
the k— dependent piece in the self-energy, the subleading 
logarithmical dependence of Ag/g could not, indeed, be 
anticipated in a direct perturbation theory with the bare 
spin susceptibility. 



3. higher-order diagrams 

We next need to understand what happens when we go 
to higher-orders in the self-consistent perturbation the- 
ory for fermionic self-energy and the spin-fermion vertex. 
The issue is indeed whether higher order diagrams scale 
as higher powers of A or log A. We argue that higher- 
order terms contribute higher powers of log A but not 
higher powers of A. Our reasoning is based on the fact 
that in the second-order theory, the vertex correction and 
the regular part of the fermionic self-energy are 0(log A), 
and 0(A) term in the second-order fermionic self-energy 
emerges only as a result of a proper regularization of the 
double pole in the integrand for dY,/dfl. Higher-order 
self-energy and vertex correction diagrams contain more 
fermionic propagators with different momenta, and one 
can easily make sure that there is no phase space for a 
double pole. Accordingly, higher-order vertex correction 
diagrams just contribute higher powers of log A, while 
higher-order diagrams for the fermionic self-energy cither 
simply contribute higher powers of log A, or account for 
log A corrections to E on (k, ft). We verified this argument 
by explicitly computing next order self-energy and vertex 
corrections. 



4- spin polarization operator 

A closely related issue is how the fermionic self-energy 
and vertex renormalization affects the form of the spin 
polarization operator. We recall that above we evalu- 
ated the fermionic self-energy and Ag/g assuming that 
that the spin polarization operator has the same form 
as in the direct perturbation theory. We now need to 
verify to which extent this assumption is correct. We 
show that the renormalization of II(Q,w) from its free- 
fermion form, Eq. (EOf) is only due to vertex corrections 
and to the regular part of £(k, fi), while a much larger 
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S an = Xfl does not affect the spin polarization opera- 
tor. To demonstrate this we re-evaluate II(Q, ui), Eq. 
(H), with G(k,fi) = fi + S on (0) - e fc . The computa- 
tion is rather straightforward as the inclusion of S an into 
fermionic propagator modifies the quasiparticle residue 
Z and the quasiparticle mass (defined as vp = pp/m) 
as Z = 1/(1 + A), m* = m(l + A). Substituting the 
renormalizcd form of G(k, f2) into (H3) and performing 
elementary manipulations, we find that Z and m* /m ap- 
pear in II(Q,w) in the combination Z 2 (m* /m) 2 which is 
independent of A. In other words, as long as we restrict 
with only S on (k, f2), II(Q,cl>) remains the same as for 
free fermions. 

The inclusion of S re / (k, £1) into the fermionic propa- 
gator and the corrections to the spin-fermion vertex does 
affect the form of II(Q, w), but the corrections obviously 
hold in powers of log A. 



C. summary of Sec II] 



We now summarize what we obtained in this section. 

i We found that the direct perturbation theory for 
fermionic self-energy and spin-fermion vertex holds 
in powers of A oc ()I(vf£,~ 1 ) and does not converge 
at strong coupling. 

ii We found that for typical bosonic energies which 
mostly contribute to fermionic S(k, tt) and to ver- 
tex renormalization, the bosonic self-energy also 
becomes relevant for A > 1 and accounts for the 
change in the bosonic dynamics from a ballistic one 
to a relaxational one. 

iii We performed the self-consistent perturbative ex- 
pansion in which we used the relaxational form of 
the spin susceptibility as an input, and found that 
it holds in powers of log A, plus there exists the 
anomalous, nonperturbative piece in the fermionic 
self-energy which still scales as A. At strong cou- 
pling, this piece much larger than the regular self- 
energy and obviously should play a central role. 

iv To logarithmical accuracy, the anomalous piece in 
the fermionic self-energy is given by the second- 
order diagram, higher-order terms just add extra 
powers of log A. 

v We found that as long as we restrict with only O(A) 
piece in the self-energy and neglect vertex correc- 
tions and the renormalization of the spin-fermion 
vertex, the spin-polarization operator remain the 
same as for free fermions, i.e., the self-consistent 
perturbation theory is justified. 



IV. N -> oo LIMIT, T = 

These results of the previous section imply that ne- 
glecting logarithms, one can construct a fully self- 
consistent new "zero-order" theory. In this theory, spin 
excitations are purely relaxational and vertex corrections 
are absent. At the same time, the (anomalous) fermionic 
self-energy is strong and progressively, as A increases, de- 
stroys fermionic coherence. 

This new "zero-order" theory indeed makes sense only 
if one can specify the limit when it becomes "exact" 
and also construct a controllable perturbative expansion 
around it We notice in this regard that all log A terms 
contain a factor 1/N where, we remind, N(= 8) is the 
number of hot spots in the Brillouin zone. The presence 
of 1/N is a direct consequence of the fact that the log A 
terms appear due to a change from a ballistic to a re- 
laxational spin dynamics. The prefactor for logarithm 
should then be inversely proportional to the spin damp- 
ing rate as when this rate is zero, the self-energy and 
vertex corrections are much larger and scale as 0(A). A 
spin fluctuation with momentum Q has N/2 channels to 
decay into particle-hole excitation near hot spots, hence 
the damping rate scales as N, and the prefactor for the 
logarithm scales as 1/N. Below we formally treat N as an 
arbitrary number. Then one can define the N ~ » oo limit, 
when logarithmical corrections can be totally neglected. 
In this limit, our new "zero-order" theory becomes exact. 
At finite but small 1/N, the perturbation theory around 
our new vacuum is the expansion in 1/iVlogA . 

Instead of introducing an artificially large number of 
hot spots, one can also extend the spin-fermion model to 
a large number of electron flavors M (= 1 in the physi- 
cal case) and expand in 1/(8M). This later expansion is 
more appealing from physics perspective as the extension 
to large N requires a rearrangement of the Fermi surface 
to a large number of "hot" segments which is difficult 
to visualize. In both cases, however, the main idea is 
to enhance the effect of the bosonic damping and there- 
fore reduce the strength of regular self-energy and ver- 
tex corrections compared to the anomalous term in the 
fermionic self-energy. The computations are identical in 
both cases, and below we will just label the expansion as 
"large N" theory. We discuss some formal aspects of the 
1/N expansion in Appendix Kill . 

The limit N = oo bears some similarity to the Migdal- 
Eliashberg limit for electron-phonon problem [|Z8| — |80f1 al- 
though in our case bosonic modes are not independent de- 
grees of freedom. Nevertheless, the smallness in 1/N has 
the same consequence as the smallness of m/M where m 
is the electronic mass, and M is the ionic mass: one can 
neglect vertex corrections and also the momentum de- 
pendent piece in the fermionic self-energy. Below we will 
occasionally refer to the N — oo limit as the Eliashberg- 
type theory. 
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We now proceed with the more detailed discussion of 
the N = oo limit. 



A. derivation of the Eliashberg-type equations 

So far our analysis of the fermionic self-energy was re- 
stricted to the limit of vanishing fermionic frequency fl. 
In this limit, we have at N — ► oo, E(k, Q.) — Xfi. We have 
also checked that for this E(k, O), the spin polarization 
operator is the same as for free fermions. 

We now derive the full form of the fermionic self-energy 
at N — oo, valid for arbitrary frequencies. The input for 
the derivation is the knowledge that vertex corrections 
and the corrections to tk both scale as 1/N and vanish 
at N = oo. Without them, we have, quite generally, 

G(k,«) = w — ^ , (33) 

where E(O) = O + E(faifi), and fcn is the momentum 
component along the Fermi surface. The dependence 
of the self-energy on fcy is a subleading effect which we 
will be discussing below. This effect is relevant at the 
lowest frequencies but becomes progressively weaker as 
$7 increases. We explicitly verified a'posteriori that at 
N — > oo, typical frequencies are such that the dependence 
of G(k, $7) of intermediate fermions on fey can be totally 
neglected. We therefore just ignore it in the derivation 
below. 



a) 



a 



£gl °" O 



FIG. 9. The diagrams for the fermionic self-energy and 
spin polarization operator at N = oo. 

The full diagram for the self-energy at N — oo is pre- 
sented in Fig p|a. Another diagram in Fig |p3 is for the 
full spin polarization operator. We cannot a'priori use 
the free fermion form for II(Q, u>) as so far we have only 
demonstrated that it is preserved if we use the small fi 
form of the fermionic self-energy What we have to do 
is to derive and solve a set of two coupled equations for 
E(fc||fi) andn(Q,u/). 

The evaluation of the diagrams is straightforward. 
Let's first consider fermions near hot spots when 
E(fc||fi) « E(fl). Substituting Eqs. (|§) and (|l|) into 
the diagram for the self-energy and integrating over mo- 
menta along the Fermi surface at k; is + Q, we obtain 

s (^) = ■JTZ^ / d<jjdq x — 



(34) 



8 ^ 2 j - [g + £ 2 (l-II(Q, W ))]Va 

E(u>) - E(w + fi) 
(v F q x - Y,(uj))(v F q x - E(w + 0)) 



where as before q = q — Q. 

We already know that the "regular" part of the self- 
energy vanishes at N — > oo, so let's focus on the "anoma- 
lous" part which should reduce to E an (f2) in the O — ^ 
limit. This piece comes from the range of frequencies 
where the two fermionic poles are in different half-planes 
of q x . Since signE(w) = signer, this u> range is sandwiched 
between — Q < lo < (for definiteness we set tt > 0). The 
next step is to integrate over q x . At small fi, we neglected 
q x as its inclusion would yield higher powers of Q. For 
arbitrary f2, we cannot use this argument. Still, it turns 
out that even for high frequencies, keeping q x in the spin 
susceptibility give rise to corrections which are at least 
small by 1/N 2 . We explicitly show this below in a sep- 
arate subsection. At this moment, we just neglect q x in 
the spin propagator and explicitly integrate over q x in 
(|34|). Evaluating the integral, we find that the fermionic 
self-energy is canceled out, and E(f2) takes the form 



E(Q) = 2A / 
Jo 



du> 



Vi-n(Q,w) 



(35) 



This cancellation of the self-energy for intermediate 
fermions implies that self-consistent, FLEX-type calcu- 
lations in the normal state are not necessary as the dom- 
inant piece in E(fi) is captured already by using free- 
fermion form of the propagator of intermediate fermions. 
We next evaluate the spin polarization operator for 
arbitrary E(fi). Substituting the spin propagator ( |33| ) 
into (fL9[), replacing, as before the momentum integration 
by the integration over dekdek+Q, and integrating over 
energies, we obtain 



n(QV> 



Wsf 



du 1- 



E(w)E(w + w') 



y^HVE^w) 



(36) 



Eqs. ( |35j ) and ( |36[ ) form the set of two coupled inte- 
gral equations for the fermionic self-energy and the spin 
polarization operator. 



B. solution of the Eliashberg-type equations 



A more careful look at Eqs. ( |35| ) and (|36|) shows that 
the solution of the set is straightforward. Indeed, we 
clearly see that the magnitude and the functional depen- 
dence of E(w) is totally irrelevant for II(Q, u>), all that 
matters is the fact that sign S(w) = signw. The fre- 
quency integration in (Bq) is then straightforward, and 
performing it we find that 



n(Q,«) 



w s f 



(37) 



i.e., for arbitrary E(O) it preserves exactly the same form 
as for free fermions. Substituting this result into ( pq ) we 
immediately find that 
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£(51) = 



2A51 



1 



1 — i— - 



(38) 



This expression can also be re- written in the scaling form 



S(f2) = XuJsfg 



UJ 



U> s f 



g(x) 



2x 



i + yr^TN 



(39) 



At small 51 we indeed recover the previous result £(51) = 

A5i. 

The general explanation why the fc— independent 
fermionic self-energy £(51) does not affect the polariza- 
tion bubble at a finite momentum was given by Kadanoff 
| p3[ . He pointed out that the expansion of the fermionic 
energy to first order in momentum deviation from kp is 
equivalent to imposing an approximate Migdal sum rule 
on the spectral function A(k, uj) — (1/n) ImGo(k,uj) 



J<kkA(k,u) = l 



(40) 



Expressing Go(k, w) in terms of -A(k, uj) via Kramers- 
Kronig relation and making use of (Eoh , one finds 



n(Q,o/) 



Wsf J_ 



aw — ; — U) 
duj 



-0(uj ,: 



(41) 



where f(u>) is the Fermi function. Eq. (EOJ) then follows 
from the fact that /(uj) is 1 at ui = — oo and at u> = +oo. 

Eqs. ( p7| ) and (pq) are the central results of A = oo 
analysis. We see that the scale for nonlinear effects in 
frequency in £(51) is set by the typical spin relaxation 
frequency oj a f = A-Kv x v y ^ 2 / (N g) (the 1/N factor in u> a f 
can be eliminated by res caling v — > vN and g — » gN 
as we discuss in Appendix XIII ). This u s f is obviously a 
measure of the deviation from the QCP. At scales smaller 
than u) s f, spin susceptibility is nearly static, and £(51) 
is expandable in 51. In this regime, we should gener- 
ally expect a Fermi liquid behavior with S"(J1) oc u> 2 . 
However, above U) s f, the system should cross over into 
a quantum-critical regime, where the system behavior is 
determined by the fixed point at £ = oo. In this regime, 
spin excitations become massless diffusons, and fermionic 
self-energy should modify such that to eliminate the de- 
pendence on £. 

We indeed find this behavior in our E(f2). We see from 
( p8J ) that for |fi| < u) s f, the self-energy behaves as 



E(fi) =AO 1 + i 



101 



4u> 



sf 



(42) 



For 1 51 1 > ut s f, however, it crosses over to a different 
behavior 



£(ft) = (i\n\ui) 1/2 signil 



(43) 



We see therefore that above w s f, both S'(S1) and E"(51) 
scale as ^/|0|. The normalization energy 



= 4A 2 



J sf 



9g 2v x v y 
2nN vl 



(44) 



does not depend on £ as it should be in the quantum- 
critical regime. This Q sets the upper cutoff for the 
quantum-critical behavior (its dependence on N can 
again be eliminated by rescaling g — > Ng. 

The vQ form of the self-energy at the antiferromag- 
netic transition was first obtained by Millis [|l[ . The yfui 
form also emerges in the quantum-critical models with 
disorder gj. 

Substituting the quantum-critical form of the self- 
energy into the fermionic propagator, we obtain that at 
k sa k^ s and 101 < u> 



G(k,n) 



1 y^|fT|cJsign51 - 



Efc 



mi 



/« 



(45) 



We see that in the quantum-critical regime, the system 
behavior is qualitatively different from that in a Fermi 
liquid: there is no pole in the fermionic propagator at 
real frequencies. Instead, G(k, O) has a pole along an 



imaginary frequency axis, at |0| 



<S1 



el/u>. This 



non-Fermi-liquid pole gives rise to a broad maximum at 
51 = ef . We also verified that in the crossover region 
51 ~ ^sf, the pole in the fermionic propagator gradu- 
ally moves, with increasing 51, from real to imaginary 
frequency axis 

At frequencies larger than u>, the bare to term in the 
fermionic propagator wins over the self-energy, and the 
spectral function recovers the peak at 51 = e&. Still, this 
behavior is not a Fermi-liquid one as the width of the 
peak scales as y/fi. 




FIG. 10. The behavior of ImG(k, 51) at w af = @. The 
inset shows the "quasiparticle" dispersion inferred from the 
position of the maximum in the spectral function (indicated 
by arrows). 

The behavior of ImG(k, 51) at w s f — is shown in Fig. 
|l0| (see Eq. (EH)). The inset shows the "quasiparticle" 
dispersion inferred from the position of the maximum in 
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the spectral function. We see that the destruction of the 
Fermi liquid behavior accounts for the effective flattening 
of the "quasiparticle" dispersion at frequencies below w. 
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FIG. 11. Real and imaginary parts of the fermionic 
self-energy, Eq.(B9(), vs frequency. 

In Fig. nil we present both real and imaginary parts 
of the fermionic self-energy vs frequency at a finite u> s f . 
We clearly see a Fermi-liquid behavior at the smallest 
frequencies, and the y/U behavior well above uj s f. By 
numerical reasons, the crossover region in between these 
two limiting regimes is rather wide - it stretches between 
0.5w s / and 6 — 8w s /. In the crossover regime, E"(k, u>) 
is, to a surprisingly good accuracy, a linear function of 
frequency. This linearity, however, does not follow from 
any theory considerations, and is just a hidden numerical 
property of the self-energy, Eq. (p8|). This linearity, how- 
ever, is relevant for the explanation of the photocmission 
and optics data in cuprates. 
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In Fig [T^ we present the forms of the spectral func- 
tion 7mG(k, fi) for different values of efc and finite 0J s f- 
Comparing this result with Fig [10| we see that the quasi- 
particle peak sharpens up at frequencies comparable to 
uj s f, as it indeed should as the system crosses over to 
the Fermi-liquid behavior. Notice, however, that even 
at lu ~ Lu s f , the width of the peak is comparable to its 
amplitude. 

We next discuss to which extent the full anomalous 
self-energy depends on the momentum along the Fermi 
surface. Repeating the same calculations which lead to 
( p8| ) but for a finite k = k k^ s , we find that S(k, Ci) = 
S(fc|i, fi) can be rewritten as 



E(*,n) = A(*)- 



_><> 



1 



1 



|0| 



(46) 



f (k) 



QJ(£> M 



i.e., it still has the same form as at a hot spot, but with 
k dependent 

A(fc)=A/(l + (fc0 2 ) 1 / 2 , u> 8f (k)=u sf (l + (k0 2 ) (47) 

We see that away from the hot spots, the effective cou- 
pling gets smaller, and the crossover frequency uj s f(k) 
increases. Still, however, at frequencies which exceed 
u> s f(k), the system displays the same non- Fermi liquid 
behavior as at hot spots, i.e., E"(fc,f2) first scales lin- 
early with fl, and crosses over to a V" form at larger 
frequencies. 

In the Fermi liquid regime, the slope of E" depends on 
k as (1 + (&£) 2 ) ^ 2 - I n the non-Fermi liquid regime, this 
dependence gets weaker and completely vanishes in the 
yVL regime where the overall factor in E(fe, 51) becomes 
X(k)\/uJ s f(k) — Xy/uJsf- This can be seen already from 
(|38|). This implies the momentum range where E(fc,0) 
weakly depends momentum (i.e., the "size" of a hot spot) 
depends on frequency. For 57 < lu s j, this range is obvi- 
ously constrained by |fc|£ < 1, but for larger frequencies it 
increases and eventually extends to \k\£ < (lill/ws/) 1 / 2 . 
In particular, deep in the quantum-critical regime, i.e., 
for frequencies comparable u = 4A 2 w s /, the weak mo- 
mentum dependence of the fermionic self-energy extends 
to |k — kfj S | < g/vp, i.e, the "size" of a hot spot remains 
finite even at £ — > oo. This result implies that in the 
quantum-critical regime the momentum variation of the 
self-energy does not play any significant role and can be 
safely neglected. 



C. the accuracy of the Eliashberg-type theory 



FIG. 12. The behavior of ImG(k, Q) at a finite w s /. For 
definiteness we used A = 1.7 in which case Gj = 11.56w s /. 
The hig h-frequency behavior of ImG(k, Q.) is the same as in 
Fig.©. 



We now have to go back and verify the accuracy of 
the Eliashberg-type Eqs. (B5|) and (p6h. We begin with 
the fermionic self-energy. Recall that in deriving ( |3q ) 
we neglected the dependence on q x in the spin propa- 
gator. Near the singular pole, q 2 ~ (E(w) / 'v f) 2 ■ To 
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estimate the relative strength of this term, let's evalu- 
ate it with S(w) from ( J3q ) obtained by neglecting this 
contribution. At u> < Lo s f, S(w) w Xuj, hence typ- 
ical q 2 ~ X 2 uj 2 /v 2 f = (Lo^)(Lu/(AN 2 uj s f)) are negligi- 
ble compared to ory. Clearly, q 2 term can be omit- 
ted. In the quantum-critical regime, lu ^$> t-u s f, we have 



<n 



jioj/vp = {ijj"f){a>v x v y /Nv 2 F ) 2 . We see that the term 



p — yjj j jyuu x u y /i\ up 

which wc omitted in (35) and, hence 



has the same 

frequency dependence as the damping term in the spin 
susceptibility, but has an extra factor (1/N) 2 and there- 
fore can be neglected in the N — > oo limit. 

To understand how 1/N expansion works in terms of 
numbers, we evaluated the overall factor for £(w) in 
the quantum-critical regime by solving Eq. ([34]) self- 
consistently for a physical N = 8. Substituting E(fi) = 
A(i|r2|w) 1 / 2 signcj as an input into (p4[), we indeed ob- 
tained the same form in the output. Solving for A, we 
obtained A — 0.94 which is very close to A — 1 - the 
result without 1/7V corrections. 

We also compared the N — oo result with the full 
second-order expression for £(f2) at arbitrary fi. The 
calculations are presented in the Appendix A. The result 
is particularly simple for v s = oo in which case we have 
from fll6g|) 



E(il) ~ 2Xu s fa In 



where 



iy/Kn - 1 + y/Kn + 1 



iy/Kn-l+An + y/K n + 1-As 



(48) 



Kl 



■■0,f 



2i to 

An = ^ (49) 



and, we remind, a = (vp£~ 1 )/u> s f — (XN/3)(vp/v x v y ). 
Our N = oo result, Eq. (|3^), is the limiting from of 
(Eq) at a — > oo. The expansion to a linear order in 1/a 
yields a correction to (f58|) which at low frequencies is 
logarithmical in A and coincides with Eqs. (E2) and (24). 
In Fig. O we plot the full second-order result for the 
self-energy at N = 8 and our N — oo result which is its 
anomalous part. For defmiteness, we set v x — v y and 
A = 2 in which case a — 32/3. We see that the two 
curves are very close to each other up at least to ft = Co, 
which, we remind, is the upper cutoff frequency for the 
quantum-critical behavior. We consider this agreement 
as another evidence that the restriction with only the 
anomalous contribution to the self-energy works well for 
the physical N = 8. 




FIG. 13. Real and imaginary parts of the full fermionic 
self-energy at finite N as & function of frequency. Dashed 
lines rep resent the self-energy without the regular part (see 
Fig. Q. 

A minor comment is in order here. In the quantum- 
critical regime, the second-order expression for T,(lo) con- 
tains a leading fi 1 / 2 term and a subleading ft log Q term. 
On the other hand, performing computations with the 
full fermionic propagators, we found a subleading log- 
arithmical term in X(k/ ls ,f2) but no filogf2 term. A 
careful look into this discrepancy shows that the reason 
why the subleading (regular) piece in the second-order 
self-energy has il log ft form is because for free fermions, 
dGo 1 (\c,w)/duj = -dGo\\s.,w)/de k = 1. In the self- 
consistent calculations, we still have dG~ 1 (k,u>)/dek = 
— 1, but now 9G~ 1 (k, oj)/dui — dJ^(cu)/duj. In the 
quantum-critical regime, dY,{ui)/dui oc ui^ 1 ' 2 . The ap- 
pearance of the extra to 1 ' 2 in the denominator, trans- 
forms rHogfi term into VO/iV 2 , exactly as we found. 

Finally, we also verify that the renormalization of e k 
and of the spin-fermion vertex scale as log X/N even when 
we use quantum-critical forms of the self-energy for in- 
termediate fermions. Lets, hrst analyze £(k, 0). Since 
the fermionic self-energy in (133) does not affect e k , we 
can straightforwardly expand S(k, 0) in e k +Q- Perform- 
ing then a simple power counting, we find that in the 
quantum-critical regime, the inclusion of the self-energy 
of intermediate fermions only gives rise to 0(1/N 2 ) cor- 
rections to the second-order result. 

The same result holds for vertex renormalization. Sub- 
stituting (|3^) into ([hi]) we obtain 



>f) 



.9 (2^) 3 



1 



d 2 qduj — -„ 

4 z - i\(jj\'y + q z 

1 1 



£(w) - {v x q x + v y q y ) £(w) - {-v x q x + v y q y ) 



(50) 



We immediately see that if we neglect the momentum 
dependence of the spin susceptibility, then the fermionic 
self-energy is totally eliminated after the integration over 
d 2 q, and the vertex correction remains exactly the same 
as if we use free fermions (see Eq. ([H])). If we keep the 
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momentum dependence in x(q, u>), the functional form 
of Ag/g does not change, but the coefficient changes by 
1 + 0{1/N 2 ). 



D. vertices at small momentum transfer 

Before we proceed with the 1/N expansion, let's con- 
sider in more length the issue of vertex renormalization 
for different bosonic momenta q. Earlier we found that 
the particle-hole vertex with momentum transfer Q is 
only weakly (as log A/TV) renormalized by spin-fermion 
interaction. Mathematically, the smallness of the vertex 
correction is due to the fact that the two internal fermions 
in the vertex correction diagram had different directions 
of Fermi velocities, and hence one can independently in- 
tegrate over efc and €k+Q ■ The vertex correction then can 
roughly be written as 



-~4 I du,N 2 (w) X (Q,w) 
9 v F J 



(51) 



where N(u>) = f de^Gfe, w) = — i7rsign(w). The fre- 
quency dependence then comes only from x(Q) w ) oc - 
Vp/(Ng 2 uj), and frequency integration yields Sg/g ~ 
(L/N) where L is the logarithm of the lower cutoff. 

This reasoning is valid for all external q, not neces- 
sary Q = [it, 71") (for q ^ Q, the logarithm is cut also 
by |q — QD- However, it breaks up when one consid- 
ers a particle- hole vertex with zero momentum transfer. 
For such vertex, the Ward identity implies pd that at 
vanishing bosonic frequency Sg/g = 1 + <9£(k, oS)/du) 
for a scalar vertex (the one which does not depend on 
fcrmionic momentum and conserves spin projection) and 
Sg/g — 1 + 9£(k, cj)/dek for a current vertex (which has 
an extra dek/dk). In our situation, 9S(k, uj)/dtk is small 
in 1/N, but c?£(k, u>)/doj is not. Moreover, as £(u>) 3> u> 
at strong coupling and u) < u), the renormalized scalar 
vertex well exceeds the bare one. 

That the renormalization of the scalar vertex with zero 
momentum transfer does not contain 1/N can indeed be 
obtained in perturbative calculations. First, we derive 
the Ward identity for the scalar vertex. It follows from 
the fact that the integral equation for the full vertex 



g tot (k,LU,n)=g 



3ig 



d qduj' 



g tot (k + q,a/,») X (q,^-iy) 

(Efc+g(w' + 0) - efc +g )(Sfc+g(w') - € k +q) 



(52) 



where as before T, u — u + E w , is equivalent to the integral 
equation for Sfc(oj + Q) — Efc(w): 



Sfe(w + f2) — Sfe(w) = Q 



Zig_ 

8tt 2 



d qdui' 



(£fc +g (a/ + rt)- T, k+q (uj')) x(oj - a/, q) 
0t, k+q (u>' +fi) - e fe + 9 )(S fc + g (w') - e k+q ) 



(53) 



Comparing ( |52| ) and (^) we immediately find that 



g tot (k,u;,n)=g 1 



Sfc(u; + fl) - Sfc(^) 

n 



(54) 



n 



the r.h.s. of (ph yields 1 



In the limit 

aE(k,w)/aw. 

The very fact that g tot (k,uj,il) is related to 
<9£(k, uj)/du! implies that the large renormalization of the 
scalar vertex with zero momentum transfer is due to the 
same anomaly which accounts for large S(u>). Indeed, 
consider the lowest-order diagram for the vertex renor- 
malization. When external bosonic frequency is strictly 
zero, the two fermionic propagators have poles in the 
same half-plane of frequency, and the integration over e k 
does not vanish only due to the fact that the spin suscep- 
tibility also has poles at |q — Q| oc (Nuj) 1 ' 2 . Obviously, 
the resulting vertex correction is small in 1/N. However, 
at finite O, the poles in the fermionic propagators are 
split, and there is a nonvanishing contribution from the 
frequency range where the split poles are in different fre- 
quency half-planes. Evaluating this contribution in the 
same way as for the fermionic self-energy, we indeed find 
that 



^ = 1 + A 



dw' 



1 



£(Q - (J) + E(w') (1 + |w + <y|/u; 8 /) 

(55) 



It is clear from ( |55| ) that the large renormalization of g 
is due to the same anomaly which gave rise to a large 
S(w). Collecting the anomalous contributions from the 
higher-order diagrams and summing them up, one should 
indeed recover Eq. (|4J). 

Note that this procedure is not as straightforward 
as one might expect. Indeed, consider the limit when 
Q, lj <C u s f- I n this limit, Eq. ( p5|) yields 



Ag 
9 



= 1 



A 



1 + A 



(56) 



As for vanishing fi, the full g tot = g(l + A), one may 
suggest that the perturbation series are geometrical, i.e. 
1 + A/(1 + A)- > 1/(1-A/(1 + A))= 1 + A. However, the 
actual situation is not like this: the explicit calculation 
of the two-loop vertex correction yields 



Ag_ = 2 

9 \2loop 7r 



A 



1 + A 



(57) 



i.e., the next order diagram contains an extra n/2 
compared to what is would be if the series were ge- 
ometrical. This peculiarity of the perturbation series 
emerges because Efc(w) does possess some non-singular 
fc-dependence along the Fermi surface. This momentum 
dependence is indeed fully accounted for in the Ward 



19 



identity but complicates perturbation series in some non- 
singular way. This is just another example that momen- 
tum dependence of the fermionic self-energy complicates 
calculations but does not affect the physics. 



E. summary of Sec IV 



Let' now summarize what we have at the moment. 

i We find that at strong coupling, one can construct a 
new, fully self-consistent "zero-order" theory which 
explicitly includes the most divergent 0(A) piece in 
the fermionic self-energy, but neglects vertex cor- 
rections and the renormalization of the Fermi ve- 
locity which both scale logarithmically with A. The 
divergent piece in the self-energy comes from the 
second-order diagram. Higher order self-energy and 
vertex correction diagrams only account for totally 
irrelevant, small numerical corrections to this term. 

ii We demonstrated that this "zero-order" theory be- 
comes exact in the formal limit when N — * oo. In 
this theory, fermions are strongly affected by in- 
teraction with spin fluctuations and display non- 
Fermi liquid, quantum-critical behavior at frequen- 
cies larger than a typical spin relaxation frequency 
to s f. In the quantum-critical regime, the fermionic 
self-energy interpolates between linear in u> and -Juj 
forms. At the same time, there is no feedback from 
fermionic incoherence on spin excitations which re- 
main diffusive with the same diffusion constant as 
if they were made from free fermions. 

Our next step is to develop a perturbation theory in 
1/N. Obviously, our interest is to understand the role of 
logarithmical vertex and self-energy corrections. 



V. THE PERTURBATION THEORY IN 1/N 

A. one loop RG analysis 

We recall that the logA/JV) terms give rise to the 
two new features in the theory: vertex corrections which 
renormalize both fermionic and bosonic self-energies, and 
static fermionic self-energy S(k, 0). To the lowest order 
in 1/N, we have from ( |3l| ) and (135) 



9 


Q(v) . . 

= N log A, 


Ae fc = 


12 V x Vy 

- e k+Q »r 2 

■nN vi. 



log A 



(58) 
(59) 



We see from ( |58| , |59| ) that the 1/N corrections to the 
vertex and to the velocity of the excitations are almost 
decoupled from each other: the velocity renormalization 
does not depend on the coupling strength at all, while 
the renormalization of the vertex depends on the ratio 
of velocities only through a non-singular Q(v). The ab- 
sence of the coupling constant in the r.h.s of the Eqs. 
(|58|,|9|) is, we recall, a direct consequence of the fact 
that the dynamical part of the spin propagator is ob- 
tained self-consistently within the model. Due to self- 
consistency, the overall factors in Aefc and Ag/g are 
g(w s ft; 2 ). Since the fermionic damping is produced by 
the same spin-fermion interaction as the fermionic self- 
energy, uj s f scales as 1/g, and the coupling constant dis- 
appears from the Eqs. ( |58|]5g| ). 

The fact that the lowest-order corrections diverge at 
A = oo obviously implies that the 1/N expansion breaks 
down near the QCP, and one has to sum up the series of 
the logarithms. We will do this in a standard one-loop 
approximation by summing up the series in (l/iV)log£ 
but neglecting regular 1/N corrections to each term in 
the series (and, indeed, we neglect regular 1/N 2 correc- 
tions from FLEX- type diagrams). We verified that in 
this approximation, the cancellation of the coupling con- 
stant holds even when g is a running, scale dependent 
coupling. This in turn implies that one can separate the 
velocity renormalization from the renormalization of the 
vertex to all orders in 1/N. 

Separating the corrections to v x and v y and performing 
standard RG manipulations, we obtain a set of two RG 
equations for running v x and v^ 



K 

AL 



12 



R\2,,R 



(<)X 



irN (v*) 2 
12 



■ «) 2 
«) 2 k? 



nN (v*) 2 + «) 



R\2 



(60) 



where L = log£. The solution of these equations is 
straightforward and yields 



v* = v x Z;v* = v v Z- 1 ;Z = 



24Ln 



TtN V x 



1/2 



(61) 



where, we recall, 



and 



are the bare values of the 



velocities (the ones in the Hamiltonian) . 



where Q(v) = (4/7r)tan 1 (y x /v y ) interpolates between 
Q = 1 for v x = v y , and Q — 2 for v y — ► 0. 
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FIG. 14. The nesting of the Fermi surface as the result 
of the renormalization of the Fermi velocity. Solid line - the 
Fermi surface. Dashed line - the shadow Fermi surface (shifted 
by antiferromagnetic Q). At £ = oo, the velocities at hot 
spots become perpendicular to the magnetic Brillouin zone 
boundary. This creates a "bottle neck effect" : the Fermi sur- 
face easily disconnects at hot spots immediately below the 
magnetic transition. 

We see that v y vanishes logarithmically at £ — > oo. 
This implies that right at the QCP, the renormalizcd 
velocities at k^ s and k/ ls + Q arc antiparallel to each 
other, i.e. the Fermi surface becomes nested at hot spots 
(Fig. Ill) . This nesting in turn creates a "bottle neck ef- 
fect" immediately below the criticality as the original and 
the shadow Fermi surfaces approach hot spots with equal 
derivatives (see Fig. |l4|). Quite generally, the nesting at 
hot spots is the first step in the evolution of the Fermi 
surface towards hole pockets |84[] . If the nesting occurred 
at some finite £, (as one can judge by formally extend- 
ing the lowest order result, Eq. (|5Jj) outside the range of 
its applicability), then the system would develop strong 
SDW precursors already in a paramagnetic phase [g5| g8| . 
It turns out, however, that this process is precluded be- 
cause nesting at finite £ (i.e., v R = at finite v R ) would 
imply that Lu s f renormalizes to zero, but renormalized 
Lj s f is the overall factor in the velocity renormalization. 
The only way to avoid this negative feedback effect is 
to consider the formal limit of small N, when the spin 
damping can be neglected p(J . In this limit, the feedback 
effect on the velocity renormalization is absent, and the 
Fermi surface evolves towards hole pockets already at a 
finite £. 

Another feature of the RG equations ( pQ ) is that they 
leave the product v x v y unchanged. This is a combination 
in which velocities appear in u> s f. The fact that v x v y is 
not renormalized implies that without vertex renormal- 
ization, w s /£ preserves its form, i.e. spin susceptibility 
still has a simple diffusion pole. 

We next consider vertex renormalization. Using again 



the fact that guj s f does not depend on the running cou- 
pling constant, one can straightforwardly extend the 
second-order result for the vertex renormalization, Eqn 
(pa), to the one- loop RG equation 



dg R 
dL 



Q{v) R 

N 9 



(62) 



where g R is a running coupling constant, and Q(v) is 
the same as in fl5q) but with renormalized velocities v x 
and v R instead of v x and v y . Observe that the coupling 
constant increases as £ — ► oo. 

At the QCP, the dependence on £ under the logarithm 
transforms into the dependence on frequency and exter- 
nal q — Q: 



L -> (1/2) log 



max(M,(q-Q) 2 /7)' 



(63) 



Rl R 

J y I u x 



Using the fact that for £ — > oo, i 
Nw(Vy/v x )/(24:L) and expanding Q(v) near v R = 0, we 
find Q(v) » 2(1 - (2/tt)v r /v£) = 2 - N(v y /v x )/(6L) 
Substituting this result into (|62| ) and solving the differ- 
ential equation we obtain in the limit L ^> N 
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g exp 



(2L/N)L- V y/ &V * 



(64) 




FIG. 15. 
tion. 



The polarization bubble with the vertex correc- 



We see that the logarithmic factor L appears in the 
exponent, i.e., series of logarithmical corrections even- 
tually yield power-law dependence of the full vertex on 
frequency and momentum. This result contradicts the 
assertion that the quantum critical point is described by 
a purely bosonic "0 4 " theory at the upper critical di- 
mension d + z = 4 as in the latter case all logarithmi- 
cal corrections should sum up into the powers of loga- 
rithms. We explain the origin of the power-law behavior 
below, and in the next subsection compare our results 
with the bosonic c/> 4 theory in some detail. Here we just 
note that the very fact that L appears in the exponent 
implies that we need to know precisely in which com- 
bination 7|w| and (q — Q) 2 appear in L. To address 
this issue, we explicitly computed the polarization bub- 
ble with a vertex correction at a nonzero momentum and 
frequency, Il2(q, ui). The corresponding diagram is pre- 
sented in Fig. [TR We normalized n 2 (q, to) such that 
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X (q, «) = Xo/(r 2 + (q - Q) 2 + tM + n 2 (q, W )). The 
computation of the diagram is presented in Appendix C. 
The result is 



n 2 (q,w) 
Q(v 



N 



((q-Q) 2 + 7 M)log 



^jLO 



(q-Q) 2 + 7H 



(65) 



We see that 7w s / and (q — Q) 2 are renormalized in ex- 
actly the same way. This result implies that the spin 
diffusion coefficient does not undergo logarithmical renor- 
malization, i.e., the dynamical exponent z remains equal 
to two despite singularities. Earlier two of us missed the 
renormalization of the (q — Q) 2 term and erroneously 
concluded |61|] that z becomes larger than 2. 

That both u> and (q — Q) 2 terms in the bosonic prop- 
agator undergo logarithmical renormalizations requires 
some extra explanations. We argued earlier that to lead- 
ing order in 1/N, the spin polarization bubble has a 
regular piece, which comes from high fermionic energies 
and is absorbed into a bare static susceptibility, and an 
anomalous 7|w| term which comes from fermionic ener- 
gies which are smaller than u>. If we linearize fermionic 
dispersion near kp and integrate first over momentum, 
and then over frequency, we only obtain 7|u>|, the high- 
energy piece just does not appear. It is then tempting 
to conclude that the logarithmical vertex renormalization 
should yield u> log u> term but no (q— Q ) 2 log | q— Q | term. 
One should bear in mind, however, that for fermions with 
a linearized dispersion, the polarization operator is lin- 
early divergent at high energies. Due to divergence, the 
result of the momentum/frequency integration depends 
on which integration is performed first. Indeed, evalu- 
ating the particle-hole bubble for free fermions by first 
integrating over frequency and then over momentum, we 
do obtain, in addition to 7|a->| an extra regular piece, of 
order of the upper cutoff in the momentum integral. 

Consider now what happens when we evaluate the 
same bubble with the logarithmical vertex correction. 
Since we are interested in the terms which scale cither 
as ojL or (q — Q) 2 L, we can subtract the result for 
lu = (q— Q) 2 = 0. The remaining momentum/frequency 
integral diverges only logarithmically, and hence the re- 
sult of the integration no longer depends on which in- 
tegration is performed first. This implies that while "/to 
term in the "zero-order" spin propagator is an anomalous 
piece which cannot be obtained by expanding in u, the 
logarithmical corrections to the spin propagator are not 
anomalous and emerge due to the fact that d — z = 2 
case is marginal. [Alternatively speaking, these correc- 
tions can be found in a standard RG procedure which as- 
sumes that typical internal energies are much larger than 
the external ones.] This reasoning explains why (q— Q) 2 
and j\u>\ are renormalized in the same way. 

We now discuss why the series of logarithmical vertex 
corrections yield power-law momentum and frequency 



dependence of the renormalized vertex. Mathematically, 
this is a consequence of the fact that Ag scales as g 3 /"f- 
If 7 was independent of g, then dg would be proportional 
to 5 3 , i.e., dg 2 R /dL ~ {g 2 R ) 2 or l/{g R ) 2 oc L. This would 
result in a logarithmical behavior of gn. In our case, how- 
ever, bosonic dynamics is made by fermions, and 7 oc g 2 . 
This yields dg R /dL ~ <?#, i.e., logg/? oc L. 

The power-law renormalization of g R can be under- 
stood as follows. Suppose we perform calculations at ar- 
bitrary d. The vertex correction contains the product of 
two fermionic propagators and one bosonic propagator. 
As we discussed above, typical momenta in the fermionic 
propagators are much smaller than in the bosonic prop- 
agator. Then one can perform a 2D integration over 
dekdek+Q keeping the corresponding momentum depen- 
dence only in fermionic propagators. This momentum in- 
tegration transforms the two fermionic propagators into 
the constant densities of states. The vertex correction 
then reduces to the evaluation of the frequency integral 
of x(g, lo) in dimension d—2. The corresponding integral 
is logarithmically singular if d — 2 is the lower critical 
dimension in the problem. As the lower critical dimen- 
sion is the upper critical minus two, the vertex correction 
becomes logarithmically singular if d is the upper critical 
dimension in the problem. This is precisely what we have 
in our case. Now, at the lower critical dimension, there is 
no restriction as to what the series of logarithms should 
produce, and they can sum cither into the inverse power 
of a logarithm (as in Heisenberg magnets) or into the ex- 
ponent of a logarithm (as in XY model). It is therefore 
not surprising that we did find power-law renormaliza- 
tion of ga at the upper critical dimension. 

Note by passing that from general point of view, 
the spin-fermion model is characterized by two different 
space-time dimensions. One is obviously d + zb = d + 2, 
where zb = 2 is the bosonic dynamical exponent. The 
other space-time dimension is associated with the anoma- 
lous piece in E(k,0). This piece is produced solely by 
fermions (bosons can be treated as static). For fermions, 
the dynamical exponent zp = 1 (typical fermionic mo- 
menta and typical fermionic frequencies scale with each 
other) . Accordingly, the space-time dimension associated 
with the anomaly is d + zp = d+ 1. For d = 2, it is equal 
to 3 and is below the upper critical dimension. In partic- 
ular, one can easily check that for arbitrary d, the anoma- 
lous self-energy at £ = 00 scales as S(f2) ~ f2( d_1 )/ 2 and 
accounts for non- Fermi liquid behavior for all d < 3. 

We now return to Eq. J64J). Substituting L = 
(l/2)log(a>7/(7|u>| + (q — Q) 2 )), we obtain that at the 
QCP, the running coupling constant diverges as 
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R ccg (7M 



I log (7 k 



(q-Q)T 1/JV 

+ (q-Q) a )|-*«/e«. 



(66) 



Substituting this result into the spin polarization opera- 
tor we find that at the QCP, 
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X - i (q,c)cx^j ( 7 | w | + ( q _Q)^ 

ex [ 7 |u/| + (q - Q) 2 ]^ 1 | Iog(7|w| + (q - Q) 2 )| 



-Vy/3V X 



(67) 



This result implies that vertex corrections give rise to 
anomalous scaling dimension r\ for the dynamical spin 
susceptibility, 



X - 1 (q, W )(x( 7 |c| + (q-Q) 2 ) 1 -''. 



(68) 



Up to logarithmic corrections, rj = 2/N. For N = 8, 
Tj = 0.25 and 1 - r) = 0.75. 

We also emphasize that the renormalization of the 
bosonic propagator has little effect on fermions as the 
spin susceptibility appears in the fermionic self-energy in 
the combination with two side vertices, i.e., as gfjX^i <*>). 
We see from ( |67| ) that this combination retains its bare 
form (j\oj\ + (q— Q) 2 ) -1 . Accordingly, the only singular 
correction to fermionic self-energy comes from the loga- 
rithmical renormalization of vp: 



s(^)cx|fi| 1 / 2 |i gn| 



-1/2 



(69) 



B. a comparison with </f> 4 theory 

The results of the previous subsection imply that the 
2D fermionic system at the antiferromagnetic QCP is 
not described by the effective bosonic theory at the up- 
per critical dimension. This contradicts (f> 4 theory of the 
QCP developed mostly by Hertz [q5[ and Millis |5g] (for 
recent developments, sec |9l]]). Hertz and Millis (HM) 
departed from the same spin-fermion model as we study, 
but integrated out fermions in the effective action and 
obtained the effective Lagrangian for bosonic (spin) de- 
grees of freedom. They conjectured that the expansion of 
the action in powers of the bosonic variables S qtl) yields 
<p A theory with the dynamical exponent z = 2. 

S e ff = T^{r + q 2 - iw)Sq )U S-q,- u 
aj,q 

+b 2_, Sq l ,u 11 S q2 ^ 2 S q3 ^. i S qi ^ i 8(T, i qi) 5(SjWj) + .... (70) 

where q measures a deviation from (tt, tv), r is a distance 
from the critical point, and b is a constant. At d = 2, the 
<^> 4 theory is at the upper critical dimension [d + z = 4), 
and the bosonic self-energy is marginally irrelevant. This 
implies that for r = the self-energy only adds fractional 
power of log(max(w, q)) to the bosonic propagator. One 
can check this by explicitly computing the self-energy 
order by order in b. 
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FIG. 16. a) Diagrammatic representation of the four-fold 
interaction vertex between spin fluctuations, b) and c) Two 
ways of contracting the external legs of the four-fold vertex. 

The bare bosonic propagator in the HM theory is the 
same as our N = oo result. We, however, dispute the 
assertion that b can be approximated by a constant. In- 
deed, if this approximation is valid, then the fist-order 
correction in b involves one bosonic loop and is infrared 
finite. We, however, found in the previous subsection 
that it is logarithmically singular. 

The four-boson vertex is diagrammatically given by 

At external u>i 



the graph presented in Fig. 16 



Qi = 0, 



it becomes (we skip spin indices for a moment) 



b = NfT ]T G 2 (q', W ')G 2 (q' + QV). (71) 

q',0?' 



Since 



J de' q G 2 {c(,uj')=0, 



(72) 



this vertex vanishes if one linearizes the fermionic disper- 
sion near the Fermi surface. The inclusion of the curva- 
ture of the fermionic dispersion (i.e., nonunivcrsal, lattice 



effects) makes the momentum integral in (71) and hence 
b finite |5q |. The fact that b is produced by high-energy 
fermions then implies that it is insensitive to any modifi- 
cations of the low-energy fermionic dynamics near quan- 
tum criticality and can be treated as a constant even at 
the critical point. 

The potential problem with this reasoning is that al- 
though the momentum integral in ( |72| ) indeed vanishes 
for a linearized dispersion, the integrand contains a dou- 
ble pole. At finite external tOi and qi, the double pole 
splits into two single poles which in a small frequency 
range of width loi are located in different half-planes. 
The integration over this range does not produce a small- 
ness in u>i as the two poles are only separated by a 
max(u>i,VFqi)- This implies that even for linearized 
fermionic dispersion, the value of b is generally finite and 
depends on the ratio of Wj and vf<^. 

The evaluation of the four-boson vertex for finite ex- 
ternal frequencies and momenta is straightforward. For 
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simplicity, we restrict with free fermions. One can eas- 
ily verify that the inclusion of the fermionic self-energy 
does not affect the results below. Introducing the exter- 
nal momenta and frequencies as in Fig 16 a,. we obtain 
(see Appendix C, Eqn. (p89h) 



b = 



Ng \UJ2 + U)\ + \i02 — lu\ — \uj\ + u>\ — \oJ\ — U)\ 

8ttv x v v [i(ui +lu 2 ) - € qi+q2 \[i{u>i - lo 2 ) - £Q+ q 



1-92J 



(73) 



We see that generally b scales inversely with the incom- 
ing frequency, but the coefficient depends on the relative 
magnitudes of external loi and vfQi- This in particular 
implies that the bosonic self-energy depends on which ex- 
ternal legs of the four-boson vertex are contracted. There 
are two topologically different self-energy diagrams to 
first order in b. They are presented in Fig. Rq b and 
c. If b was independent of momentum and frequency (as 
is assumed in the HM theory), both diagrams would be 
equivalent and yield b J d 2 gcL;x(q, c<j). This integral is 
nonsingular in 2D and only accounts for a nonessential 
renormalization of £ . In our case, the result is different 
as we already know: the first diagram is equivalent to 
adding fermionic self-energy to the particle-hole bubble 
and vanishes if we restrict with linearized fermionic dis- 
persion near the Fermi surface. The second diagram is 
equivalent to adding a vertex correction to a particle-hole 
bubble and is logarithmically singular. 

For completeness, in Appendix C we explicitly com- 
puted both diagrams by choosing the appropriate exter- 
nal frequencies in b, and convoluting b with the spin prop- 
agator. This computation indeed yields the same results 
as we just described, namely the diagram in Fig. pL6| b 
vanishes while that in Fig . |l6| c yields the same logarith- 
mical renormalization as (|65|). 

We see therefore that the explicit computation of the 
first-order bosonic self-energy in the effective </> 4 model 
yields a logarithmically singular piece which comes from 
a singular momentum/frequency dependence of the <f) A 
vertex at smallest energies. This singularity is missed if 
one restricts with a constant four-boson vertex, as it is 
assumed in the 4 theory. 

A very similar computation has been performed by 
Lercher and Wheatley p3J ]. They also argued that the 
four-boson vertex is singular at vanishing external mo- 
menta and frequency, and depends on the ratio between 
momentum and frequency. Our result for the vertex 
slightly differs from theirs, but the conclusions are iden- 
tical. 



C. summary of Sec. M 

Let us now summarize our 1/N results. 



i We found that the singular corrections to the Fermi 
velocity cause nesting behavior at hot spots at 
£ = 00, but have little effect on the spin dynam- 
ics. These corrections are, however, only logarith- 
mically singular at £ — > 00. The corrections to the 
vertex on the other hand do not affect Fermi ve- 
locity and fermionic self-energy, but give rise to a 
nonzero anomalous scaling dimension 77 = 2/N of 
the bosonic propagator. 

ii We demonstrated that the co-existence of the loga- 
rithmical renormalization of the Fermi velocity and 
the anomalous dimension in the spin channel is due 
to the fact that the corrections to the spin-fermion 
vertex are structurally equivalent to the renormal- 
ization of the order parameter in dimension d— 2. In 
the latter case, the logarithmical corrections are not 
necessary summed up into a series of logarithms. 

iii Our results disagree with purely bosonic 4 the- 
ory of the QCP. We argued that the difference is 
due to the fact that in our model, the four-boson 
vertex is not a constant but rather has a complex 
momentum/frequency dependence at the lowest en- 
ergies. This is another indication of the presence of 
anomalies in the spin-fermion model. We explicitly 
demonstrated how our results can be obtained in 
the effective <fi theory with full four-boson vertex 
made by fermions. 

iv We also argued that from general perspective, the 
spin-fermion model has two different space-time di- 
mensions, d + 2 = 4 and d + 1 = 3. The first is as- 
sociated with the RG-type physics, and the second 
is associated with the anomalous piece in the self- 
energy which comes from low energies and cannot 
be obtained within RG procedure. For 2 < d < 3, 
the bosonic dynamics at criticality is Gaussian, but 
fermionic self-energy scales as E(J7) oc fi( d_1 )/ 2 and 
accounts for non-Fermi liquid behavior at the crit- 
ical point. 



VI. LARGE N THEORY AT FINITE T 

So far, our analysis was restricted to T = 0. We now 
discuss how the fermionic self-energy and the spin polar- 
ization operator are modified at finite T. 

A conventional wisdom would be that at finite T, 
the system in the quantum-critical regime possesses ui/T 
scaling, i.e., fractional powers of frequency are replaced, 
at small w, by the same powers of temperature. We will 
see that this effect, indeed, takes place. We will also find, 
however, that the fermionic self-energy in the quantum- 
critical regime contains an extra term which does not fit 
into lu/T scaling. The strength of this extra term de- 
pends on the order of limits A — > 00 and N — > 00. 
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As before, we first present the results of the compu- 
tations at N — oo and then extend the theory to finite 
N 



A. N — » oo limit 

1. spin polarization operator 

We begin with the spin polarization operator II(Q,w) 
Replacing the integration over Matsubara frequencies in 
Eq. (pq) by a summation as 



2tt "* ^ 



we obtain 



J n \ 



^(n n ) 2 (Q n +UJ m ) 2 ) ' 



(74) 



(75) 



where, as usual, uj m — lixTvn, f2„ = ■nT{2n + 1). The 
summation over n is elementary and yields 



n„ 



■kT 

\0Jm\ 
U> s f 



Z^ 



{1 - sign ([2n + 1] [2(n + m) + 1])} 



(76) 



We see that at finite T, the spin polarization operator 
preserves exactly the same form as at T = 0. In partic- 
ular, n(Q,0) = 0, i.e., at finite T, the spin correlation 
length does not acquire an extra correction from low- 
energy fermions. 



2. fermionic self-energy 

As the spin polarization operator does not depend on 
T, the fermionic self-energy is still given by Eqs. ( 122 ) 
and (p3). At N — > oo we have 



E(kh 8 ,fl! m ) = iirTX ^ 



signo>„ 



1 



k 



= iirTX 



\m\ 



1 + i £^ 



sign£! r 



(77) 



The extension to momenta away from a hot spot proceeds 
in the same way as at T = 0: A and uj s f are just replaced 
by X(k) and w s /(fc) given by (fl7|). The full analysis of 
self energy as function of frequency and temperature is 
performed in the Appendix D. For |O m |, T <C u> s f(k), the 
self-energy obviously has a Fermi-liquid form 



s(k,fi m ) = ? A(fc) [n m + 



4u, f {k) 



signTl r 



(78) 



For T <C ^sf(k) and arbitrary fi m , the thermal piece in 
the self-energy preserves T 2 dependence, but the prefac- 
tor depends on frequency: 



£ T (k, O m ) 



■ ^T 2 

%-. rrrr 



4iJ a f(k) \0J s f(k) 



a, 



signw + 0(T A ) 



(79) 



where F(0) = 1, and F(x > 1) — > 1/3. The full expres- 
sion, including T 4 terms is presented in ( |213j ). 

In the opposite limit |Q m |,T 3> to s f(k), the system 
behavior is quantum-critical. In this limit, the contri- 
butions to the sum in Eq. ([77j) from all n ^ m can be 
explicitly expressed in terms of Ricmann Zeta function 
(see Eq. (203)). Combining this with the n = m contri- 
bution from classical, static spin fluctuations, we obtain 



£(k,O m ) = mTX{k) 



c 



(80) 



where £ is a Zeta function. For m 3> 1, the second term 
scales as y/Tm, i.e., as y/u. 

In real frequencies, the imaginary part of the self- 
energy takes the form (see Appendix D) 



S"(fi) = -Im J 



dio 



1 



s /(n/2T,w/2T) (81) 



w,/ 



where f(x,y) — tanh(x — y) + l/tanh(y). In the Fermi 
liquid regime, we obviously have £"(fi) oc T 2 , with fl 
dependent prefactor, as in (|79[). In the quantum-critical 
regime, 



?'lkMT) = ,TMk) !- (^ 



1/2 



Dig 



(82) 



The scaling function D(x) is given by (214). In the two 
limits, D(x > 1) = s/x, and D(x -C 1) = -1.516 + 
0.105x 2 . 

The real part of S(k, Q) is obtained from ( p2j ) by 
Kramcrs-Kronig transformation. It contains the same 
y/T dependence as (82) with another scaling function of 
fi/T, but does not have a AT term. This can be eas- 
ily understood as static spin fluctuations account for a 
scattering with a finite transferred momentum and zero 
transferred frequency, and therefore act as impurities. 
Not surprisingly, their contribution to the self-energy is 
purely imaginary. 

We see from Eqs. (p0|), ( JB2[ ) that the piece in the self- 
energy which comes from nonzero bosonic Matsubara fre- 
quencies (i.e., from n ^ m in (|77J)) yields a conventional 
scaling form of S(k, il, T) in which the vf2 behavior at 
ft ^$> T is collaborated by y/T dependence at T 3> fl, 
and the prefactor does not depend on the spin correla- 
tion length 
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On the other hand, the n = m term in the sum in < Jn\ ) 
yields a linear in T contribution to the self-energy with 
the prefactor which scales as £ and diverges at the QCP. 
This term obviously violates the scaling. This violation 
is a direct consequence of the fact that in our theory 
£(T) = £(T = 0). For the scaling behavior, one would 
needed £, 2 {T) — £ 2 (0) oc T. We discuss this issue in a 
separate subsection below. 
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The results for the fermionic self-energy at 



Finally, at 51 — 0, E"(T) has a simple form which 
interpolates between Fermi liquid and quantum-critical 
regimes. Taking the 51 = limit in (pyl) we obtain 



E"(fi = 0) = ttTAL ( — 



where 



LfcO 



1 

7r.r 



dylm 



VI - *2/ sinh(y/x) 



(83) 



(84) 



At x > 1, L[x) ~ 1, i.e., E"(51 = 0) « ttTA. In the 
opposite limit of small x, L(x) = 7rx/4, and we recover 
the Fermi liquid behavior. 

In Fig. [17|we plot the iV — * oo result for the self-energy 
at finite T. We see that at small frequencies, ImE(fe, 51) 
is linear in T starting already from T > 0.5cj s /. At 
higher frequencies, particularly in the region 51 3> T, 
ImE re *(k, 51) is linear in 10 with almost T— independent 
slope. A finite temperature only gives rise to an extra 
constant term in ImE re *(k, 51). 



B. beyond the N — » oo limit 

We now consider to which extent we can trust the 
N — * oo results for the fermionic self-energy and the 



spin polarization operator, when we apply them to the 
physical N = 8. We will show that there are two prob- 
lems at finite N which are only partly related: one is 
how strong are vertex corrections, the effect of keeping 
the self-energy of intermediate fermions, and the momen- 
tum dependent piece in the fermionic self-energy (i.e., the 
renormalization of the Fermi velocity) , and another one 
is to what accuracy one can factorize the momentum in- 
tegration in deriving the equation for the fermionic self- 
energy. [By factorization we mean that in the diagram 
for fermionic E(fc,0) one integrates over the momentum 
component along the Fermi surface in the bosonic prop- 
agator and over the momentum component transverse 
to the Fermi surface only in the fermionic propagator. 
The resulting self-energy is then the convolution of the 
effective ID "local" spin susceptibility and the fermionic 
density of states] . 



1. accuracy of the factorization of the momentum 
integration at finite T 

Let's first discuss the second issue. Consider the 
second-order diagram for the fermionic self-energy. At 
T = 0, the factorization of the momentum integration in 
this diagram, or, equivalcntly, the possibility to neglect in 
the spin susceptibility the momentum component trans- 
verse to the Fermi surface, is well justified for all rele- 
vant frequencies as the correction to the self-energy from 
keeping the momentum dependence in the susceptibility 
has an extra power of 51 in the Fermi liquid regime and is 
small by 1/N 2 in the quantum-critical regime. For T^0, 
analogous result holds for the n y= m contribution to the 
self-energy, and the proof almost exactly parallels the one 
for T = 0. The situation is, however, different for the 
n = m contribution to ( |77| ) . This term comes from static 
spin fluctuations, and hence typical (q — Q) x ~ \ui\/vf 
has to be compared to £ -1 rather than to (^lu) 1 / 2 (x axis 
is chosen transverse to the Fermi surface at k/j S + Q). 

The new frequency scale u,f£ _1 is related to to and u> s f 
as 



vfC 



N 



NX 



v x v y 12A 



J sf 



V X Vy 



(85) 



We see that the location of this scale compared to Q 
depends on the ratio N/X. In the formal N — ► oo limit at 
finite A, «f£ _1 S> w. Then, even at w ~ w, \u>\ ^C vpt; _1 , 
and the factorization of momentum integration is well 
justified. Alternatively, however, one can keep N finite 
and study the system behavior at A — ► oo and finite T. 
In this limit, i>_f£ _1 *S u>, and the corrections to the 
formal N — oo result for the static, thermal self-energy 
are relevant at frequencies \u>\ ~ uj(N/X) <C u). In reality, 
for N = 8 and A > 1, vp^ 1 < ui, i.e., the modification of 
the static piece in the fermionic self-energy is a modest 
effect for ui < u>. Note also that this modification is 
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only relevant in the quantum-critical regime as in the 
Fermi liquid regime, \u>\ < ui a f, and hence \uj\/(vf£ _1 ) < 
w-z/^r 1 ) = 0(1/N). 



this limit, the overall factor in the fermionic self-energy 
scales as Aup£ _1 ~ 5, i.e., it no longer depends on £. This 
is indeed what one should expect in the truly quantum- 
critical regime. 




FIG. 18. Function f(x) defined in (|89|). Dashed line is the 
asymptotic form 21og2a;/(7ra;), valid at large x 

The exact form of S(O m ) at a hot spot, can be ob- 
tained from the general formula for the self-energy, Eqs. 
( f22| ) and (23), extended to a finite T. We separate the 
contributions from dynamical and static spin fluctuations 
and write E(fi m ) = £ st (£l m ) + ^dyn(^m)- As we said, 
Tidy n can be safely computed at N = 00 as the correc- 
tions scale as regular powers of 1/N. The result is the 
second term in Eq. ([77]): 
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^dyn(^m) = JTrTX £ 4 



1 



=signfi r 



(86) 
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In the quantum-critical regime this reduces to the second 
term in Eq. @ 



^dyn \* 'r. 



m 



1/2 



<(i.i 



(87) 



The static contribution has to be evaluated without ap- 
proximations and reduces to 



■2 S < (S2 r , 



i ATtt/ fi 



where f2 r , 



l^ml/^F^ * an d 

2 cos -1 fi r 



/(a») = - 



1 - ni 



(88) 



(89) 



The limiting values of f(x) are /(0) = 1, /(l) = 2/tt, 
/(a; 3> 1) sw 21og2a;/(7rx). The plot of f(x) is presented 
in Fig. [18]. We see that £ s t(O m ) coincides with the 
N = 00 result i7rA only at vanishing Vt m . At finite fi m , 
it is smaller than for N — 00, and for f2 m 3> 1, £ s t(O m ) 
decreases as l/ft m , up to logarithmical corrections. In 
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FIG. 19. The modified second-order fermionic self-energy 
- the sum of Eqs. (|9C|) and (J92J). It includes the full 
contribution from thermal, static spin fluctuations, and the 
N — > 00 contribution from dynamical spin fluctuations. 
We used v x — v y — vf/V%, and A = 1, in which case 
vf^ 1 = 4^/3 = 16w,//3. 

In real frequencies, £"(k, ls ,fi) = £' s ' t (fi) + E'^JO), 
where YI' d (fi) is the same as in (|8l| ) but with fi(x, y) = 
tanh(cc — y) + 1/ tanh(y) — 1/y instead of f(x, y): 



-'dyn 



(SI) 



A 



-Im 



(luJ 



=/i(n/2T,w/2T). (90) 



For oj > o> s /, this reduces to the second term in (p2|): 



-:U( () '-(^f) 



1/2 






(91) 



where -D(x) is given by (214). 

The static contribution is obtained by transferring 
(pq), (p9[) onto real axis: 



£" t (0) = 2TAIm- 



log 



o 2 



v/i + fi 2 



\/l + n 2 



(92) 



For Cl < 1, £" t (ft) oc (T/fijlogO/VF*;- 1 . We see that 
the prefactor depends on £ only in the logarithmical term. 
This is indeed the same result as we found earlier by eval- 
uating the thermal self-energy in Matsubara frequencies. 
In Fig. O, we plot E"(k/j S ,0) obtained by combin- 
ing the the N — 00 form of S^n (Eqn. (pQ)) and the 
modified form of E" t (Eqn. (J92)). Comparing this figure 
with Fig. O where we plotted the formal N — 00 result, 
we observe that the main difference is in the tempera- 
ture variation of the self-energy at high frequencies. At 
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N = oo, the plots at different T differ by a constant. In 
the modified expression, the thermal piece in £"(k, 17) 
decreases with frequency, and the curves at different T 
progressively come closer to each other. 
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FIG. 21. The photoemis- 

sion intensity 7(17) oc JmG(k,fi)ni?(fi) for various e^. For 
definiteness we used A = 1.7 and T = 0A3u s f. 



FIG. 20. The full second-order result for the E" (fe.fi) 
obtained without factorizing the momentum integration in 
both static and dynamics pieces in the self-energy. In the 
insert the same data for T — 1.3 (solid line) is compared with 
the approximate form of the self-energy from Figjlfl (dashed 
line). vf£ -1 is the same as in Figjl9l 

To verify how accurate is our "modified N = oo" form 
of the self-energy we present in Fig. pfl the result for the 



full £", which we obtained in Appendix D (Eqn (218) 
with T? R A = 17 — ui ± iS) without factorizing the momen- 
tum integration in both, static and dynamical parts of 
£". Comparing Figs. |l^ and [2(] we find that the full 
result is rather close to our approximate one - the only 
difference is that now he curves at various T progressively 
come close to each other with increasing frequency, and 
cross only at very high frequencies. 

In Fig. ^l] we present the set of the results for the pho- 
toemission intensity /(17) oc iroG(k, !7)7i^ (17) for vari- 
ous efc. They were obtained using the spectral functions 
from Fig. |2fj. For definitness, we used A = 1.7 and 
T = 0.43oj s /. Comparing these forms with the corre- 
sponding results at T = (see Fig. |l2J) , we see that a 
finite temperature gives rise to an additional broadening 
of the photoemission intensity at small 17 ~ w s /, and 
plays virtually no role for 17 3> <jj s f. 



2. the role of the self-energy of intermediate fermions 

In the above analysis we restricted with the second or- 
der self-energy diagram. Now we consider what happens 
if we keep the self-energy of intermediate fermions (i.e., 
evaluate the self-energy diagram as in Fig. (Q)). Recall 
that at T = 0, the self-energy of intermediate fermions is 
not small in 1/JV but nevertheless doesn't play any role as 
the momentum integration can be factorized, and S(w) 
is expressed in terms of the density of states of intermedi- 
ate fermions. For the linearized fcrmionic dispersion, the 
latter is independent on the self-energy and is the same 
as for free fermions. This implies that the full expres- 
sion for the self-energy coincides with the second-order 
result. At finite T, this is no longer the case as the mo- 
mentum integration in the static piece in £(17) cannot be 
factorized. 

In principle, the self-energy of intermediate fermions 
should include both T = and temperature contribu- 
tions. We will show below that the temperature piece 
in the self-energy for intermediate fermions becomes rel- 
evant for T 3> Nto s f when also vertex corrections and 
the renormalization of e& cannot be neglected. To sep- 
arate thermal and quantum effects, we focus first on 
T <C Nu> a f when the vertex corrections and the correc- 
tions to efc are the same as at T = 0, i.e., are small in 
1/N, At these temperatures, the self-energy of interme- 
diate fermions can be approximated by its T — value. 
The latter, in turn, is well approximated by the N = oo 
result, E T =o(17) = 2A17/(1 + (1 - ifl/u.f) 1 / 2 ). 
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FIG. 22. The full result for the E"(fc, Q) (Eqn. |218|). The 
form of £"(fi) at low temperatures is close to the second-order 
result plotted in FigMJ. The initial decrease of E"(0) at 
higher T is probably the artifact associated with the neglect 
of vertex corrections. i>f£ _1 is the same as in Fig.n9l 

The full expression for the static part of E" is the same 
as in Eq. (|92|), but now fl = (0 + Et=o(^))/vf£ _1 - 
In the quantum-critical regime, Q < Q, Er = o(fl) ex- 
ceeds $7, and Cl becomes iil/Q where ft — (vp^ 1 ) 2 /u> = 
N 2 uj s f(vp/v x Vy)/36. Then, when Q, < Q < Q, the 
static piece in the self-energy behaves as E" t (f2) ex 
TN^ui/fl log fl/N 2 u} s f. We see that up to a logarith- 
mical prefactor, E" t (f2) acquires the same scaling form 
as the dynamical piece in the fermionic self-energy. In 
other words, at high enough frequencies, we eventually 
recover il/T scaling. 

In Fig. p2 we plot the full expression for E" obtained 
in the Appendix D, (Eqn. (218)) without factorizing 



the momentum integration in both static and dynamic 
parts and keeping the T = 0, N — oo self-energy for 
intermediate fermions. Comparing this figure with the 
second-order result, Fig.gfJ, we see that at low T, the two 
forms of E"(f2) are close to each other, as they indeed 
should be, as at T — keeping the self-energy of inter- 
mediate fermions only yields regular 1/N corrections. At 
higher temperatures, £'P" me in Fig. [2^ differs from the 
second-order result at high frequencies, as we found an- 
alytically, but it also initially decreases with frequency 
due to a stronger downturn renormalization of the static 
piece then in the second order expression. 

This initial decrease poses a dilemma in selecting the 
theoretical form of E" for experimental comparisons. 
From theory perspective, the T = self-energy of in- 
termediate fermions is not small in 1/N and therefore 
cannot be neglected in evaluating E" t (0) for which the 
momentum integration is not factorized. The restoration 
of the scaling behavior at finite T due to this self-energy 
is also a good sign. On the other hand, one can include 



the self-energy of intermediate fermions but neglect ver- 
tex corrections and the momentum dependence of the 
self-energy only as long as T < Ncu s f (see below). Up to 
which temperature this can be safely done for a physical 
N = 8 is a'priori unclear. It looks to us that the initial 
drop in E"(f2) in Fig. 02 at T > w s / is an artifact of 
overestimating the reduction of E" t (f2), and should be 
compensated by vertex corrections. This reasoning im- 
plies that for T > u) s f and N — 8, the second-order result 
for the self-energy, Fig. EOl may be closer to reality than 
Fig. p2l which explicitly includes the self-energy of in- 
termediate fermions but neglects vertex corrections. In 
Sec [X| we will be using the second-order result for the 
comparisons with the photoemission and optical data. 
We stress, however, that for T < uj s f, the difference be- 
tween the two results is minor for all experimentally rel- 
evant frequencies. In particular, we verified that optical 
conductivities obtained using the second-order and the 
full self-energy are almost identical for T < ui s f. 

We now return to the theoretical discussion and con- 
sider what happens at T > Nuj s f. As we already said, 
this temperature range is of little use for experimental 
comparisons with cuprates, at least near optimal dop- 
ing. However, it is important for the understanding of 
what is the correct form of the self-energy at a finite T 
at the QCP, when ui s f — > 0. Recall in this regard that 
so far we still have E' s ' t (£7 = 0) = ttTA cx £. On the 
other hand, from general scaling arguments pjj we ex- 
pect that Eg t (fi) in the quantum-critical regime should 
be independent on £. 



3. the renormalization of the Fermi velocity 

We begin with the velocity renormalization. To esti- 
mate it, we need to analyze the strength of the static 
piece in the self-energy ^k,n m at k ^ hh s and e^+Q ^> 
fl rn . Performing the same computation as in the previous 
subsection we obtain after simple manipulations, 



E(k, Si m ) — E(k/ ls , \l r 



2T\(k) log [(e k+Q + sJel+Q + vU' 2 ) /(«**" ')] (93) 

At small ejc+Q , expanding under the logarithm and using 
the fact that e/c = v x k x +v y k y and ek+Q = ~v x k x +v y k y , 
we obtain 

v* = v x (l + 8) 



5 = 



v v (l-S) 

2v x v v 



T 



3^ 



Nuj 



«./■ 



(94) 



We see therefore that velocity renormalization becomes 
relevant at T ~ Nu> s f, as we anticipated. 
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4- vertex corrections 

The result for the vertex renormalization parallels the 
one for the self-energy. We can split Ag/g into the contri- 
butions from static and dynamic spin fluctuations. The 
contribution from dynamical spin fluctuations scales as 
(l/7V)logT and is small for T < uj that we are inter- 
ested in. The contribution from static spin fluctuations is 
larger due to the absence of the large juj term in the spin 
susceptibility. Integrating over momenta of two fermionic 
propagators in the vertex correction diagram in the same 
way as at T = 0, we obtain 



(Ag 
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IT 



Nuj 



tan — 



of 



(95) 



This expression is valid for small external u> < vf£, ■ 
At larger frequencies, (Ag/g) st decreases and eventually 
scales as (gT/uj 2 MlogT\. 

We see from (g5|) that vertex corrections also become 
relevant at T > Nuj s f- 



5. the thermal self-energy of intermediate fermions 

Finally, we show that at T > Nu) s f we also have to 
keep the thermal piece in the self-energy for intermediate 
fermions, The easiest way to estimate this temperature is 
to consider the limit Q = 0. At vanishing frequency, the 
T = piece in the fermionic self-energy vanishes, and we 
are left with only £ st (T). At N = oo, S at (T) = irTX. 
At finite N, this expression is valid if £ st (T) < Vp£~ l , 
or T < Nu> s f. At larger T, the thermal piece in the self- 
energy for intermediate fermions cannot be neglected, at 
least at small frequencies. 



C. self-consistent theory 

We are not familiar with any controlled way to com- 
pute the static piece in the self-energy at finite TV, finite 
T, and £ — ► oo such that T ^> Nu> s f. As an exercise, 
we performed the FLEX-type computation of £ st (f2,T) 
at vanishing Q. Specifically, we included the full static 
fermionic self-energy into the self-energy diagram, but 
neglected vertex corrections and the renormalization of 
the Fermi velocity. We solved the self-consistent equa- 
tion, obtained S st (T), and then re-calculated the vertex 
and velocity corrections with the new fermionic propaga- 
tors. We found that the self-energy no longer diverge at 
the critical point, and that the vertex and velocity cor- 
rections are at most logarithmical even when T 3> Nu> a f . 
While this analysis is by no means a controllable one, it 
describes the physically appealing effect that the diver- 
gence in the self-energy and in the vertex are eliminated 



by thermal dressing of intermediate fermionic quasipar- 
ticles. 

The FLEX-type computation of £ st (T) reduces to the 
solution of the self-consistent equation in the form 



%t(T) 



3T 
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where S st = £ st /(vj?£ 1 ) The renormalization of the 
Fermi velocity and the vertex correction depend on 

t st (T) as 
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For simplicity, we set v x = v y = vp/\/2. The integral in 
(30) is explicitly evaluated in the Appendix D (see Eqn 
( 21 6| ) ) . The integration in ( J97J ) can also be performed 
exactly, but the results are too lengthy to present here. 

At T -C Nco s f , the self-energy in the denominator can 
be neglected, and we obtain £ st (T) = ttTX. For these T, 
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In the opposite limit, solving (Bq) with logarithmical ac- 
curacy, we obtain 



S Bt (T) 



= *(™\ 
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Analyzing these results, we find that self-consistent so- 
lution of the thermal problem yields the results which 
are qualitatively similar to that at T = 0. Namely, 
the thermal piece in the fermionic self-energy scales as 
T 1 / 2 , up to a logarithmical prefactor, the velocities at 
\z.hs and k/, s + Q logarithmically tend to become antipar- 
allel to each other, and vertex correction remains finite 
at criticality. The self-consistent solution indeed does not 
capture full logarithmical physics, nevertheless it clearly 
demonstrates that higher order diagrams eliminate the 
unphysical O(A) divergence in the fermionic self-energy 
and yield VT behavior at the QCP. Alternatively speak- 
ing, the self-consistent, finite N solution recovers uj/T 
scaling in the fermionic variables. From this perspective, 
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the only real difference between the self-energy at T = 
and T ^ in our large N analysis is that at T ^ 0, the 
width of the intermediate regime where E" oc T scales 
with N [u a f < T < Nu a f), while at T = 0, the width of 
the analogous linear regime in T,"(uj) is, strictly speak- 
ing, of order 1 around T ~ uj s f. Curiously enough, for 
physical N = 8, the widths of the linear regimes in uj and 
in T are comparable as the one in uj occurs in between 
ui a f and 6 — 8oj s f- Whether this is more than just a co- 
incidence is unclear to us. In any event, a more accurate 
treatment of thermal fluctuations is clearly called for. 

Note, however, that while uj/T scaling in the fermionic 
self-energy is probably recovered in the self-consistent so- 
lution, still, there is no uj/T scaling in the dynamical 
spin susceptibility. Indeed, the thermal correction to the 
spin-fcrmion vertex remains finite at criticality, while at 
T = 0, the corrections were logarithmical and eventually 
yielded an anomalous dimension rj. While, as we said, 
logarithmical physics is not precisely captured by self- 
consistent solution, there are no reasons to expect T 1 ^ 71 
behavior as the anomalous power at T — emerged due 
to the presence of anomalies which are purely quantum 
effects. 



D. a temperature variation of £ 

Finally, we discuss is some more detail the temper- 
ature variation of the correlation length. The issue is 
whether there exist a universal temperature correction 
to the correlation length. We recall that in our theory 
£(T) = £(T = 0). This independence of T is the conse- 
quence of the fact that in the computation of the spin 
polarization operator, we linearized the fermionic disper- 
sion near the Fermi surface (see Appendix C). Expanding 
fermionic energies beyond the linear order in k — kp, one 
indeed obtains thermal correction to £~ 2 , but the lat- 
ter is nonuniversal and is of order (T/A) 2 , where A (our 
upper cutoff) is comparable to the fermionic bandwidth. 
The nonlincarity of the fermionic dispersion also gives 
rise to a constant term in four-boson vertex b ex 1/A. 
For a constant b, the first-order bosonic self-energy in 
2D, £;, ~ b T yi r-l J d 2 qx(<i, uj) is logarithmically singular 
at T ± §|j E 6 ~ bTjd 2 q/{q 2 + ^ 2 ) ~ 6Tlog£. 
Cutting the divergence at criticality by 3D effects, we do 
find that at the QCP, £ _2 (T) ex bT. This correction is 
larger than T 2 term, but still it is nonuniversal and scales 
inversely with A. 

Another widely used approach to the quantum-critical 
behavior near the antiferromagnetic instability explores 
the assumption that x(Q) w ) satisfies a temperature inde- 
pendent sum rule constraint j d 2 qx(<i, w) = const p9[ . 
One can easily demonstrate that for z — 2 dynamics 
at low-energies, the constraint implies that at the crit- 
ical point, £ -2 = aT with a model independent, uni- 
versal a [jn|[nj. This contradicts the analysis based on 



the spin-fermion model. The most likely reason for the 
discrepancy is that the spin-fermion model is a finite U 
version of the Hubbard model, while the hard constraint 
is valid only at infinite U. Indeed, if U is finite, the no- 
double occupancy constraint is not exact, and there is a 
finite probability for two electrons to occupy the same 
site. The spin of this doubly occupied site can be either 
zero or one, i.e., the "length" of the spin at a given site 
is not fixed. From this perspective, it is likely that the 
sigma-model approach captures the physics at very large 
U, while our theory describes the physics at moderate 
U < A, when lattice effects (such as bT correction to 
£~ 2 ) are irrelevant. 



E. summary of Sec. V] 



We now summarize what we obtained in this section 

i We found that the spin polarization operator at 
finite T has precisely the same form as at T = 0. 
There is no universal T dependent correction to the 
spin correlation length. 

ii We found that the N = oo theory for electronic 
self-energy at finite T is rather peculiar: the T 
dependent corrections from scattering on dynam- 
ical spin fluctuations (the ones with finite Matsub- 
ara frequency) have the same dependence on T as 
of frequency, i.e., they obey Cl/T scaling. In the 
quantum-critical regime, this correction scales as 
y/Tu>. At the same time, the piece in the self-energy 
from thermal scattering on static spin fluctuations 
yields impurity- like contribution E(T) = iirXT and 
still scales with the correlation length £. The still 
presence of £ in the quantum-critical form of the 
self-energy is a serious artifact of the N = oo the- 
ory at a finite T. 

iii We demonstrated that this artifact disappears 
when calculations are performed at a finite N. 
First, we obtained the full form of the self-energy 
in the second-order perturbation theory, when in- 
termediate fermions are treated as free particles. 
We found that at large frequencies, the static piece 
in the self-energy becomes independent on £, up 
to logarithmical factors. Then we demonstrated 
that for T <JC Nuj s f, one can obtain in a control- 
lable way the full E"(J7), by explicitly including 
the T = self-energy of intermediate fermions. At 
frequencies which exceed fl = N 2 oj s f(v' F /v x Vy)/36, 
this self-energy obeys ft/T scaling. We also found 
that for T > w s/ , this E'P" me (£7) has a dip at low 
frequencies which is likely an artifact of overesti- 
mation of the reduction of the static part of the 
self-energy at small but finite fi. We argued that 
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the second-order result for the self-energy may be 
closer to reality at these T. 

iv We found that at small tt, the controllable calcula- 
tions of fermionic S(f2) definitely break down when 
T exceeds Nu> s f. For these temperatures, we per- 
formed FLEX-type calculations and found that the 
self-energy recovers \ITlo form, up to logarithmic 
corrections. These calculations are not controlled 
as the vertex renormalization and the renormaliza- 
tion of the Fermi velocity are 0(1). We argued, 
however, that these renormalizations affect loga- 
rithmical factors but not VT dependence. 

v Our finite N analysis indicates that at finite N, 
the fermionic self-energy likely possesses il/T scal- 
ing. At T, $7 > u> s f, E(fi, T) is linear in Q and 
T at intermediate energies and crosses over to 
max(fi 1 / 2 , T 1 / 2 ) behavior at high energies. The pe- 
culiarity of large N is that the width of the interme- 
diate range along T axis is by a factor of N larger 
than that along frequency axis. 

vi We found, however, that there is no ui/T scaling 
in the spin susceptibility, i.e., tu 1 ^ 71 behavior at the 
QCP at T = is not collaborated by the same tem- 
perature behavior of the spin correlation length. 
We argued that this difference is due to the fact 
that the anomalous dimension r\ emerges due to the 
presence of anomalies which are quantum effects. 



VII. OPTICAL CONDUCTIVITY 

In this section we use the results for the self-energy and 
study in detail the behavior of the optical conductivity. 

The diagonal component of the optical conductivity is 
given by [H 



lo 1 1 i r 

(Jii{iu) = CTl(w) + l(J2{u) = -T-7T- d9 

47T ZTV J 



IUM) 
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(101) 



where uj v i is the plasma frequency, 9 is the angle along the 
Fermi surface, and II CT (0, lo) is fully renormalized current- 
current correlator at zero momentum transfer. 
Without vertex corrections, H a {9,uS) is given by 



n„ 



Q 



E 



d 2 k 



■G(k,iuj n + iuj m ) G(k,iuj m ) 
(102) 



The normalization factor Q has to be chosen such that 
in an ideal Fermi gas II (us ) = 1. The vertex corrections 
to H a in our model are due to irreducible corrections 
to the current vertex. These corrections obviously scale 
with dH/ckk PQ] and can be safely neglected. There are 
also corrections to Ii a due to extra interactions which 



are neglected in the spin-fermion model. Of particular 
relevance here is the residual p— wave type four-fermion 
interaction. This interaction give rise to RPA-type series 
of particle- hole bubbles p9,p3] . How strong these correc- 
tions are is a'priori unclear and has to be determined by 
experimental comparisons. We discuss this point in some 
detail in Sec [x|). 

For a conventional dirty metal with S(k, u>) = 
isignuj/2T, Eq. (102) yields the Drude formula a(u>) — 
(cjp;/47r)r/(l — iuir). In our case, however, the fermionic 
self-energy strongly depends on frequency, and the be- 
havior of <i{ui) is more complex. 

For experimental comparisons, it is convenient to ex- 
press the conductivity in terms of the generalized Drude 
model (lO | 



cr(w) = 



V 



t(u,T) 



A-k m*(uj, T)/m — iujT(uj, T) 



(103) 



where the effective relaxation rate t -1 ^) and the effec- 
tive electron mass m*(ui,T) depend on temperature and 
frequency. Both r and m* can be expressed in terms of 
<j\ and 02 as 
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Sometimes, the effective relaxation rate is also defined 



as 
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We now obtain the explicit expression for a(ui,T). 
Subs tituting the fermionic Green's functions into 
( |l02| ), using the spectral representation G(k,iuj m ) — 
(l/V) J dx ImG R {k. 1 x)/(x ~ iuj m ) to convert to real fre- 
quencies and retarded variables, and integrating over ek, 
we obtain after a simple algebra |93,M 



n«(0,n)= du> 



uf(uj + O) — 71f(w) 
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where S w = E(0,w). 

The relaxation rate is related to 11^ by 



V) = 



wlm(n£(0,w)) 



(108) 



(Im(n«(0, W ))) 2 + (Re(n«(0,o,))) 2 

= J y~ . . , means averaging over the Fermi 



where (. . .) 
surface. 

The angular dependence of the fermionic self-energy 
is due to the fact that both the coupling A and spin re- 
laxation frequency u) s f depend on 9 (see Eqn. (47); for 
small 9, k introduced in (|J7|) is related to 9 as k = kp9). 
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One can easily make sure that the integral over 9 is de- 
termined by 9 = O(l), i.e., the conductivities at low ui 
and T predominantly come from the momenta away from 
hot spots. Strictly speaking, this implies that very close 
to criticality, our forms of the self-energy obtained by 
expanding near hot spots, are inapplicable for the com- 
putations of the conductivity at small u> and T. Still, 
however, even at £ = oo, our formulas for the self-energy 
can be used to compute conductivity at frequencies above 
u s f(6 m ax), as at these frequencies, the whole Fermi sur- 
face behaves as a hot spot. At small frequencies and 
£ — ► oo, the critical behavior is affected by the instability 
in the pairing channel (see Sec. [X below), and has to be 
modified anyway. 

Below we focus on the forms of conductivities at finite 
£ when pairing fluctuations can be neglected. As we al- 
ready said and will discuss in more detail in Sec [X], at 
these £, the self-energy dependence on 9 is modest and 
u> s f{6max) is comparable to u) s f{9 = 0). In this situation, 
the 9 dependence of the self-energy affects the numbers 
but not the functional forms of the conductivities. In 
the calculations below, we just neglect the 9 dependence 
of £ without further discussions. We now consider the 
behavior of conductivities at various frequencies and tem- 
peratures. 



A. T = 

At T = 0, there are three different frequency regimes: 
i) u> <C w s /, ii) u) > Q and iii) u> 3> u s f. In the first 
regime, the fermionic self-energy has a Fermi-liquid form, 
Yj R {u>) = Xuj + i\oj 2 /4iu s f, and conductivities indeed also 
behave as in a Fermi liquid 
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FIG. 23. Real and imaginary parts of conductivity, effec- 
tive relaxation rate, and effective mass as functions of u in 
the Fermi-liquid regime. For defmiteness, we used A = 1. 

In Fig. g3| we plot the two conductivities, the effective 
relaxation rate and m* /m at ui < u> s f, obtained by solv- 
ing (107) with some angular- independent, "average" A 
and uj s f . The Fermi-liquid forms are clearly visible at the 
lowest frequencies. We see, however, that the deviations 
from the Fermi-liquid behavior become visible already at 
u> ~ L) s f/2. This is a consequence of the fact that the 
fermionic self-energy rapidly deviates from Fermi-liquid 
form at ui > w s //2. 




FIG. 24. Frequency dependence of a\, <72, 1/t, and m* /m 
at T = at intermediate frequencies for A = 1. 

At intermediate u> 3> u> > ^ s f, the fermionic S(w) can 
be approximated by a linear function of frequency up to 
~ 6 — 8uj s f. This results in the linear in to behavior of 



:« 



1/r and \jo\ In Fig. |i|we show the results for the con- 
ductivities, 1/r and m* jni in this region. The linearities 
are clearly visible. 
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FIG. 25. Frequency dependences of ax, 02, 1/t, and m*/m 
at T = as functions of lo/lu at w> u) s f- For simplicity, we 
set uisf = 0. 



Finally, at w » uj sf , E fl (w) = x/lu/2\lu\(lu + i\u\) be- 
comes independent of 9 i.e., the whole Fermi surface acts 
as a hot spot. In this limit, we obtained the exact ex- 
pression for H a {uj). 
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where a = y/(ito — S)/ui. This expression can be simpli- 
fied at w <C w and lu S> w (i.e., at a < 1 and a ^> 1, 
respectably). At w < w the bare Vt term in the de- 



nominator of (107) is small compared to the fcrmionic 
self-energy, while at lu ^$> lu, the self-energy is smaller 
than f2. At u) <C lu, i.e., n«l, we have 
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where A = 2- V2Tog (y/2 + 1) 
limit, (D<w, i.e., a 3> 1, we have 
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The behavior of the conductivities, the relaxation rate 
and the effective mass at lu 3> w s j and arbitrary lu/lu is 
shown in Fig]23. To emphasize the dependence on u>/u), 
we plotted the results right at ui s f = 0. Two features 
in Fig.|25| worth mentioning. First, l/r(w) scales as y/uj 
and virtually does not change between lu < lu and lu > lu 
(the two limiting results for 1/r are amazingly close to 
each other as to a good accuracy 1/A — 4/3!). Second, 
l/<7i(w) initially (at lu < ui) increases as \J~lu~, but then 
remains linear in to over a substantial frequency range 
and crosses over to w 3 ' 2 behavior only at very large lu/lu. 




FIG. 26. Frequency dependence of ImH a (u>) and 

ReU a (u>) at T = and cj s/ = 0. 

This behavior can also be extracted from the form of 
the polarization bubble presented in Fig. |26|: we see that 
Im~n. a {tjj) first increases as t/lu, but then nearly satu- 
rates at a value near 0.2 and very slowly decreases with 
increasing lu. It is amazing that ImW a (Lu) remains flat 
up to frequencies which by an order of magnitude exceed 
the lower boundary for y/u> behavior in the fermionic self- 
energy. 

We now assemble the T = results. For A 3> 1, 
lu = AX 2 Lu s f 3> uisf, and the scales u> and ui s f are well sep- 
arated in the sense that the upper boundary of the linear 
behavior of the fermionic self-energy (6 — 8w s /) is located 
below lu. In this situation, a^ (lu) is quadratic in lu at 
the lowest frequencies, becomes linear in watw~ u> s f/2, 
then crosses over to y/uj behavior at u> ~ 6 — 8ui s f, then 
again becomes linear at lu ~ ul, and eventually behaves as 



at the highest lu > lu. The last crossover is, however, 



w 3 / 2 .. .... ....... 

only suggestive as it occurs at very high lu ~ 10 2 <I> when 
our low-energy theory is clearly inapplicable. The inverse 
(72 behaves as lu in the Fermi-liquid regime, changes slope 
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but remains linear at w > u> s f then crosses over to to 1 / 2 
behavior at 6 — 8w s /, and at higher frequencies behaves 
as ijj. In the "y/u?' regime, real and imaginary compo- 
nents of the conductivity are close to each other. The 
relaxation rate 1/t(u>) is quadratic in frequency at the 
smallest frequencies, becomes linear in w at w ~ uj s f/2, 
and crosses over to y/uj form at lu ~ 6 — 8w s /. This 
behavior is schematically shown in Fig. p7[ 
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much with A, the ratio of cr^/vi indeed depends on A. 
This dependence may be used as a guide for selecting A 
which best fit the data. 
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FIG. 28. Real and imaginary parts of the conductivity for 
various A. 



FIG. 27. A sketch of the sequence of crossovers for the two 
conductivities and the relaxation rate for A ^> 1. The shaded 
regions are the extended crossover regimes. For A ~ 1, the 
crossover region at lo > ui a f extends over u>. In this situation, 
the two regions of linear behavior of crT" 1 merge. We found 
numerically that the slope of crj~ 1 (w) does not change much 
between the two crossover regimes, i.e., c-f 1 ^) remains nearly 
linear in lo over a very wide frequency range. 

We emphasize that that the behavior of conductivities 
atw> u> s f is fully universal and depends only on the 
ratio co/uj but not on the coupling A (see Eq. (110)). In 
particular, at, e.g., u> ~ w/2, 
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For A ~ 1, the upper cutoff for the linear behavior of 
fcrmionic self-energy exceeds u). In this situation, there 
should be no intermediate a/c<7 regime. Numerically, we 
found that n(w) does not change much with frequency 
between cu s f and u) > Q, i.e., the inverse conductivities 
and the relaxation rate remain linear in u> with little 
change of the slope over a wide frequency range. 

For completeness, in Fig. GH we present the results for 
the two conductivities for various A. We see that while 
the functional forms of the conductivities do not change 



1. Sum rule for optical conductivity 

The optical conductivity satisfies so-called /-sum 
rule |)5| associated with the conservation of the number 
of carriers 
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This sum rule has a practical application as it is used to 
extract the plasma frequency uj p i from the experimental 
data. The issue we consider is up to which frequency 
one has to integrate to satisfy t he s um rule. We will 
demonstrate that the integral in (114) is poorly conver- 
gent, and one has to integrate up to an unexpectedly high 
frequency to get the right value of u> p i . 

The results of the numerical integration of the conduc- 
tivity cti(w) up to some frequency u are presented in Fig. 
|2<4 We see that the convergence is very slow. This is re- 
lated to the fact that over a very substantial frequency 
range, the imaginary part of the polarization bubble re- 
mains almost a constant (see Fig. |2^), and hence <7i(u>) 
scales as 1/uo in which case the integral logarithmically 
depends on the upper cutoff in frequency integration. 
Only at very high u> 3> Co, the conductivity crosses over 
to oj~ 3 < 2 behavior, and the integral over <j\ converges. 
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FIG. 29. The sum rule (114) as a function of the upper 
limit of integration. We used T — 0, A = 1. 

In cuprates, the value of the plasma frequency is nor- 
mally obtained by integrating the measured cr\(uS) up 
to about 2 - 2.5eV (HJJgf-fl . We will argued below 
that this scale roughly corresponds to about 10o) Fig |2^ 
then indicates that the value of the plasma frequency 
extracted from these measurements is somewhat smaller 
than the actual value of uj p i , the difference decreases with 
increasing A. 



B. finite T 

We next consider the behavior of the optical conduc- 
tivity at finite T. We begin with the limit u> — > 0. 



1. finite T, lo^O 

At finite T, JmS(w) becomes finite at u> = 0, and 
Yi a {ui) at the lowest frequencies can be obtained by just 
expanding in u> in ( [lOTf ). Performing the expansion, we 
obtain 



w. 



P l 



1 



^T) = -^- 



1 



(IlO 



4tt 8T J_ oq ImS^ cosh 2 (w/2T) 



r-\T) 
m* /in 



lO 



^aV 1 



4ir 
16T 



1 + dRe^/dto 



du) 



[ImS^] 2 cosh 2 (w/2T) 



(115) 



Obviously, the typical frequency in these integrals is of 
order of temperature. At small T <C ujsj, the self energy 
has the Fermi-liquid form (see Eq. (|7q)). The straight- 
forward calculations then yield 



t-Ht) 



4tt 6T 2 
6T 2 



( w «/> 



6 ((1 + ^K 2 /) / 6 \ 



(116) 



where £(z) is Riemann zeta function. 




FIG. 30. Temperature dependences for at, 02, 1/r, and 
m* /m at finite T < u) a f and u) — 0. We used A = 1. 

To obtain the conductivity at T> uj s f we need as an 
input the self energy at u > w s /. The dominant con- 
tribution to its imaginary part comes from static spin 
fluctuations. At T -C Nuj s f which we assume to hold, 
this static piece is independent of frequency and just re- 
duces to XttT. Substituting this result into ( |115| ), we 
immediately obtain 



ciCO 



u) 
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8Att 2 T 
■ _1 (T) «2AttT 
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We see that both crj" and r _1 are linear in T, and the 
slopes scale with the coupling constant. The computation 
of the mass ratio requires more care as it involves _ReS^. 
For intermediate ui s f <C U ) <C w, 3Re£ fl /<9w 3> 1. Sub- 
stituting -ReE^ into (|llq ) and evaluating the frequency 
integral we obtain 



m*/m sw 0.48 V^/T 



(118) 



In the opp osite limit a; 3> w, <9Re£ fl /<9aj <c 1. The 
integral in (115) is then trivially evaluated and yields 



m* /m rs 1 



(119) 
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FIG. 31. Temperature dependences 

for cri(T), p(T) = l/cri(T), 1/r, and m* /m at u = and 
T> Lj 3f . We used \ = 1. 

In Fig.[H] we present the results for the dc conductivity, 
resistivity, relaxation rate and the mass ratio at T > w s y . 
We see that although the asymptotic "high temperature" 
behavior of the conductivity and the relaxation rate, Eq. 



(117) is reached only at rather high temperatures, the 
linear in T behavior of a^[ and r _1 (and, accordingly, 
the linear in T behavior of the resistivity) begins already 
at a low T ~ u) s f and over a wide temperature range can 
be well approximated by AT — B where A is somewhat 



smaller than in (117), and B is a small positive number. 



2. finite T, finite ui 




FIG. 32. Frequency dependence of the conductivities, in- 
verse a i , mass ratio and the relaxation rates 1/r and 1/r* at 
different temperatures. We used A = 1.7. The results were 
obtained with the second-order self-energy (see Fig 



Finally, in Fig. [32] we present the results for the con- 
ductivities, the mass ratio, and the relaxation rates 1/r 
and 1/r* when both frequency and temperature are fi- 
nite. We see that for any T, the inverse conductivity <j\ 
is still linear in to over a substantial frequency range, and 
the slope is almost independent of T. Observe also that 
the conductivity curves for various T cross each other 
such that at small to, conductivity decreases with T, while 
at large u it increases with the temperature. The con- 
ductivities in Fig. B3 are obtained with the second-order 
fcrmionic self-energy (its imaginary part is plotted in Fig. 
pfjl). We verified that if we use the "full" self-energy (218), 



we obtain almost identical results for the conductivities, 
mass ratio and relaxation rates. 



C. summary of Sec. VII 



We now summarize the main results of this section 

1. At T = 0, real and imaginary parts of the opti- 
cal conductivity, <j\(uj) and 02(0;), and the effective 
scattering rate l/r(w) undergo a series of crossovers 
schematically shown in Fig J27J. In particular, for 
moderate couplings, crj~ and u^ are linear in lo 
over a substantial frequency range. The relaxation 
rate 1/r is also linear in 10 but over a smaller range 
of frequencies, and then crosses over to y/uj behav- 
ior. 

2. We found that the /—sum rule for u\ is recovered 
only if the integration in I(u>) = J? ai(x)dx is ex- 
tended up to extremely high frequencies. If the in- 
tegration is performed up to frequencies about 10 
times larger than u, I(w) saturates (see Fig. |2^) at 
some intermediate value which depends on A. The 
implication of this result for experiments is that 
ojpi extracted from integrating up to a particular 
cutoff is somewhat smaller than the actual plasma 
frequency. 

3. At finite T, the most important features are the 
development of the peak in 02 (w) at to comparable 
to T, and the crossing between <j\{lo) curves at dif- 
ferent T: the conductivity at the lowest T is the 
largest at the smallest u>, but become the smallest 
at larger u> (see Fig. [32(). The relaxation rate on 
the other hand do not display this kind of behavior. 

4. The effective mass m* is by a factor 1 + A larger 
than the bare electron mass, but the ratio m* /m 
decreases with increasing frequency and tempera- 
ture. 
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VIII. OTHER OBSERVABLES 

In this section, we briefly discuss the behavior of other 
observables - the fermionic density of states, the NMR 
relaxation rate and the Raman intensity. 



A. fermionic density of states 

The fermionic density of states (DOS) is given by 

N(uj) = -y\ImG(k,u;)\ (120) 

k 

Assume for simplicity that the Fermi surface is circular. 
Then we can split momentum integration into the inte- 
gration over d9 along the Fermi surface, and over de k — 
vpdk±, where fen is a momentum component along the 
Fermi surface. Substituting G _1 (k, lj) = w+S(fcii , lj) — e k 
and integrating over de k first we obtain that 



N(u) = const 



(121) 



independently of the strength of the fermionic self-energy 
It is only essential that the self-energy does not depend 
on e k . 

For more complex Fermi surface, N(lj) does acquire 
some frequency dependence which, however, is com- 
pletely unessential from physics perspective. Also, Eq. 



(121) is indeed only valid at frequencies much smaller 
than the fermionic bandwidth, when one can linearize 
the fermionic dispersion near the Fermi surface. At fre- 
quencies comparable to bandwidth, N(lj) curves down 
and eventually vanishes. 



B. neutron scattering and NMR relaxation rate 

The neutron scattering measures the dynamical struc- 
ture factor S (q, w) = (l + 2n.B(w))x // (q, lj). Substituting 
the expression for the susceptibility, Eqs (jlq), (po|), we 
obtain 



S(q,w) = 



w(l + 2n B (w)) 



x e 

w sf (1-K 2 (q-Q) 2 ) 2 + (^M/) 



(122) 



At lo ^> T, the Bose factor is irrelevant, 5(q, lj) at a given 
q has a simple x/(l + x 2 ) form where x — w/(w s /(l + 
£ 2 (q ~ Q) 2 )) is a rescaled frequency. 

The NMR relaxation rate 1/Xi is proportional to the 
momentum integral of S(q, lj) weighted with the momen- 
tum dependent formfactor F(q) p7[: 



Tf 1 ^) «£>(q)S(q,w) 



(123) 



The typical frequencies in the NMR experiments are in 
the range of few MHZ (~ W'^K), hence T > u> for all 
reasonable T. Then 1 + 2ns (lj) ~ 2T/u>, and 



7T 1 oc2Tlim Vf(q) 



*"(q,w) 



LO 



(124) 



The form factor F(q) depends on the local environment 
of the isotope atom. For Cu NMR, F(q) d oes n ot vanish 
at Q = (ir,n). Substituting x"(q, w) into (124), one c an 
easily check that the momentum integration in (124) is 
then confined to small |q — Q| ~ £ -1 , and hence 



1 



^(Q)xo 



Attlj 



sf 



(125) 



As the product F(Q)xo can be extracted independently 
from Knight shift data p7| , the measurements of 1/T± 
allow one to determine w s f which, we remind, is one of 
the two input parameters in our theory. 



C. Raman intensity 

Another observable sensitive to the form of the 
fermionic self-energy is the intensity of the Raman ab- 
sorption. In a Fermi gas, a transferred photon energy lj 
can be absorbed by a metal only if a; is smaller than vp\q\ 
when q is a transferred momentum (this is the Fermi 
golden rule) [ 100 1 . As photon momentum is vanishingly 
small, the absorption is possible only for extremely small 
lj ~ LJinVp/c where LJi n is the incident photon frequency, 
and c is the velocity of light. However, if fermions have 
a nonzero ImS, the absorption is possible for all frequen- 
cies: the energy extracted from photons is dissipated due 
to a finite lifetime of fermionic excitations. 

The intensity of the (non-resonance) Raman absorp- 
tion R(lo) is generally proportional to the imaginary 
part of the particle-hole bubble at zero external momen- 
tum and finite frequency, weighted with the form-factors 
which depend on the scattering geometry 1 101, 102) ] . For 
the most studied B\ g scattering, the form- factor Fr(1s) cx 



cos k x — cos k y is at maximum near hot spots, and can 
be approximated by a constant. The Raman intensity is 
then simply given by 



R(lj) oc ItiiTIr(lj) 



(126) 



The Raman bubble Hr(lj) differs from n CT (w) as side 
vertices are now scalars. As we discussed in Sec. IV, 



a scalar vertex with zero total momentum is subject to 
strong vertex corrections. For a density-density correla- 
tor with momentum independent vertices, vertex correc- 
tions fully cancel fermio nic self-energy, and as a result 
there is no dissipation [ 103 1 . [This cancellation is the 
consequence of the fact that the number of particles is 
conserved]. However, for B\ g Raman scattering, the side 
vertex has d— wave symmetry, and the renormalization of 
the vertex requires a d— wave component of the interac- 
tion. In this situation, there is no exact relation between 



:-!,x 



fermionic self-energy and vertex corrections. We explic- 
itly verified that for magnetically-mediated scattering, 
the correction to the B\ g Raman vertex is nonsingular 
even when £ — > oo. Neglecting this nonsingular vertex 
renormalization, we obtain a simple relation between the 
normal state Raman intensity and the normal state con- 
ductivity: 



-R(w) ~ u/ai(u)). 



(127) 



In particular, at T = 0, i?(w) oc w at small frequencies, 
then saturates, then crosses over to R(u>) ex a/c<7, and then 
(if lattice effects do not interfere) first saturates and then 
begins decreasing as uj^ 1 ' 2 . This last regime, however, 
is likely already masked by the development of a "two- 
magnon" peak from short-range spin fluctuations [104]. 



IX. PAIRING VERTEX 
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FIG. 33. The diagram for the full vertex in the parti- 
cle-particle channel. The bare (unshaded) vertex is infinitesi- 
mally small. The wavy line is the spin propagator. Since the 
spin-mediated interaction is retarded, the full vertex depends 
on frequency uj. 

Finally, we briefly discuss the development of the pair- 
ing instability in the spin-fermion model. It is cus- 
tomary for the analysis of the pairing problem to in- 
troduce an infinitesimally small particle-particle vertex 
gp p (k,—k) = g pp (k) and study its renormalization by 
magnetically mediated interaction. The corresponding 
diagram is presented in Fig |33|. The temperature where 
the renormalized vertex diverges is obviously the onset of 
pairing. The magnetically-mediated pairing interaction 
is repulsive for s— wave pairing, but is attractive in the 



d x 2_ y 2 channel [|17j,|25|. Accordingly, we will be searching 
for a d— wave solution for the full pairing vertex g pp (k). 
For momenta near hot spots, the d— wave constraint im- 
plies that g PP (k + Q) = —g PP {k). 

We demonstrate in this section that the renormaliza- 
tion of the pairing vertex is not a 1/A effect, and hence 
the pairing instability is possible even in the limit when 
A — ► oo, when our Eliashberg-type theory becomes ex- 
act. 

For a particle-hole vertex with a nonzero momentum 
transfer, the 1/A smallness of the vertex renormalization 
comes from the fact that for large N, the momentum de- 
pendence of the spin susceptibility can be neglected, and 
the vertex correction is a convolution of the two fermionic 



densities of states and the dynamical spin susceptibility 
at momentum Q. The latter contains N in the denomi- 
nator, and this gives rise to the overall 1/N factor in the 
correction to a particle-hole vertex. For zero momentum 
transfer, the 1/N does not appear, but as we discussed 
in Sec IV D , this is the consequence of the fact that the 



momentum integral is singular and has to be properly 
regularized. 

For a particle-particle vertex, the situation is different 
as now vertex renormalization involves fermions with op- 
posite momenta and frequencies, and hence at a given 
frequency, the poles in the two fermionic propagators are 
located in different half-plains of e^. As a result, there is 
no double pole in the momentum integral. The latter in 
turn implies that at N — ► oo, one can again factorizing 
the momentum integration i.e., neglect in the spin sus- 
ceptibility the momentum component transverse to the 
Fermi surface. Still, however, one has to integrate the 
susceptibility over the momentum component along the 
Fermi surface. This last procedure changes the power 
of 1/N in the denominator such that eventually 1/N is 
eliminated from the equation for g pp {k) (see below). 

Two of us and Finkel'stein have demonstrated explic- 
itly [|40| that for relevant k near each of hot spots, the 
momentum variation of g pp (k) is irrelevant and can be 
neglected. The momentum integration then can be per- 
formed explicitly, and at T — we obtain 



5g p 



oc 



3 5 2 



9 PP Sn 2 v F 



dWr, 



\u m - iS(w m )| 



X L{u m -n m ) (128) 



where fl is the external frequency (g pp = g pp (£l, — J7)), 
and Xl(v) = 7TXo£/(i + M/t^/) 1 / 2 is the "local" ID 
susceptibility obtained by integrating over momentum 
component along the Fermi surface and neglecting q — Q 
transverse to th e Fe rmi surface . Substituting this sus- 
ceptibility into (12S), we obtain 
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1 + (2A) + y/l + \uj m \/usf 



(129) 



We see that the vertex correction is positive, i.e., the 
addition of a spin-fluctuation exchange enhances the 
particle-particle vertex. Rescaling frequency as x = 
w/u) a f we immediately find that 5g PP / g pp depends only 
on A— independent effective coupling A. The origin of 
this result is rather transparent in the quantum-critical 
regime uj m ^> w s /. Here xl oc ^/u> s f/u> m scales as 
(1/A) 1 / 2 . Simultaneously, 



oJ„ 



iE(u m ) ~ y/uj^uJ ~ \uj m . / 



(130) 



also scales as (1/A) 1 / 2 . Taking the ratio of the two ex- 
pressions we obtain A— independent result 
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Another issue is how large is the pairing instability 
temperature. For A < 1, the answer is rather obvious as 
the pairing problem is very similar to the BCS phonon 
problem with u> s f playing the same role as the Deby e fre - 
quency for phonons. Indeed, when A < 1, u> s f in (12£) 
provides a natural u pper cutoff for the frequency inte- 
gral. Extending Eq. (129) to finite T in a standard way 
and evaluating the frequency sum with logarithmical ac- 
curacy, we obtain 



s fJ, 



pp 



tfpp 



A uj sf 

TTA log ^ 



(132) 



Obviously then, Ti ns oc w s /exp ( 1+A )/ A . This is similar 
to the McMillan formula for phonons [ 105 1 . The 1 + A 



factor in the numerator accounts for the mass renormal- 
ization at frequencies below uj s f. 

Eq. (132) is formally valid for all couplings A provided 
that the pairing problem is confined to a region u) < Lu s f 
where the system has a Fermi-liquid behavior. However, 
we see from (129) that the pairing kernel preserves 1/oj 
form also at frequencies larger than u/ s f and crosses over 
to 1/ui 3 ' 2 f orm only at lo > Q. For u> ^> u m ^> u) s f, we 
have from (p9|) 
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Alternatively speaking, the kernel for the pairing prob- 
lem remains 0(1/uj) not only in the Fermi-liquid regime, 
but also in the quantum-critical regime where the system 
does not behave as a Fermi-liquid. This implies that Ti ns 
in principle can be of order to, i.e., substantially higher 
than in the McMillan formula. This is not guaranteed, 
however, as the pairing kernel in ( |129| ) depends on both 
internal and external frequency, and the logarithm in the 
r.h.s. of ( |133| ) is cut by fl m . In the latter case, the lad- 
der series are not geometrical, and there is no a'priori 
guarantee that they yield an instability at T ~ Q. 

The full analysis of the pairing problem requires extra 
care and is beyond the scope of the present paper. We 
just cite the result |40(| : the instability temperature Ti ns 
does scale with Co at strong coupling, and saturates at 
Tins ~ 0.1 7Q at £ = oo. Num erically, the saturation of 
Ti ns was detected in Ref. |l06| . The behavior below Ti ns 
is rather involved, as we discussed in the Introduction 
(see Fig. §). 

Note in passing that the studies of the pairing insta- 
bility at the verge of the magnetic transition yield rather 
nontrivial results also when the pairing is mediated by 
ferromagnetic spin fluctuations [107] . 

The fact that the pairing instability temperature scales 
with u> which, we remind, is the upper cutoff frequency 



for the quantum-critical behavior, implies that our nor- 
mal state analysis which neglects pairing channel is valid 
only in a restricted frequency and temperature range. 
Still, however, numerically Ti ns <C u, and hence there 
exists a wide region of frequencies/temperatures where 
normal state quantum-critical analysis is valid. A related 
issue is how strong are the pairing fluctuations above 
Tins- We show below that the pairing susceptibility has 
an overall factor 1/iV and hence pairing fluctuations only 
affect fermionic self-energy in a narrow range near Tj„ s . 
Outside this region, they contribute only 1/N corrections 
to the fermionic self-energy and hence do not affect our 
N — oo analysis. 




FIG. 34. The diagram for the correction to the electronic 
self-energy due to fluctuations in the pairing channel. 

The corrections to the electronic self-energy due to 
pairing fluctuations are shown diagrammatically in Fig 
|3J. Consider first the shaded regions in Fig. f54| Each 
time we add an extra bosonic line to the ladder, we also 
add two fermionic lines. As we demonstrated above, this 
addition does not contribute extra powers of 1/N. The 
1/N factor comes from the unshaded part of this dia- 
gram. It contains a pair of fermionic propagators, and 
a pair of bosonic propagators, i.e., it contains an extra 
bosonic propagator compared to the building block in 
the shaded part of the diagram. The bosonic propagator 
scales as 1/N as there are N channels for spin decay. We 
explicitly verified that this 1/N appears as the overall 
factor in the diagram, i.e., 



Ep P (fi) oc 



N 



(134) 



where T is the enhancement factor from the shaded part 
of the diagram. The divergence in T implies pairing insta- 
bility and obviously affects fermionic self-energy How- 
ever as long as T < N, S pp (0) has a prefactor of 1/N, 
and can be neglected in the limit N — > oo. 

A final comment. The pairing instability at a finite 
Tins at ^ = oo obviously implies that there exists a dome 
on top of the magnetic quantum critical point. Inside this 
dome, the quantum-critical theory has to be reconsidered 
as the fermions near hot spots become gaped. In other 
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words, the quantum-critical behavior which causes the 
pairing instability at a finite Ti ns by itself gets affected 
as a feedback from the pairing. In particular, we explic- 
itly verified that the correction to the spin-fermion vertex 
is saturated at finite value below Ti ns , i.e., it is no longer 
logarithmically singular. This, in particular, implies that 
the anomalous magnetic response in the quantum-critical 
regime, which we found in Sec. M exists only at interme- 
diate energies and cannot be extended to the lowest w 
(which fall inside the dome). The quantum-critical be- 
havior inside the dome is a subject of a separate study 
and we will not discuss it here. 



At strong coupling, the region near the QCP is divided 
into the Fermi liquid regime and the quantum-critical 
regime where the system behavior is the same as at crit- 
icality. The quantum-critical behavior extends roughly 
up to frequencies comparable to the spin-fermion cou- 
pling constant g. The crossover from the Fermi liquid to 
the quantum-critical behavior on the other hand occurs 
at energies of order w s f ~ ^ 2 9 <C g. We found that 
there is a single crossover energy for both electronic and 
magnetic properties of the system. 

We now list the catalog of our key results. We first 
list the results at T = and then show how they are 
modified at finite T. 



X. SUMMARY AND EXPERIMENTAL 
COMPARISONS 



A. T = 



In this paper, we presented the full scale calcula- 
tions for the normal state properties of the spin-fermion 
model. The model involves low-energy fermions inter- 
acting with their own spin collective excitations which 
are assumed to be peaked at the antiferromagnetic mo- 
mentum Q = (71", tt). We argued that this model is the 
low-energy version of the lattice, Hubbard-type models 
for strongly interacting fermions, provided that the spin 
fluctuations are the only low-energy collective bosonic 
excitations. The model indeed makes sense only if the 
spin-fermion coupling does not exceed fermionic band- 
width which we assume to hold. If this condition is not 
satisfied, the separation between low-energy and high- 
energy excitations becomes problematic. 

We demonstrated that near magnetic instability, the 
model falls into the strong coupling limit, where the con- 
ventional perturbation theory does not work. The cor- 
rections to the Fermi-gas behavior are the strongest for 
fermions near hot spots - the points at the Fermi surface 
separated by antiferromagnetic Q. The presence of hot 
spots is essential for our study - without them the critical 
theory at T — would have had a different dynamical 
exponent z = 1. 

We developed a controlled way to perform calculations 
in the strong coupling limit by expanding in 1/N where 
N(= 8 for the physical case) is the number of hot spots 
in the Brillouin zone. Equivalently, one can extend the 
model to a large number of electron flavors M and expand 
in 1/M. 

Our major results are summarized in Figs, nil and 
p0| for the fermionic spectral function, and Fig. |32| for 
optical conductivity. For the dynamical structure fac- 
tor, our major result is Eqn (p7[). We demonstrated that 
near the QCP, the dimcnsionless spin-fermion coupling 
A ex £ is large, i.e., the system falls into the strong cou- 
pling regime. This strong coupling behavior can also be 
reached at intermediate £ when the spin-fermion interac- 
tion g increases. 



1. Fermi-liquid regime 

In the Fermi-liquid regime, we found for the fermionic 
self-energy 



£(ft,fi) = flX(k)+iX{k) 



n|o| 



(135) 



The momentum-dependent coupling constant X(k) and 
and u) s f(k) are given by 

A(fc) = A/(l + (kf) 1/2 , "./(*) = w./(l + (kf) 

(136) 

where k is the deviation from a hot spot along the Fermi 
surface (k = e^q/vf). These Fermi liquid forms are 
valid when Q < u a f(k). Obviously, the self-energy cor- 
rections are the largest for fermions near hot spots. For 
the same frequencies, we found that the spin susceptibil- 
ity near the antiferromagnetic momentum Q can be well 
approximated by its static form 



x(q,w) 



x e 



l + (q-Q) 2 £ 2 



(137) 



In this limit, <J\{u>) weakly depends on u>, and 02 (w) 
scales as w _1 . 



2. quantum- critical regime, intermediate frequencies 

At intermediate frequencies co s f(k) < Q < 6 — 8u> s f, we 
found that to a surprisingly good accuracy, the imaginary 
part of the fermionic self-energy is linear in frequency: 



S"(fc, n) « 0.3A(fc)(n - 0.7w s/ ) 



(138) 



Note, however, that this behavior is not an intermedi- 
ate asymptotic which width could be controlled by some 
parameter of a system, but rather a peculiar behavior 
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in a wide crossover region between truly Fermi-liquid 
and truly quantum-critical (high frequency) regimes. 
The real part of the self-energy E'(fc, 51) roughly follows 
lj log M behavior which is the Kramers-Kronig transform 



of (138), but this form is only approximate as the linear- 
ity in S" exists only in a limited frequency range. 

At frequencies larger than cj s f, one also has to use the 
full expression for the dynamical spin susceptibility as 
neither the gap in the static susceptibility nor the dy- 
namical piece can be neglected 



x(q,w) 



xoe 



l + (q-Q) 2 e-i\to\/u> sf 



(139) 



Finally at these frequencies, both <7i(lo) and 02 (u>) scale 
as 1/uj. 



3. quantum- critical regime, high frequencies 

At high frequencies, 51 ~ Co 3> U) s f{k), the whole Fermi 
surface behaves as a hot spot. We found that in this 
regime, the fermionic self-energy scales as y/uj: 



£(ft) = e^V^P|signn 



(140) 



This form of the self-energy implies that fermionic 
Green's function G~ 1 (k,Q) ex (zQ — ef./Cj). Compar- 
ing this with the Fermi liquid result G~ 1 (fc, 51) ex 51 — 
e fc/(l + <Kfc))j we see that as the system crosses over to 
the quantum-critical regime, the pole in the fermionic 
propagator gradually moves from the real frequency axis 
onto the imaginary axis. The pole along imaginary axis 
implies that fermions are completely overdamped and can 
only propagate diffusively. 

Simultaneously, the dynamic spin susceptibility ac- 
quires a strong frequency dependence which also gives 
rise to the diffusive behavior near q — Q: 



x(q,w) 



\0 



(q-Q) 



*M7 



(141) 



where 7 = (ujsf^ 2 )^ 1 is independent of £. Alternatively 
speaking, both fermionic and spin excitations are diffu- 
sive at w ~ u > uj s f. 

Finally, in the true quantum-critical regime 01 (w) in- 
terpolates between a; -1 / 2 at u> < u>, and w~ 3 / 2 at the 
highest frequencies. The imaginary part of conductiv- 
ity interpolates between ui^ 1 ' 2 for lo < Co and lo^ 1 for 
oj > Co. We found that the crossover region for 02 (w) is 
rather narrow, but that for (Ji(lo) is very wide, and in a 
wide range of frequencies lo > Co, Air/ to 2 1 01(10) is close to 
Q.2UO- 1 . 



4- quantum- critical point 

Finally, we analyzed the forms of the self-energy and 
the spin susceptibility when u> s f — > 0, and for fermions 



near hot spots the quantum-critical behavior extends 
down to the lowest frequencies. In this limit, there ex- 
ists a wide frequency range where on one hand lo 3> ui s f , 
and on the other hand to <C Co, i.e., a truly low- energy 
quantum-critical regime. We found that in this regime, 



Eqs (14C) and (141) are modified by subleading, \og(co/Co) 
terms in the perturbation series. We studied the effects 
of extra logarithms in the one-loop RG theory assuming 
that the number of hot spots N is large and neglecting 
nonlogarithmic corrections in 1/N (N — 8 in real situa- 
tion) . We found that the spin dynamics is still described 
by the dynamical exponent z = 2, but the dynamical 
susceptibility acquires an anomalous dimension: 

X (q,^„)cx(7k m | + (q-Q) 2 )- 1+ '' 

[log( 7 K„| + ( q -Q) 2 )]-^/ 3 ^ (142) 

where 77 = 2/N = 0.25, and and v y and v x are the 
components of the bare Fermi velocity near hot spots: 
£fc = v x k x + v y ky]e k +Q = ~v x k x + v y k y . The appear- 
ance of a finite n is the consequence of the fact that the 
spin-fermion model contains anomalies which give rise 
to a complex structure of the four-boson vertex at the 
QCP. In this respect, our theory differs from the effec- 
tive bosonic <p A theory, which is marginal in d — 2 since 
z + d = 4. 

The fermionic self-energy is not affected by the anoma- 
lous dimension r\ and differs from (140) only by a loga- 
rithmical factor: 



E(51 m ) oc |51 ?) 



1 1/2 



log 51 r , 



-1/2 



(143) 



Finally, the Fermi velocity is also logarithmically renor- 



malized in the quantum-critical region 



i£ = v x I 1 



2ALv y V /2 ,. _ 
^~ x ;Vy ~' !l 



2AL v y 
nN v x 



-1/2 



(144) 



where L = | log[min(£ -1 ,7a>)]| This renormalization im- 
plies that as £ diverges, v x — > 0, and the velocities at khs 
and k/j S + Q become antiparallel to each other, i.e., the 
Fermi surface becomes "nested" at hot spots. This "nest- 
ing" is the first step in the transformation from a large 
Fermi surface to a small one, consisting of hole pockets. 
If this effect occurred at a finite £, then one might ex- 
pect a subsequent topological transition in which a large 
Fermi surface disconnects into hole pockets and the rest. 
We found, however, that the "nesting" occurs only right 
at £ = 00, and is only a weak logarithmical effect. This 
is a consequence of the fact that the low frequency spin 
excitations are diffusive rather than propagating. If the 
Fermi surface did not contain hot spots or the damping 
was just weaker, one could expect more stronger nesting 
effects and, accordingly, some spin-density-wave precur- 
sors in the normal state. 
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The logarithmical modification of the fermionic self- 
energy gives rise to logarithmical modifications of the 
conductivities: both real and imaginary parts of the con- 
ductivity behave as a(u)) cx | logwl 1 / 2 /-^. 

We next list the results at finite T. 



B. finite T 

1. Fermi liquid regime 

The Fermi liquid regime for T dependent terms in the 
fermionic self-energy is confined to nT < 0.5to s f(k). For 
these temperatures we obtained E(fc,f2) = E(fc,f2, T = 
0) + £ T (fc,ft) where £(fc,fi,T = 0) is given by dl35|), 



and £r(fc,£7) predominantly affects the imaginary part 
of the self-energy: 



s T (fc,n) 



2t^2 



n 2 T 



Au af {k) 



signal F 



n 



u> 8 f{k) 



0(T 4 ) 

(145) 



where F(x) is a smooth function of x with the limits 
F(0) = 1, F(x > 1) -> 1/3. The full expression, includ- 
ing T 4 terms, is presented in ( [213J ). 

We found that the spin susceptibility is not affected 
by a finite temperature except for a regular, (9((T/A) 2 ) 
correction to the correlation length. The conductivities 
are affected and change to <ji(uj) cx T~ 2 and <T2(u>) oc 
uj/T 4 at u> < T 2 /uj s f. The resistivity indeed obeys a 
Fermi liquid form p(T) oc T 2 . 



2. quantum- critical regime, intermediate temperatures 

At intermediate temperatures, 0.5w s / < T < Nuj s f, we 
found that the thermal contribution to the self-energy is 
predominantly linear in T In the formal limit N 3> 1, 
there are two different types of linear behavior. One oc- 
curs at relatively high T, when on one hand, T < Nui s f, 
and on the other hand, T is larger than 6 — 8ui s f- In this 
regime, the dominant piece in E-r(k, 0) is the impurity- 
like contribution from the scattering by static spin fluc- 
tuations. It yields Ey(k, 0) = iixTX^. The contribution 
to Ey (k, 0) from dynamical spin fluctuations is smaller 
and scales as \/T. In explicit form (see Eqn fl213| )) 



££(fc,0) = nT\(k) - 1.516 ( ^ 



1/2 



(146) 



At T < 6 — 8uj s f , the scattering from dynamical spin 
fluctuations also yields a linear in T piece in the fermionic 
self-energy. In this temperature range, we found, to a 
good accuracy 



At finite tt the self-energy is more involved and is gen- 
erally given by Eq. ( |218| ). The dynamical piece in the 
self-energy is well approximated by the N — oo result, 
Eqns. ( pfj| ) and (|9l|). In the static piece, the corrections 
to the N = oo result iirTX become relevant above typ- 
ical frequency evolves is f2 = N 2 ui s f(vp/v x v y )/36. At 
N = oo and finite A, this fj exceeds Q = 4X 2 cj s f. How- 
ever, at finite N and large A, O C w. We found a con- 
trollable way to compute the full fermionic E(f2) at T <C 
Nui s f and arbitrary f2/f2. We found that at uj > il > O, 
the static piece in the fermionic self-energy decreases as 
E'jt(fi) oc TAr v / ^7niog(fi/(A^ 2 w ;s/ )). At these frequen- 
cies, the full self-energy (sum of dynamical and static 
pieces) obeys the scaling form E(Q,T) oc y/Tf(fl/T), 
up to logarithmic corrections. For completeness, we pre- 
sented both the full result for the self-energy, Fig. |2^, 
and the second-order result, Fig g(]. We argued that al- 
though the derivation of the full self-energy is justified 
for T < Nto s f, by numerical reasons, vertex corrections 
may become relevant by numerical reasons already at 
T > uj s f. We conjectured that at these temperatures, 
the second-order result may be closer to reality than the 
"full" expression which includes self-energy of intermedi- 
ate fermions but neglects vertex corrections. 

The dynamical spin susceptibility is the same as in 
( |139| ), and the conductivities scale, at vanishing u>, as 
a i oc 1/T, <72 oc uj/T 2 . The behavior of conductivities 
at finite frequencies is rather involved and we refrain 
from discussing the limits. The full result is presented 
in Fig. |32|. 



3. quantum- critical regime, high temperatures 

At high temperatures, nT > Nuj s f, vertex corrections 
and the corrections to the Fermi velocity cannot be ne- 
glected, and our controlled computational scheme breaks 
down. In this regime we obtained the results by neglect- 
ing vertex corrections without justification and perform- 
ing self-consistent FLEX -type calculations. We then ver- 
ified that the vertex correction, evaluated using the full 
fermionic propagators remain 0(1), i.e., at least they do 
not diverge. We found that the contribution to Ey(k, 0) 
from scattering on static spin fluctuations scales down 
from ttTX and at the highest T behaves, up to logarithms, 
as vT with £— independent coefficient, i.e., in the same 
way as the contribution from the scattering on dynamical 
spin fluctuations. Combining the two contributions, we 
obtained 



E'^(fc) 



- 8 (£) 



1/2 



log 



ST 



Nuj 



1.516 



«/ 



1/2 



(148) 



',(fc,0) = 0.4 7rTA(fc) 



(147) 



This evolution of the thermal piece in the self-energy 
affects the conductivities which scale at ui — ► as 
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a 1 (T)(xT- 1 / 2 ,(j 2 (T)(xuj/T. 



4- quantum- critical point 

At u> s f — > 0, the temperature range irT > Nuj s f ex- 
tends down to the lowest frequencies. We found at van- 
ishing to s f and a finite T that the fermionic self-energy 
scales as 



S^ocT^llogTl 



(149) 



The angle between the F ermi velocities at k^ s and k/ ls +Q 
scales as 1/| logT| (see ( |l00| )). Both of these results are 
identical to what we o btaine d at T = 0, if we substitute 
T by SI (see Eqs. (|l43| ) and ( |l44[ )). Equal powers of log- 
arithms at T = 0, fl ^ and T ^ 0, Q = mean that the 
system at the QCP possesses uo/T scaling for fermionic 
variables. This is quite expected, in view of the fact that 
the fermionic self-energy is not affected by the anoma- 
lous exponent in the spin susceptibility. However, as we 
said, the self-consistent FLEX solution is uncontrolled, 
and therefore the equivalence between $1 and T in our 
solution is not the proof that the oj/T scaling actually 
exists. 

On the other hand, our theory definitely yields no lu/T 
scaling for the dynamical spin susceptibility. At T = 0, 
we found the anomalous exponent r\ (see ( |142[ )). At finite 
T and u = (q — Q) 2 = 0, we found that the susceptibil- 
ity behaves regularly, as x oc 1/T 2 . This behavior is 
certainly modified due to (f> 4 interaction between bosonic 
modes. This effect is, however, beyond the scope of our 
theory as the </> 4 interaction is produced by high-energy 
fermions and scales as 1/A. The anomalous behavior at 
T = on the other hand is universal and survives in the 
limit A — > oo. 



C. comparison with experiments 

We now compare our key results for the spectral func- 
tion and conductivity with the experimental data on op- 
timally doped high T c cuprate Bi22l2 for which both 
photoemission and conductivity data are available. We 
remind that the two inputs in our theory are the coupling 
constant A and the overall scale Q — 4\ 2 uj s f. Alterna- 
tively, we can reexpress A as A = 3u.f£ /(16w a /) and 
use i>f£ -1 and u> s f as input parameters. 

The value of the Fermi velocity can be obtained from 
the photoemission data at high frequencies, when the 
self-energy corrections to the fermionic dispersion be- 
come relatively minor. The three groups which reported 
the MDC data on Bi2212 for momenta along zone di- 
agonal PJ73,p4] all agree that that the bare value of 
the Fermi velocity along the diagonal is rather high: 
2.5-3eFA or 0.7-0.8ey a where a ~ 3.8A is the Cu-Cu 



distance. We used the t — if model to relate this velocity 
with that at hot spots. Using the value of the velocity and 
the experimental facts that the Fermi surface is located 
at k sa (OAir/a, OAir/a) for momenta along zone diagonal 
and at k sa (w/a, 0.2ir/a) for k x along the zone boundary, 
we found t ~ 0.2 - 0.25eV, if w -0.36* and ft sa -l.lt. 
These numbers roughly agree with other studies [ p3[ . 
The hot spots are located at k/ ls = (0.167T, 0.847r) and 
symmetry related points, and the velocity at a hot spot 
is approximately a half of that at zone diagonal. This 
yields v F « 0.35 - OAeVa. 

The values of w s f and £ can in principle be extracted 
from neutron scattering data on 5(q, lj) oc w/((1 + (q — 
Q) 2 £ 2 ) 2 + (u/u sf ) 2 ), and from NMR data. We are not 
familiar with the detailed analysis of the normal state 
neutron and NMR data for Bi2212 and will rely on the 
data for near optimally doped YBCO. The NMR anal- 
ysis |7(| yields w s / ~ 20meV and £ ~ 2a. The neutron 
data in the normal state are more difficult to analyze 
because of the background which increases the width of 
the neutron peak and masks some frequency dependence. 
The data show |75fl that the dynamical structure factor in 
the normal state is indeed peaked at q = Q = (it /a, n/a), 
and the width of the peak increases with frequency and 
at w = 50meV reaches 1.5 of its value at u> = 0. A 



straightforward fit to the theoretical formula, Eqn(139), 
yields w s / ~ 35 — 40meV and £ ~ a which are predictably 
larger than u) s f and £ extracted from NMR. We will use 
NMR values u) s f ~ 20meV and £ = 2a for further esti- 
mates. 

Combining the results for vf, £ and ui s f, we obtain 
A ~ 2. This in turn yields w ~ 0.3eV. As an in- 
dependent check of the internal consistency of the es- 
timates, note that our recent analysis of the supercon- 
ducting state |39|1 yields the resonance neutron frequency 
at ui res ~ 0.25<D/A, i.e., at uj res ~ iOmeV. This is quite 
consistent with the data. 

Away from hot spots, the effective coupling decreases 
as X(k) = A(l + {ek+Q/vF^- 1 ) 2 ) 1 ' 2 . This formula is in- 
deed valid only in some vicinity near hot spots as we 
didn't include in the theory the variation of the Fermi ve- 
locity along the Fermi surface. Nevertheless, experimen- 
tally, even for k = kp along zone diagonal (the furthest k 
point away from a hot spot), ekp+Q/vFikp)^ 1 ~ 1.4 [pj. 
This shows that the coupling does vary along the Fermi 
surface, but this variation is modest: A is reduced by at 
most 1.7 when we move from hot spots towards zone diag- 
onal. Actual reduction can be even smaller as vf by itself 
increases as one approaches zone diagonal. For u) s f(k), 
our theory predicts that it increases with the deviation 
from hot spots. Note, however, that uj s f oc sin^o, where 
<j>o is the angle between velocities at k and k + Q. In 
our theory, we assumed that this angle does not change. 
In reality, <j>o angle tends to 7r as k approaches the zone 
diagonal, and this reduces Lo s f. 
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FIG. 35. a. The theoretical result for the effective velocity 
of the quasiparticles v F — Vf/(1 + dT,'(uj)/du)). For definite- 
ness we used u) s f = 20meV, A = 1.7 and vf = 3eVA. b. Ex- 
perimental result for the effective velocity, extracted from the 
MDC dispersion M along zone diagonal. Observe the bump 
in the frequency dependence of the velocity at 70 — 80meV in 
the data and at about 3 — &u) s f in the theory. 

Since for both A(fc) and w s /(fc) there are competing 
effects which we cannot fully control, we believe that 
effective uj s f(k) and A(fc) should just be obtained from 
the fit to the photoemission data. In particular, cu s f 
can be extracted from the MDC (momentum distribu- 
tion curve) measurements of the electronic dispersion 
uj+Y!{uj) = 6fe. In Fig p5|we compare our (1+dT,' (u>) / dui) 
with the measured variation of the effective velocity 
vp{uj) of the electronic dispersion along zone diagonal ||. 
We see that the theoretical dispersion has a bump at 
lj ~ 3uj s f(kdiag)- Experimental curves look quite similar 
and show the bump at - 70 - 80meV [g|||4|. This 
yields uj s f{kdiag) ~ 25meV, i.e., almost the same as near 
hot spots. 

Note in passing that although ek+Q/vp^ 1 docs not 
vary much when fc moves along the Fermi surface, the fact 
that the Fermi velocity is very high implies that ek F +Q is 
roughly Up-\/2*0.227r/a « 0.8ey, i.e., it is comparable to 
a bandwidth. This implies that one certainly cannot ne- 
glect the curvature of the Fermi surface in the theoretical 
analysis. In other words, the Fermi-surface is very differ- 
ent from a near square which one would obtain for only 
nearest neighbor hopping. Furthermore, the fact that 
the Fermi velocity is large implies the physics at energies 
up to few hundred meV is confined to a near vicinity of 
the Fermi surface, when one can safely expand e^ to a 
linear order in k — kp. Finally, van-Hove singularities 
(which we neglected) do play some role [108, 109| but as 
eo )7r ~ 0.34£ ~ 85meV, we expect that the van-Hove sin- 
gularity softens due to fermionic incoherence and should 
not substantially affect the physics. 

We now briefly compare the experimental and theoret- 
ical forms of the fermionic spectral function and optical 
conductivity. 



1. spectral function 

We first use our form of the fermionic self-energy to 
fit the MDC data which measure the width of the pho- 
toemission peak as a function of k at a given frequency. 
In Fig. [56] we compare our results for Afc = T,"(k,fl)/vp 
with the measured || Afc vs frequency at T ~ 100K 
and temperature at Q — > J73| . For definiteness we used 
A = 1.7 and uo s f = 20meV. The slope of Afc is chiefly 
controlled by A. We see that A < 2 yields rather good 
agreement with the data on both, frequency and tem- 
perature dependence of the self-energy. On the other 
hand, the magnitude of our S" is smaller than in the 
data. To account for the values of Afc, we had to add 
a constant of about 50meV to £". The origin of this 
cons tant term is unclear. It may be the effect of impuri- 
ties 1 11C], and it also may be the effect of other scattering 
channels which we ignored ||92| . It is essential, however, 
that the functional dependence of £"(f2T) is captured by 
spin-fluctuation scattering. 
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FIG. 36. A comparison between theoretical results and the 
photoemission MDC data. For the Lorentzian lineshape of the 
MDC signal, observed in experiments, the MDC hwhm equals 
to E" /vf- Upper panel - the results for the MDC hwhm vs 
frequency at a given T. The experimental results are taken 
from bl. Lower panel - the MDC fwhm vs temperature at 
Q, — 0. The experimental results (right figure and the points 
on the left figure) are taken from J73[ . 

In Fig B7^ we present our results for the hwhm 
of the EDC (energy distribution curve) which mea- 
sures fermionic 7fe(f2) = Ak(£l)np(n) as a function of 
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frequency at a given k (Ak(Q) = (l/Tr)ImG(k,£l)). 
For Lorentzian lineshape, the EDC hwhm is given by 
£"(ft)/(l + E'(ft)/fi) §. The data are taken at T = 
115-ftT ||. Again, the theoretical slope reasonably agrees 
with the experimental one. Some discrepancy is asso- 
ciated with the fact that the experimental lineshape is 
not a Lorentzian and hence the measured hwhm is not 
exactly E"(0)/(1 + E'(fi)/Q). 
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FIG. 38. The theoretical and experimental results for the 
real and imaginary parts of optical conductivity. The data 
are from [B8l. 



FIG. 37. The theoretical result for E"(fi)/(1 + £'(fi)/fi). 
The points are the experimental hwhm of the EDC dispersion 
from 0. 



2. conductivity 



In Fig. In Fig. 38 we compare our theoretical results 



for the conductivity with the experimental data for o\ 
and (72 at different temperatures Pq] , For definiteness 
we used the same A = 1.7 and u> 8 f = 20meV as in 
the fit to the photoemission data along zone diagonals. 
We checked that the change in A affects the ratio 02/01 
at high frequencies, but does not change the functional 
forms of the conductivities. We adjusted the plasma 
frequency to obtain the agreement between dc conduc- 
tivity and Sigma" extracted from the MDC photoemis- 
sion data using vjr ~ 3eVA. This adjustment yields 
Lo p i ~ 20000cm -1 . This value is somewhat larger that 
Upi ~ 16000cm -1 obtained experimentally by integrat- 
ing o\ up to about 2 — 2.5eV piyp^. This discrepancy is 
consistent with our theoretical result that the sum rule is 



satisfied only at extremely high frequencies (see Fig. 29 ) 



We see that theoretical <Ji{u>) and 02 (w) capture the 
essential features of the measured forms of the conduc- 
tivities. In particular, the curves of 01 at different tem- 
peratures cross such that at the lowest frequencies, the 
conductivity decreases with T which at larger frequen- 
cies it increases with T (the same behavior has also been 
detected in |10jj99|| . The imaginary part of conductiv- 
ity decreases with T at any frequency, and the peak in 
02 (w) increases in magnitude and shifts to lower T with 
decreasing T |98| , |99 , 1 1 1 1 . At frequencies above 1500cm -1 
both 01 and 02 weakly depend on T and are comparable 
in the amplitude. 
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FIG. 39. The theoretical and experimental results for 
U"(u) = 47TC7iw/Wp,. The data are from J98| . 

In Fig p9 we compare experimental and theoretical re- 
sults for the imaginary part of the full particle-hole po- 
larization bubble n^(cj) = AiraiLu/u!^ Theoretically, at 
T = 0, Hff(u>) saturates at a value of about 0.2 indepen- 
dently on A and remains almost independent on frequency 
over a very wide frequency range (see Fig. E6h. We see 
that the theoretical value of H a does not change much 
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with T. Experimental data also clearly show a near sat- 
uration of 11^ at a value close to 0.2. We consider this 
agreement as a strong argument in favor of our theory. 



1.6 



1.2 - 



i °- 8 

a 



0.4 



Theory 



Experiment 



^^^ 


^ 


r^^^ 


uf // T=100K 




/ // T=200K 


- — 


f // T=295K 


frequency (cm 1 ) 


** frequency (cm -1 ) 



1500 



1500 



3000 



Lavagna, D. Khveschenko, G. Kotliar, H. von Lohneysen, 
A. Millis, M. Norman, C. Pepin, D. Pines, A. Rosh, S. 
Sachdev, Q. Si, O. Tchernyshyov, A.Tsvelik, J. Tu and J. 
Zasadzinski for useful conversations. We are also thank- 
ful to D. Basov and J. Tu for sharing their unpublished 
results with us and to M. Norman for bringing ref. B3| 
to our attention. The research was supported by NSF 
DMR-9979749 (Ar. A and A. Ch.), and by the Ames 
Laboratory, operated for the U.S. Department of Energy 
by Iowa State University under contract No. W-7405- 
Eng-82 (J.S). This work was partly done while A.C. was 
on a sabbatical leave at the Rutgers University and the 
NEC Research Institute. The hospitality of both places 
is acknowledged with thanks. 



FIG. 40. The theoretical and experimental results for 
1/r* = Lutji/a2- The data are from [B8|. 

The agreement between our theory and the experi- 
ment is, however, not a perfect one. In Fig. En we show 
theoretical and experimental results for 1/r* = Luai/a^- 
The advantage of comparing 1/r* is that this quantity 
does not require one to know what the plasma frequency 
is. We see that while both experimental and theoretical 
curves are linear in frequency, the slopes are off roughly 
by a factor of 3. This discrepancy is related to the fact 
that in our theory, at high enough frequencies, Pi' a {uj) 
is roughly 3 times larger than II" (w) (see Fig. |2fj) and 
hence <Jil &\ ~ 3, whole experimentally <j\ and 02 nearly 
coincide at high frequencies. Also, at zero frequency, the 
theoretical value of 1/t*(u> = 0) at, say, T = 300K is 
larger than in the data @,Q . From Fig |Zj] the ratio of 
theoretical and experimental l/r*(f2 = 0) is roughly 1.5. 

The discrepancy in 1/r* is important for the theory as 
it indicates that either the averaging over the Fermi sur- 
face and vertex correction within a bubble, or the RPA- 
type corrections to conductivity pq,B3] do play some role. 
Still, however, Figs. Bq and B3 indicate that the general 
trends of the behavior of the conductivities, at least near 
optimal doping, are reasonably well captured by the spin- 
fluctuation theory. 

In this paper we collected all what we know at the 
moment about the normal state properties of the spin- 
fcrmion model. The next step is to consider what hap- 
pens below the pairing instability. We plan to present 
the detailed account of these results in the near future. 
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XII. APPENDIX A 



In this Appendix, we evaluate several integrals pre- 



sented in Sec. Ill 



A. direct perturbation theory 

We first compute the second order fermionic self- 
energy and the vertex correction in a direct perturbation 
theory for the spin-fermion model. "Direct" means that 
the calculations are performed with the bare value of the 
spin susceptibility. 

We begin with the perturbativc fermionic self-energy 
E(k,n m ) = (fi m - ek+Q) l(k hs ,0), where l(k hs ,0) is 
given by (§) 

/(kfc.,0) - 7^3 J d 2 cidLo n 



£" 2 + q 2 + (uj/v s r 

x l — — 1 (150) 

[iuj m - v F q x y 

Here q = q Q, and ir— axis is chosen along vp at k + 
Q. Rescaling the momentum and frequency q x = £ x, 
q y — £ _1 y, and uj — £ v s z, transforming to the spherical 
coordinates and integrating over r = ^Jx 2 + y 2 + z 2 , we 
obtain (rj — vp/v s ) 

tt/2 2tt 

«*.,„)- * / * / — 

2-kt] J sine J (irj-i 





tt/2 



= -\o 



(ir) 1 cot ( 



sin 6 cos 6d9 



1 3/2 



{ [l + (?7 2 -l)sin 2 #]" /z 1 + V 

Similar analysis works also for the vertex correction 



, (151) 



We have from (14) 



fJ 



1 



Ag g f ,2 - , 

8tt 3 7 £-2 + q2 + (tJ / Ws) 2 



1 



1 



iu - {v x q x + v y q y ) iuj + v x q x - v y q y 



(152) 
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Rescaling the momentum and frequency in the same way 
as before, we obtain from (|152|) 



susceptibility is justified when a 3> 1 (see also below). 
The integration over y is straightforward and yields 



Ag A vp 

g 6ir 2 v s 



dxdydz 



1 



1 + x 2 + y 2 + z 2 [iz — ay) 2 — j3 2 x 






(153) 



l(k hs ,0) 



where a = v y /v s and (3 = v x /v s (a 2 + (i 2 = r/ 2 ). Intro- 
ducing spherical coordinated and integrating over r we 
obtain 



A5 
9 



Xv F f n/2 d0 



6ttv 

\vf r /2 _ 

SVsJo [ gin 2 

A vp 1 



cos( 



(i tan 9 — a sin < 
d9 cos 9 sin 9 



(i 2 cos 2 cj) 



a 2 cos 2 9] ^sin 2 9 + {a 2 + /3 2 ) cos 2 9 
a (Vl-a 2 - 0) 



3 v s (3^/T^a 2 v/(/3 2 + a 2 ) (1 - a 2 ) - (i 



(154) 



At u s — > oo, i.e., for a purely static bare spin suscepti- 
bility, Ag/g becomes 



Aff _ Xvf_ . x v^ 
3 3u x u„ 



(155) 



For «* = «„, we have ^ = ^ i n (i + ^2) « 0.415A. 
In the opposite limit v s <C vf, we obtain from ( 154 ) 



A5 
5 



A 



tan 



3 v x v y 



(156) 



B. renormalized perturbation theory 

We next repeat these computations using the relax- 
ational form of the spin susceptibility. We first re- 
evaluate I(k/j S ,0). From Eq. (E3) we have 



I(k hs ,0) = 3ge 



1 



d qduj m 



(iUr, 



VFq x ) 2 



(157) 



The key difference with the direct perturbation theory 
is that now the uj 2 term in the susceptibility becomes 
subleading at low frequencies. Let's first neglect this 
term. Introducing, as before, q x £, — x, q y t; — y, and 
also uj m — Lu s ft, we then rewrite Eq. ( 157 ) as 



/(k hfl ,0) 



3V X Vy 

2tt 2 Nv 2 f 



(x — it/a) 2 



dxdydt 



1 



l+x 2 +y 2 + \t\ 



where a — (vf£, 1 /w s /) = (NX/3) (v F /v x v v ). One can 
easily make sure that the neglect of the oj 2 term in the 



3v x v y 

2irNv 2 F 



dxdt 



(1 



|t|) 1/2 



(x — it/a) 2 



(159) 



The last term in (159) is a double pole. The integral 
over x then does not vanish only because the first term 
contains a branch cut. Deforming the contour of integra- 
tion over x to include the branch cut along imaginary x, 
and integrating over the branches of the branch cut, we 
obtain 



/(kfc.,0) 



l2V X Vy 

ttNv 2 f 



;l/a 



dt- 



1 



1 + * 



(160) 



Evaluating the integral with logarithmical accuracy, we 
obtain 



IregQths,0) = - 



12VxVy 

r- r " log A 

nNv 2 F 6 



(161) 



This is the result we cited in the Eq. ([24|). 

With little more efforts, one can evaluate l(kh s ,0) 
for arbitrary a and v s /vf which are two parameters for 
7(k^ s ,0). To avoid presenting very long formulas, we 
only consider the two limits v s — > oo (i.e., a purely static 
bare spin susceptibility), and vf = v s . In the first case, 
performing the same computations as above but keeping 
a arbitrary, we obtain after some algebra 



l(k hs ,0) = -- 






ydydx 



"" Jo Jo (1 +x + ay)(x + y 2 y/ 2 
The straightforward integration yields 

1 + flogf 



IregV^hs-i U) 



-A 



1 



(162) 



(163) 



This is the result we cited in the Eq. fl25[). 
Similar considerations for v s — vf yield 



iQth.,0) 



4A 



( 7T— a— v 4 — a 2 cot l ■ I 



(164) 



At a — > 0, l(khs,0) — ► —A/2 in agreement with (151). 

We next show how to compute the full (regular plus 
anomalous) second-order £(k/j S ,J7) at arbitrary fi. We 
have 



(158) E(k hs ,fi r . 



i 



d qduj Tl 
(2tt) 3 



1 + (qO 2 



/Wet 



vW 



if2. 



VFq x 



(165) 



Introducing the same variables x,y,t as before, and also 
ft = suj s f and using the spectral representation of the 
Green's functions we obtain 



ImE(fi) 



Aw 



»S 



Ira I dt 

o 



dy 



y" + l-it + ((s + t) 2 -r) 2 t 2 )/a 2 
The integration is straightforward, and the result is 
\ijj s fa 



m^±S±^^i (166) 



ImYt(VL) = Ira- , 

yr^f vh + vt 

where t\ and ti are the roots of the quadratic equation 



(l+it)a 2 = (s + tf -rft 2 







(167) 



For p urely static bare susceptibility [v s = 00, or rj = 0), 
Eqn (166) is simplified to 



£(fe/j, s ,fi) = 2\tu s fa In 



iV^n - 1 + y/Kn + 1 



iVK n -l+Aa + y/K(n) + 1-A n 

(168) 



where 



Kl 



1 



1 



. n 



Ac 



2z 

a 2 w s / 



(169) 



(we used Kramers-Kronig relation to obtain i?e£) 

We next evaluate the lowest order correction to the 
spin-fcrmion vertex at the bosonic momentum Q. The 
expression for Ag in terms of fcrmionic Green's functions 
is given by Eq. (Q. Neglecting the subleading uJ^ in 
the spin propagator, introducing x — q x £, y — q y £, and 
t = uj/uj s f, and expanding the quasiparticle energies to 
linear order in deviations from hot spots we obtain from 

(0) 

"7" = 2^N J dxdydt l + x 2 +y 2 + \t\ 

x j (170) 

(v x x) 2 - (v y y - iv F t/a) 

we remind that a = (NX/3) (v F /v x v y ) 3> 1. Let's per- 
form the integration over y first. As both poles in the 
last term in ( |170| ) are in the same half-plane, it is con- 
venient to close the integration contour over a different 
half-plane where only the spin susceptibility has a pole. 
Evaluating the integral we obtain, 



A3 
9 



v x v v 



dt 



dx 



*N J q 7-oc y/l+X* + \t\ 

(v F x) 2 +v 2 y (l + t)+v 2 F (t/af 
t 



+ 2v y V F - y/l + X 2 + \t\ 



(171) 



To a logarithmical accuracy, we can neglect 1/a terms 
in the integrand of (171) and set the upper limit of fre- 
quency integration to be t max ~ a 2 . Introducing next 



x = (1 + i) l ' 2 p, and performing the integration over p 
first, we obtain from (171) 
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2iV 



Q 



dt 

T+t 



(172) 



where 



qC- 

Vy 



4 V x Vy 



dp 



1 



4 

— tan 

7T 



y/l+p 2 P 2 + (v y /v F ) 2 
x v x _ 4</) 



(173) 



and 0o is the angle between the directions of velocities 
at ho t spots separated by Q. The integration over t in 
( 172 ) is elementary, and we obtain 



g N " \ i\, 



log A 



(174) 



This is the result we quoted in the text. For A — ► oo, 
the logarithmic dependence on A transforms into the log- 
arithmic dependence on the external frequency: log A — > 
(l/2)logw m0 a:/M where uj max ~ v F /g. 

It is instructive to observe that the same result can also 
be obtained by noticing that the u> terms in the fermionic 
propagators are much smaller than the y/ijj term from the 
spin susceptibility, Eq. (17C). These u> terms can then 
be neglected, and the vertex correction can be straight- 
forwardly simplified to 



A5 
9 



2V X Vy 

TT Z N 



where 



dxdy 1 

l+i 2 +y 2 (xv x +yv y -i0)(xv x -yvy+i0) 



(175) 



(176) 



and compared to our previous notations, we introduced 
x = x/t 1 / 2 , y — y/t 1 / 2 . Naively, one might expect that 



the 2D integral in (176) is determined by the two poles 
at vanishing x and y. However, the contribution to the 
integral from small x,y diverges as log , and has to be 
regul arized. The extra 1 + x 2 + y 2 in the denominator of 
(176) provides such a regularization. We checked that for 
v x = v y , the regularization is irrelevant, and the result 
for Ag/g comes from the two poles at vanishing x, y. 
For other ratios of velocities, the regularization yields an 
extra contribution to the integral from the range where 
x, y = 0(1). A straightforward calculation then yields 
/ = (2ir/(v x v y ))t&Ti~ 1 v x /vy 1 and hence the same Ag/g 
as in (174). Physically, this result implies that except for 
the case when v x — v y , a part of the vertex correction 
comes from small momenta (in reality, of order u)/v F ), 
and a part comes from much larger momenta, of order 

(Lo/(L0 sf e)) 1/2 ^>u/v F . 

Above we computed Ag/g with the logarithmical ac- 
curacy. One can easily check, however, that the 3D 
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integration over momentum and frequency in (H) is 
convergent, i.e., the subleading, non-logarithmical term 
can also be determined within the low-energy model. 
This calculation is more convenient to perform in po- 
lar coordinates, by expressing et = VFqcos(cj) + </>o/2), 
£/c+Q = VFqcos(4> — 4>o/2). The integration over q is 
straightforward, and performing it we obtain 



Ag 2 1 sin <fi 



fJ 



~2 N 



Re 



(21oga + logsin<?(>/2) 



log[sin(0/2)] 

COS (f> + COS (j)Q 



(177) 



(NX/3)v f /(v x3l 



where, we remind, a = vp^ 1 /uj s f 

The log a term in (177) is indeed the same as in (174) 
The second piece is the regular contribution to Ag/g. 
Evaluating this piece numerically we find that for </>o = 
7r/2 i.e., v x = v y , which is close to the actual situation 
in cuprates, Ag/g = (l/N)(loga — 0.93). This result 
implies that for moderate A, vertex correction is rather 
small numerically: Ag/g w 0.09 for e.g., A = 1, and 
Ag/g w 0.18 for A = 2. 



vf scales with N, each running fermionic line acquires 
a factor 1/N. Similarly, since ~g scales with N, a dia- 
gram with n vertices acquires a prefactor N n / 2 (recall 
that g ex g 2 ). Finally, if a given diagram to the thermo- 
dynamic potential has to closed fermion loops, a summa- 
tion over fermionic flavors yields an additional prefactor 
N m results. As a result, a diagram with n-vertices and 
to closed fermionic loops behaves as 



D(n,m)<xN m - n / 2 , 



(180) 



We used the fact that the number of internal fermionic 
lines in the diagram is n. 





XIII. APPENDIX B 



In this Appendix we discuss the general structure of 
the 1/N expansion. Suppose for defmiteness that we keep 
the number of hot spots finite and extend the theory to a 
large number of fermionic flavors. As usual, the extension 
of the theory to large N implies an appropriate rescaling 
of the coupling constant g. In our case, however, we also 
have to rescale the Fermi velocity in the same way as 
g. Indeed, consider the lowest order fermionic loop (the 
particle- hole bubble). We found above that it is propor- 
tional to ~g~N/vp. The power of the Fermi velocity in the 
denominator accounts for the number of fermion propa- 
gators in the diagram. Evidently we would like this dia- 
gram to be of order one. Then we will have a "zero-order" 
theory with no factor of N. Usually this is achieved by 
assuming that the coupling constant is inversely propor- 
tional to N. However if we do this extension, we find 
that the anomalous fermionic self energy, which is also a 
part of the "zero-order" theory scales as 1/N as it is pro- 
portional to g/vp. (recall that the fermionic self-energy 
contains only one fermionic line and no summation over 
the flavor index). A proper extension which makes both 
particle-hole bubble and anomalous fermionic self-energy 
of order 1 is 



and 



g^gN 



vp — > vfN. 



(178) 



(179) 



We now consider the general structure of the Feynman di- 
agrams not included into the "zero-order" theory. Since 





FIG. 41. A diagram for the thermodynamic potential to 
zero order in 1/N. 

The thermodynamic potential for the "zero-order" the- 
ory is given by the diagram shown in Fig.Bll In this di- 
agram, to = n/2 and hence D = 0(1). The fermionic 
and bosonic self energies are obtained by cutting one 
fermionic (or bosonic) line, and are also of order 1 as 
they indeed should be in the "zero-order" theory. Note, 
however, that this independence of N is the result of the 
interplay between g and vp as if one cuts a fermionic line 
one looses one factor N since one closed fermionic loop 
disappears. This is compensated by the fact that there 
is one running fermionic line less in the diagram. 



^' 



-% f > 



\ _^ 



w 
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FIG. 42. The diagrams for the thermodynamic potential 
to order 1/N. 

The diagrams which constitute the expansion in j? 
are then obtained by reducing the number of internal 
fermionic loops at a given number of g. The leading cor- 
rections are shown in Fig. [l2| These diagrams have 
m = n/2 — 1 and hence are of order 1/N. Cutting 
fermionic or bosonic lines yields self energy and vertex 
corrections to order 1/N. Higher order diagrams are ob- 
tained by further reducing the number of closed fermionic 
loops. 



XIV. APPENDIX C 



x / dfc-i 



x / dci 
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1 




b + + — " 
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-i^b+- 
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- — + q- 
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C- ■ u> 



»£ C + + -T + 9+ 



(182) 



The integration over db + and dc+ is straightforward and 
yields 



n 2 (q,fi) 



2g 2 N 1 
1 



T 2 J2 db - dc 



1 \LO-Uj'\ + {b 2 _+C 2 _)/2 

[sign (a;' + CI) — signw] [signw' — sign(w + Cl)] 

[b __ q _ + li y^±n ][c _ +q ^ 



; 0J — 0j'-\-i1'\ 



(183) 



In this Appendix we evaluate the polarization bubble 
with the vertex correction. We first compute this bub- 
ble explicitly. Then we redo computations in a different 
way: we first evaluate the effective four-fermion vertex 
and then contract two of its opposite external legs (see 
Fig. Hq) and convolute the vertex with the spin propaga- 
tor. This last computation allows us to directly compute 
our theory with the </> 4 theory of critical behavior for the 
dynamical exponent z = 2. 



A. direct computation 

The diagram which we need to compute is presented 
in Fig. [lqc In the analytical form we have for £ = oo 



2g 2 N 



n 2 (q,0) = -f^T 2 £ d 2 kd 2 k 



1 



1 



7|w 



1 



€k - iuj e k + q +Q - i(u> + Cl) 

1 1 

^k'+q - i(w' + Cl) e/c'+Q - iw' 



Ux™x 



where, as before e^ = v x k x + v y k y , €k+Q 
Vyky, and k measures a deviation from a hot spot. The 
combinatoric factor 27V comes from spin summation and 
summation over hot spots, and 7 = Ng/^AnVxVy). 



Consider for simplicity the case v x 



Introducing 



UrSt Kx ~r rZy — ri-\- , rZy f^x — ^— : <^ x ' 11 — + ' y X — 

k'_, q+ = q x + q y , q- = q y -qx, and then (k+ + k' + ) = 2b + , 
k+ - k' + = b-, (k- + k'_) = 2c + , k- - k'_ = c_, and 
substituting into (181) we obtain 



n 2 ( q ,o) = - 



2g 2 N 1 



T 2 J2 db - dc 



1 



7|w-W| + (^.+^.)/2 



Shifting the momenta 6_ — g_ = b-, c_ 
rewrite (183) as 



(k-k') 2 



(181) 



= c_ we 
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(2tt) 4 iv x 
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rjnZ 



UJ.LO' 



db-dc- 



1 \ u - u '\ + ((b-+q-) 2 + (c-q + ) 2 )/2 
[sign(u/ + Cl) — signw] [signa/ — sign(w + Cl)] 
\b- + i^^+ill [c_ - i^Lji+E 



(184) 



Next step is to observe that typical 6_ and c_ are of 
order of a typical frequency, and hence can be safely ne- 
glected in the bosonic propagator, in which typical mo- 
menta scale as a square root of a typical frequency. The 
remaining momentum integration then proceeds straight- 
forwardly, and we obtain 

n 2 (q,o) = -^4 tY ^(* + ")°W-*) 

2VH ' ' (2tt) 4 4w4 ^ 7p|+g 2 

xT^(sign(w' + Cl) - sign (a/ + &)) 

U)' 

x(signu/-sign(u/ + w + 0))] (185) 

where lu = lu — lu' . At T ^ 0, the summation over lu' 
is elementary and yields 4(|a)| — Cl). Substituting this 



result and the expression for 7 into (185), evaluating the 
integral over u) up to an upper cutoff A, and re-analyzing 
the computational steps for v x 7^ v y (this adds a factor 
Q(v) — (4/7r) ta,n~ 1 (v x /vy)) we finally obtain 



n 2 (q,fi)^%l(q 2 +7|M) 



A 



log- 



Aj 



iM 



log- 



7 |0| 



-(A-2\Cl\). (186) 



The constant contribution (= A) has to be neglected to 
avoid double counting as is comes from high energies and 
has already been absorbed into a bare Xo(q, Cl). The 
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logarithmical terms, on the other hand, come from low 
energies and have to be kept. 

We see that to a logarithmical accuracy q 2 and 7|w| are 
renormalized in exactly the same way, i.e., the diffusion 
coefficient does not diverge as the system approaches a 
QCP. It does acquire a nonsingular 1/iV correction which 
is beyond the scope of our analysis. We however checked 
that the logarithmical divergence at q = is fictitious 
and is eliminated when one included the self-energy of 
intermediate fermions. 

Let's also check whether usloguj term at T = is col- 
laborated by TlogT term at a finite T, i.e., whether 
vertex corrections enforce w/T scaling in the spin sus- 
ceptibility. To address this issue we evaluate II2 (k, CI) at 
finite T and CI = q = 0. The summation over Matsubara 



frequencies in (185) is straightforward. Performing the 
summation over u' and substituting the explicit form of 
7 we obtain 



n a ( q ,n) 



on/ \ A/(2ttT) 



A/ 



P =i 



P 



¥- A + 2 - T ^4r 



(187) 



This result implies that finite T accounts for the shift of 
the bosonic frequency by 2nT, i.e., it just shifts a bosonic 
CI = 2irTn from n = to n = 1. Obviously then, in real 
frequencies, II2 (q, CI) does not contain a singular T log T 
term, i.e., there is no Cl/T scaling in the spin susceptibil- 
ity at the QCP. 



B. four-boson vertex 

We next explicitly compute the four-boson vertex in 
Fig.|l6|. In analytical form it is given by 



Ng 4 f dui'd 2 q' 



1 



(27r) 2 J iuj' - e q i i(Cl + ui + u>) - tQ +q > +qi+q 
1 



i(w' + Ul + UJ 2 ) -£qi+ qi +q 2 
1 

i(w' + ijJ + uj 2 ) ~ e Q+q/+q+q2 



(188) 



Linearizing, as before, the fermionic dispersion near hot 
spots as e' q = v x q' x + v y q' y , e q > +Q = -v x q' x + v y q' yl and 
performing the momentum integration, we obtain 

Ng A f .sign(w') - sign(cj' + U\ + oj 2 ) 
duj - 



16irv x v y 



e 9i+«2 -i(wi+ua) 
sign(o;' + u>i + u) — sign(w' + u + ll> 2 ) 

e Q+qi-l2 ~«(Wl ~W 2 ) 

The integration over u>' is also straightforward and yields 



b = 



Ng 4 \u> 2 + lj\ + \uj 2 — lu\ — \oji + w\ — \u>\ — U)\ 



8irv x v y [i(u>i +uj 2 ) - e gi+93 ][i(wi - w 2 ) - tQ+ qi - q2 \ 

(189) 

This result is quoted in the main text. 

C. First-order bosonic self-energy 

We next explicitly compute the first oder bosonic self- 
energy by contracting the external legs in the four-boson 
vertex and convoluting it with the spin susceptibility. 
The corresponding diagrams are presented in Fig. |l6| 
b and c. 

The first diagram is obtained by contracting the adja- 
cent external legs of the four-fermion vertex. In this sit- 
uation u>\ = lu 2 and qi = q2- Substituting this into Eq. 



(189) and then into the expression for the bosonic self- 
energy, we find that the self-energy vanishes due to still 
presence of double poles. This is consistent with our ear- 
lier result that the inclusion of the fermionic self-energy 
does not affect the form of the bosonic propagator. 

The second diagram is obtained by contracting the ex- 
ternal legs which are not adjacent to each other. In this 
situation, u) = q = 0. Let's choose Wi to be an external 
frequency (which we label CI for consistency of notations) 
and at first set the external q± — (this corresponds to 
external momentum q = qi + Q = Q = (ir, 71")). The 
bosonic self-energy is then given by 



n 2 (Q,fi) 

(27r) 4 u a; u a J [i(Cl + u 2 ) - e q2 ][i(Cl - w 2 ) - 
duj 2 d 2 q 2 



b(Cl,uj 2l uj = 0)x(q 2 ,w 2 

bai-ini 



duj 2 d 2 q 2 
' (2tt) 3 



-qi J 



7l^2|+q 2 



(190) 



(the combinatoric factor 2 accounts for two possibilities 
to choose the external frequency, the factor —1 comes 
from the summation over spin projections, which we re- 
installed here). One can easily make sure that to a loga- 
rithmical accuracy, we can neglect CI in the denominator 
and set it as a lower cutoff in the integration over lu 2 . 
Subtracting the contribution from high energies we then 
obtain 



n 2 (Q,n) 



16n 4 v x v y 
duj 2 d 2 q 2 

7 1^2 1 + q 2 



fil 



{iuj 2 - e q2 ){iuj 2 - e Q+q2 ) 



This integrand is exactly the same as for Ag/g in the 
Appendix A. Using the results from Appendix A we find 



n 2 (Q,n) = 7|o|%/%ogfi| 



(191) 
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We next evaluate the correction to the polarization op- 
erator when both ft and qi = q — Q are nonzero. For 
simplicity we choose v x = v y — v/\/2, and also set q = q x . 
We will restore rotational symmetry at the end of compu- 
tations. Rotating the coordinates such that e 92 = v F q2 X 
and tQ+q 2 — VFq_2 V we can express n 2 (q, fl) from ( |190| ) 
as 



n 2 (q,n) 



-2^ff 4 X0 



M-ini 



i(Sl + ccJ 2 ) - v F q 2x 
duj 2 d 2 q2 



i(fi - uj 2 ) + v F q 2y - v F q y j\lu 2 \ + q\ x + gf 



(192) 



As before, the dependence on 17 in the denominator can 
be eliminated by setting tt as the lower limit on the in- 
tegration over lj 2 . Furthermore, as typical uj 2 scale as 
q\ and are obviously smaller than q 2 , the iu 2 term in 
the fermionic propagators can be reduced to i^signf^)- 
With these simplifications, we obtain 



n a (q,n) 



4iVg 4 xo 

(27T)*t£ 



Re / dio 



-*(L 



v F q 2x 



d 2 q 2 



-i0+ + v F q 2y - v F q y 7|w 2 | + q 2x + 



'4 
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Now one can easily evaluate first the integral over q 2x 
and then over q 2y and obtain 



n 2 (q,o) 



Ng 4 X o 
4tt 2 w 4 



du> 



to 2 



mi 



7W2 



(194) 



Evaluating the frequency integral, subtracting the con- 
tribution from high frequencies, restoring the rotational 
symmetry, and expressing 7 in terms of g and v F we 
finally obtain 



n 2 (q,fi) 



Q(y) 

N 



(q 2 + 7 |fi|)| log( 7 m| + q 2 )| (195) 



To a logarithmical accuracy, this is indeed the same result 
as Eqn. (p6J). 



XV. APPENDIX D 

In this Appendix, we evaluate the fermionic self-energy 
at finite temperatures. We first consider the limit N — v 
00 at finite A, and then discuss the modification of the 
results at finite N. 



A. N — * 00 limit 

The sum we need to evaluate is given by Eq. (|77|). It 
reads 



Y,(k hs ,£l m ) = inTX ^ 

n 

Using Poisson summation formula 

n °° 

^ 00 



signw„ 






(196) 



xp 



00 ^00 

(F(x) + F(-x)) cos ^dx (197) 



P =i 



to separate T = and finite T contributions to the self- 
energy, we obtain after integrating over x in both terms 
in (lL97b 



E(fc/j S ,O m )=2iA j 



a 



in m i 



-f-P(fi m ,T)signO r . 



(198) 



The T dependent term P(f2 m ,T) is expressed in terms 
of Fresnel integrals S(x) and C(x) as 
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The Fresnel integrals are defined as 



S(Vx) = 



1 



2tt Jo \/i 



sint „, . — 1 

dt C{y/x) = 



x cost 



2ir J y/i 



■dt 



(200) 



At low T -C tj a f, one can expand Fresnel integrals in 
powers of 1/ '\/x. In this limit, Eq (19S) reduces to 



P(fi m ,T) = ^(l + 2 



>>t 



in, 



1/2N 



8 W 



«/ 



Ifir. 



0(T 6 ) 



(201) 



where /(0) = 0, /(oo) = 1. At O m <C w s /, one re- 
covers a conventional Fermi liquid result P(Q m ,T) rj 
(7tT) 2 /(4w s/ ). At f2 m > w s /, but still T < w 8 / (i.e., 
for large to), the leading functional dependence is still 
T 2 , but the prefactor is reduced by 1/3. Also notice that 
the subleading, T A corrections become comparable with 
T 2 terms starting from relatively small T ~ 0.5 w s /. 

In the opposite limit T 3> 0J s f, the contributions to 
the sum ( |196] ) from all n 7^ to can be simplified by ap- 
proximating 
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1 



-1/2 



J 'f 



1/2 



(202) 



U> s f / \JWn-mL 

With this approximation, one can sum over n explicitly: 

sign(2n + 1) 



\ /"OO 1 

Z dyn (n)^-Im / /i(n/2T, 

2 y_oo J\-%s- 



E(fi m ) = *^ T^ 1 / 2 ^ 



w/2T) (211) 



At f2 ^> w 8 /| typical cj ^> w s /, and Edyn(^) can be ap- 
proximated by 



1/2 



11/2 



C?) "(<G)-4 



^dyn(i') 



OO 7 

^=f x {n/2T,w/2T) 



(203) 



where £ is a Zeta function. The absence of the correlation 
length in this expression is a consequence of the fact that 
the product XJuJJJ = 0.5\/<^ is independent of £, One 
indeed would expect this independence in the quantum- 
critical regime. 

In real frequencies, the imaginary part of the self- 
energy at finite T can be obtained using the spectral 
representation . We have at k = k/j S and v s — > oo, 



£"(fi) = 



\v F £; 



2tt 2 
x Im 



Im 



dujcPk 



i + (kty 



W.f 



v F k x - £«(Q • 



-f(n/2T,u>/2T) (204) 



where 



/(a;, y) = tanh(x — y) + 1/ tanh(y). 



(205) 



At JV — > oo and finite £, the momentum integration 
can be factorized and wc obtain: 



E"(0) = £/m jf 



cko 



f(n/2T,cj/2T) (206) 



Note that the presence of T, R in the integrand is irrel- 
evant as when momentum integration is factorized, the 
integral over k x just gives a constant density of states, 
independent on Y, R . It is convenient to separate the 
contributions to E"(f2) from static and dynamical spin 
fluctuations. This can be done in a standard way by 
replacing f(x,y) by 



V 



(207) 



where 

h (as, y) = tanh(.x - y) + 1/ tanh(y) - 1/y, (208) 
Substituting the form of f(x,y) into ( |206| ) we obtain 

s"(n) = s;' t (n) + s d3/n (o), (209) 

where 

S' S ' 4 (Q) = TAtt (210) 

and 



4%/2 J_oo V^~ 
This result can be rewritten in the scaling form 



/ rp- 

^dynity = I — 

where D(x) is given by 
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(213) 



1 



e»- 1 



11/1 1 

y 2 \ e v+ x + 1 e«- x + 1 



(214) 



= V 1 
-1.516 



In the two limits, -D(x) behaves as D(x 3> 1) 
2.58 - 0.822/a; 3 / 2 + ..., and D(x < 1) = 
0.105a; 2 + .... 

The real part of Sd yra (^) is obtained from ( |213| ) by 
Kramers-Kronig transformation. It contains the same 
y/T dependence as ( 213 ) but a different scaling function 
of fi/T. 



B. finite JV 

We recall that the need to study a finite JV is related 
to the fact that at a nonzero temperature, some of 1/JV 
terms also scale with £ and therefore become relevant 
for physical N — 8. Here we single out and explicitly 
evaluate there terms. We, however, still will be neglecting 
regular 1/JV contributions. 

Our point of departure is Eqn ( |204|) . At finite JV, 
the momentum integration in ( |204| ) cannot be factorized. 
Let's use (207) and consider the static and dynamical 
contributions to £"(Q) separately. This separation re- 
quires care as to single out the truly static contribution 
(the one from a boson with zero Matsubara frequency) we 
not only have to replace f{x, y) by 1/y but also set u> = 
in the fermionic propagator. The static contribution is 
given by 



SM = ^€ Hm 



2tt 2 
x Im 



dujd 2 k 



2T 



1 



Ul.f 



v F k x - E R (Cl) 



(215) 



The integration over u) is straightforward. Also, since 
ImYi R > 0, we can integrate over k x by close the integra- 
tion contour over the lower half plane. Introducing then 
£ = E-^fi)/^- 1 , we obtain 
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S^(0)=rA/m/ 

J oo 



rf.r 
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2TAi> 



log 



Va; 2 + 1 + £ Va; 2 + 1 

e + vT+I 



Vi + s 2 



(216) 



As we discuss in the main text, at T <C Nu> s f, the 
fermionic self-energy in the integrand can be approxi- 
mated by its T — value, Eqn (pq). The latter, however, 
can be neglected only iff2<§;lorf2^;J7 = N 2 ui s f/9. 
For infinite N, this scale is above Hi and is irrelevant for 
our study. For finite N, however, Cl scales with w s / and 
definitely become smaller that w if £ — * oo at a finite N. 
Obviously, in this situation, E cannot be neglected, and 
the N = oo result for £^, Eq. ( |210| ) is invalid. 

In the quantum critical regime, E 2 = ifl/fl. Substi- 
tuting this into ( |216| ) we obtain the result cited in (p2|). 

The dynamical contribution to £"(Q) consists of two 
parts E dyn {il) = E^ n (Q) + S^ n (fi). The first contribu- 
tion accounts for the difference between T, R (Q — w) and 
E H (fi): 



^dyn(^) 



x Jm 



2tt 2 
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2T 



1 + (k£) s 
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v F k x - H R (VL - u) v F k x - Y, R {fl) 



(217) 



This integral is nonsingular when £ — > oo. We can then 
use a conventional TV — > oo approach and factorize the 
momentum integration. Performing it, we immediately 
find that E^ (fi) vanishes to leading order in 1/N. 

The s econd contribution to the dynamical part is given 
by (204) with /i instead of /. As fx(x, y) is non-singular 
at y — » (i.e., at oj — > 0), the leading (in 1/iV) piece in 
E 2 , can again be obtained by factorizing the momentum 
integration. This piece is then the same as in ( |213| ). 

Finally, we found that the momentum integration in 
(204) can be performed explicitly at arbitrary N by, e.g., 
transforming to polar coordinates. Furthermore, by ex- 
plicitly evaluating ImG(k, f2) and pulling the Im symbol 
for bosonic propagator out of frequency integral, we can 
obtain the full expression for E(0,T), not only its imag- 
inary part. Performing the momentum integration, we 
found 



A f°° 
E(0,T) = — : / duj f(n/2T,ui/2T) 

(Sii(n,uj) — Sa(Q, ui)) sign(f2 — lo) 



(218) 



where 

S R>A (Q,u;) 



Zl A + Iln 



x log 
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T, 2 
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(219) 



n, 



1 — iuj /ui s f and E 



2 
R,A 



^R,A 



(Q 



and Hr = II_r(w) 

^/("Ft 1 )- 

This expression is of limited use as it contains regu- 
lar \/N terms which are beyond our accuracy as we ne- 
glected 1/N vertex corrections. Nevertheless, it is useful 
for the understanding how small are regular 1/JV terms. 
The E" from Eq. (p>18|) is plotted in Fig. M. 
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